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FOREWORD 


The  process  of  learning  to  use  a  complex  computer  program  is  most 
often  frustrating.  Regardless  of  the  care  taken  to  dociuient  and  explain 
the  various  features  of  the  program,  a  new  user  is  easily  confused  by 
unfamiliar  terminology  and  descriptions  of  the  methods  employed. 

Usually,  the  developers  have  become  so  conversant  with  the  problems 
addressed  by  the  program  and  the  way  that  they  have  interpreted  and 
solved  them,  that  it  is  difficult  to  communicate  with  the  uninitiated. 
They  simply  forget  how  far  they  have  come  from  the  beginning.  Often  the 
appropriate  beginning  point  can  be  found  and  quick  learning  obtained 
through  dialogue  between  the  developer  or  an  experienced  user  and  the 
novice.  Unfortunately,  the  written  page  does  not  provide  for  such 
feedback-  One  should  keep  in  mind  that  it  is  at  least  as  difficult  to 
write  a  manual  that  adequately  describes  a  program  as  it  is  to  learn  to 
use  the  program. 
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PREFACE 


The  SEADYN  cable  dyiiainics  computer  program  was  ilevelopcct  over  a 
6-yt?iu  period  primarily  by  the  Naval  Civil  Eiigi iioeriiig  Lal'oratory  (NCdU.) 
under  Uie  sponsorship  of  the  Naval  Facilities  Engineering  Command 
(NAVFAC) .  SKAUYN  is  the  largest  of  several  cal)l<?  analysis  programs 
developed  under  the  Large  Displacement  Cable  Dynamics  project.  I'liest' 
programs  have  been  validated  by  laboratory  and  at-sea  experiments  and 
allow  for  the  confident  analysis  of  arbitrarily  configured  cable 
structures  subject  to  a  wide  variety  of  environmental  and  system  loads. 

This  report  is  one  in  a  scries  of  five  used  to  document  the  SEADYN 
computer  program.  This  set  of  reports  consists  of; 

(1)  SEADYN  User's  Manual:  Describes  the  program  input  format 
and  externai  file  re<juiremenLs,  and  briefly  discusses  use  of 
the  program.  Example  inputs  are  presented. 

(2)  SEADYN  Theoretical  Models:  Describes  the  finite-element 
formulation,  implementation  of  particular  submodels  (slrunmiing 
model,  linearizations  in  the  frequency  domain  solution,  etc.), 
and  numerical  solution  techniques  (Ref  1). 

(3)  SEADYN  Programmer's  Reference:  Describes  SEADYN  coding 
structure,  logic,  memory  usage,  and  required  system  routines 
(Ref  2). 

(4)  SEADYN  Test  Cases  Report:  Includes  detailed  input  and  corre¬ 
sponding  output  for  use  in  confirming  the  program's  operation 
(Ref  3). 

(5)  SEAFLOT  User's  Manual:  Describes  input  format  for  the 
graphics  program  that  complements  SEADYN' s  tabular  output; 
requires  the  CDC-Cybernev  UNIPLOT  graphics  support  package 
(Ref  4). 

References  5  through  7  provide  the  following  information  and  data 
on  SEADYN- re la ted  topics:  Data  on  cables,  chain,  and  other  SEADYN  input 
parameters  (Ref  5);  a  summary  of  the  comparisons  between  SEADYN 
predictions  and  measured  data  (Ref  6);  and  a  bibliography  of  NCEL 
reports  published  on  cable  dynamics  (Ref  7). 

The  user  is  cautioned  that  only  the  cable  portion  of  the  program 
has  been  satisfactorily  validated.  The  mooring-related  options  are 
still  under  development  and  are  subject  to  change.  Therefore,  the  use 
of  SEADYN  for  problems  involving  vessel  dynamics  is_not_recoinnieji\ded  at 
this  time . 

THE  U.S.  GOVERNMENT  ASSUMES  NO  LIABILITY 
FOR  ANY  LOSS  OR  DAMAGE 
RESULTING  FROM  THE  USE  OF  THIS  PROGRAM. 
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1.0  INTRODUCTION 


The  SEADYN  computer  program  simulates  the  responses  of  cable  and 
truss  type  systems  in  an  ocean  environment.  Such  structures  often 
appear  deceptively  simple.  As  a  result  the  uninitiated  is  surprised 
that  careful  modeling  is  required  to  obtain  a  solution.  Cable  structures, 
even  a  simple  catenary  span,  are  highly  nonlinear.  Adding  the  offshore 
environment  greatly  compounds  the  problem.  This  computer  program  provides 
extensive  capability  for  dealing  with  these  problems  and  includes  con¬ 
siderable  flexibility  and  variety  in  the  solution  methods.  But  capability, 
flexibility,  and  variety  have  not  been  achieved  without  a  price  --  an 
increase  in  complexity. 

SEADYN  is  a  finite  element  program  that  employs  simple  line  elements. 
These  line  elements  can  be  used  repetitively  to  describe  long  mooring 
lines,  intricate  truss  structures,  or  complex  cable  systems.  SEADYN 
includes  lumped  body  models  to  simulate  buoys,  anchors,  weights,  etc.  in 
the  structure.  These  lumped  masses  are  treated  as  point  loads  that 
represent  weight/ buoyancy  and  fluid  drag  effects.  SEADYN  also  models 
dynamic  rigid  bodies  that  represent  ships,  barges,  platforms,  and  mooring 
buoys.  These  are  distinct  from  the  previously  mentioned  lumped  bodies 
since  they  provide  for  multiple  line  attachments  and  full  six  degree-of- 
freedora  body  response  calculations.  Two  additional  unique  features  of 
the  program  are  the  ability  to  represent  variable  length  lines  (payout/ 
reel-in)  and  the  component  adequacy  checks  (tests  buoys,  anchors,  and 
lines  for  potential  underdesign). 

The  complicated  nature  of  the  nonlinear  problems  dealt  with  by 
SEADYN  and  the  variety  of  numerical  methods  built  into  the  program  make 
it  tempting  to  expound  on  the  theories  involved.  Indeed,  it  is  felt 
that  a  program  such  as  SEADYN  can  be  used  intelligently  only  when  the 
user  has  some  fundamental  knowledge  of  nonlinear  mechanics  and  numerical 
solution  methods.  Attempts  to  treat  such  a  program  as  a  magic  black  box 
will  inevitably  lead  to  costly  piles  of  paper  full  of  useless  numbers. 
Unfortunately,  those  who  succumb  to  the  black  box  syndrome  are  often 
unaware  of  erroneous  or  meaningless  results  since  they  seldom  question 
what  the  computer  has  done.  To  such  a  user,  a  successful  execution  of 
the  program  with  no  diagnostic  messages  means  correct  answers.  A  success¬ 
ful  solution  of  nonlinear  equations  requires  the  selection  of  appropriate 
solution  procedures  and  the  specification  of  incrementing  and  or  iteration 
,  parameters.  Default  values  can  be  provided  when  the  user  knows  little 
else  to  do,  but  often  the  appropriate  method  and  parameters  are  problem 
dependent.  It  is  unrealistic  to  expect  a  single  manual  to  tell  all  that 
is  known  on  the  subject.  In  the  first  place  that  is  simply  impossible; 
even  if  the  body  of  knowledge  is  static  much  of  the  know  how  is  obtained 
through  experience  and  cannot  be  given  as  a  set  of  rules  or  instructions. 
Secondly,  the  bulk  and  detail  would  be  overwhelming  to  the  novice,  and 
much  that  is  important  would  be  diluted. 
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The  approach  taken  here  is  to  briefly  outline  the  nature  of  the 
nonlinear  problem,  describe  in  simple  terms  the  pertinent  features  of 
the  nonlinear  solution  methods,  present  the  input  form  for  SEADYN,  and 
illustrate  the  modeling  process  with  some  simple  problems.  Exposition 
of  the  governing  equations  and  theoretical  details  is  contained  in 
Reference  1. 
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2.0  CAPABILITY  OVKRVIEW 


The  SEADYN  program  is  dosigiied  primarily  to  analyze  the  static  and 
dyiKimic  responses  of  underwater  cable  and  truss  structures.  The  major 
fe>atures  of  the  program  are  summarized  in  Table  2-1.  The  structural 
load  paths  are  def.ined  by  straight  lines  between  node  points.  The 
elements  represented  by  these  lines  are  assumed  to  have  negligible 
bending  resistance,  and  torsional  effects  are  ignored.  Thus  the  only 
load  in  each  element  is  tiMision  or  compression.  An  option  is  provided 
to  have  compression  treated  as  negligible.  Should  any  element  in  the 
structure  reach  a  displaced  state  representing  a  shortening  of  that 
element  beyond  its  unstretchod  length,  the  stiffness  is  set  to  zero. 

The  stiffness  is  restored  when  a  stretched  state  is  obtained.  The 
physical  limitations  imposed  by  the  water  surface  or  the  seafloor  can  be 
approximated  by  setting  limits  on  specific  nodes  that  the  user  expects 
to  be  involved.  The  positions  of  the  nodes  are  continually  monitored  in 
nonlinear  solution  step.*'  to  activate  or  deactivate  these  tension-only 
element  and  limit  conditions. 

Two  basic  material  models  are  provided  to  represent  static  load 
behavior.  One  uses  a  load-strain  tabular  format  for  specifying  general 
nonlinear  behavior.  The  other  is  a  two-pa raine ter  model  t ha tj^p resumes 
curve  fitting  of  empirical  data.  The  assumed  form  is  T  =  ae  .  Material 
damping  effects  can  be  included  through  a  one  parameter  (Kelvin  damping) 
or  a  two-parameter  (Reiii-NOAA  damping)  viscoelastic  model.  A  proportional 
damping  ojition  is  also  provided. 

The  fluid  dissipative  and  loading  effects  on  the  structure  can  also 
be  modeled.  Reynolds  number-dependent  drag  and  constant  added  mass 
coefficients  are  provided  for  the  line  elements  and  lumped  bodies. 

These  program  coefficients  can  be  selected  or  the  user  can  input  his  own 
via  a  suhroutlno.  Dynamic  fluid-solid  interaction  on  rigid  bodies  is 
treated  using  drag  and  added  mass  coefficients,  buiit-in  potential 
theory  solutions  fer  mooring  buoys  in  frequency  domain  analyses,  and 
ej.ternally  supplied  data  files  for  ships  and  platforms;  only  drag  is 
considered  in  stat.ic  analyses.  An  approximate  static  loading  function 
is  also  provided  for  ships. 

The  solution  formats  in  SEADYN  can  be  assembled  in  varying  order. 
Using  the  bas.ic  assumption  that  large  deformations  occur,  the  behavior 
at  each  motion  state  is  dependent  on  its  current  configuration  and 
possibly  on  the  previous  history.  A  staged  or  seq 'cnced  solution  form 
can  be  used  to  apply  loads  and  pe.rform  calculations.  The  configuration 
at  the  beginning  of  one  stage  (called  a  subanalysis)  is  automatically 
defined  as  the  configuration  reached  at  the  end  of  the  previous  stage, 

A  typical  sequence  might  be  as  follows: 

St  ige  i:  Establish  a  st.'ucturally  stable  reference  state. 

(Usually  this  represents  the  static  state  subjected 
only  to  gravity  loads.) 

Stage  2:  Apply  steady-state  or  static  loadings  and  special 
boundary  conditions  (e.g.,  impose  wind  and  current 
loads) . 
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Stage  3:  Apply  dynamic  loads  and/or  impose  movement. 

Stage  4:  Kecovor  a  previously  calculated  state  (RESTART  option) 
and  impose  a  new  static  loading. 

Stage  5;  Apply  a  new  dynamic  load  sequence. 

Stage  6;  Calculate  mode  shapes  and  aatural  frequencies. 

Stage  7:  Continue  dynamic  analysis  of  Stage  S ,  etc. 

The  input  data  structure  of  SEADYN  is  segregated  into  a  problem 
description  set  and  a  sequence  of  subanalysis  data  .sets.  Each  subap.alysis 
(Stages  1  through  7)  is  identified  by  a  subanalysi.s  option  (SAO)  llag- 
■A  list  of  the  SAO  flags  and  a  brief  description  of  their  purpose  follow; 

DEAD  -  Nonlinear  static  analysis  to  apply  gravity/buoyancy  and 
specific  point  loadings. 

LiVE  -  Nonlinear  static  analysis  to  apply  f  hn  li  loads  and  various 
other  load  states, 

DYN  *  Nonlinear  dynamic  analysis  (time  domain).  Allows  time 
variations  in  loads,  currents,  motions  (towing,  etc:.), 
line  lengths  (payout/ reel- in) ,  and  lumped  body  impact. 

TSSS  -  Time  Sequenced  Static  Solutions.  Approximates  nonlinear 
dynamics  by  neglecting  accelerations  and  solves  for 
configuration  clianges  due  to  payout/reel-in  and/or  imposed 
boundary  motion  as  a  sequence  of  static  nonlinear  solutions. 
Generates  a  LIVE  solution  for  each  time  step  requested. 

MODE  -  Calculates  vibration  mode  shapes  and  natural  fi;?quencies 

for  a  linearized  representation  of  the  current  state.  This 
provides  information  only.  The  mode  data  are  not  used  by  any 
other  subanalysis  now  defined  in  SEADYN. 

FREQ  -  Frequency  domain  dynamic  response  using  linearized  frequency 
dependent  solutions.  Calculates  regular  wave  response  or 
response  spectra  at  any  node  for  specified  wave  spectra 
imposed  on  ships  and  mooring  buoys. 

CHEK  -  Evaluates  the  adequacy  of  various  mooring  components  in  the 
currently  defined  state. 

The  LIVF.,  TSSS,  and  DYN  SAO  sets  allow  for  adjustment  of  drag  loading 
to  approximate  the  effects  of  line  strununing,  Tlus  is  done  in  llie  form 
of  response-dependent  drag  amplification  factors  that  are  periodic.jl ly 
calculated  for  user-identified  strings  of  elements. 


The  FRKQ  SAO  set  can  use  an  implied  LIVE  SAO  to  approximate  the 
change  in  static  position  due  to  wave-induced  steady  drift  effects. 

A  restart  capability  is  provided  to  allow  the  saving  and  subsequent 
use  of  of  the  calculated  states.  A  restart  can  be  in  the  current 
nni  (as  in  tlie  earlier  example)  or  in  some  subsequent  computer  run.  The 
restart  files  are  also  used  as  a  data  base  for  a  separate  post-processor 
graphics  program  SKAPLT',  Reference  4. 

SEADYN  provides  a  minimal  set  of  drag  coefficients,  time  variation 
functions,  current  profiles,  and  ship-load  functions.  When  these  prove 
inadequate,  alternate  definitions  can  be  provided  through  user  subroutine 
and  data  files.  Frequency  domain  dynamic  equations  are  provided  for 
mooring  buoys  and  other  bodies  that  can  be  approximated  by  a  simple 
sphere  immersed  to  the  equator.  More  complex  body  dynamics  can  be 
described  through  an  externally  generated  ship  motion  data  file. 


Table  2-1.  SEADYN  Capability  Summary 


General 

•  3D  large  displacement  response  of  cable  and  truss-type  structures 
using  the  finite  element  method 

c  Can  include  lumped  bodies,  six  degree-of-freedom  rigid  bodies 
(ships,  platforms,  etc.)  and  fluid-solid  interactions 

•  Staged  format,  sequential  analysis  for  statics  and  dynamics 

•  Treats  nonlinear  materials  with  internal  damping,  nonconservative 
loads,  and  nonlinear  constraints 

Special  Features 

•  Variable  length  lines  to  simulate  payout/reel-in  dynamics 

•  Automatic  estimation  of  drag  coefficient  amplification  to 
approximate  strumming 

•  Restart  options 

•  Wave  spectrum  analysis  with  approximation  for  steady  drift  forces 
«  Component  adequacy  checks  using  design  rules 

•  Plotting  interface 

•  Free-field  input  format 

•  Catenary  element  for  treating  bottom  interaction 
Load/Boundary  Conditions 

•  Gravity/buoyancy  loads  in  water  and  air 

•  Arbitrary  point  loads  and  flow  fields 

•  Wind  and  surface  current  loads  on  rigid  bodies 

•  Built-in  or  user-supplied  drag  functions  with  flow/response- 
dependent  amplification 

•  Arbitrary  time  variations  (built-in  or  user-supplied) 

•  Moved  boundaries 

•  Conditional  constraints  for  surface  and  bottom  limits 
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Static  Solution  Methods 

•  Sequence  of  linear  increments  (Euler's  method) 

•  Residual  feedback  method  (incremental  self-correcting) 

•  Modified  Newton-Raphson  (various  forms) 

•  Viscous  relaxation  method 
Dynamic  Solution  Methods 

•  Nonlinear  Transient  -  Sequence  of  Linear  Increments 

-  Newmark’s  p  (residual  feedback  form) 

-  Direct  Integration  Method  (a  multi-para- 

meter  predictor/corrector) 

-  Time  sequenced  static  solutions  (quasi¬ 

static) 

•  Frequency  Domain  -  (linearized  with  respect  to  the  static  state) 

-  Mode  shapes  and  frequencies 

-  Response  to  wave  spectra  -  Superposition  of 

frequency-dependent  steady-state 
solutions,  fully  coupled  ship,  buoy 
and  line  responses 
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3.0  THEORY  OVERVIEW 


The  SEADYN  program  seeks  to  represent  a  general  spatial  arrangement 
of  cable  and  truss  components  as  a  collection  of  simple  elements.  The 
present  version  of  SEADYN  allows  two  element  forms:  a  straight  line 
between  two  points  (nodes),  and  a  bottom- limited  catenary.  The  straight 
line  element  is  straight  before  and  after  deformation  of  the  structure, 
it  assumes  negligible  bending  resistance  and  uses  the  instantaneous 
distance  between  the  defining  nodes  along  with  the  unloaded  (unstretched) 
length  of  the  element  to  estimate  the  strain.  A  single  value  for  the 
element  tension  is  then  obtained  from  the  Le/ision/strain  data  given  for 
that  element.  Since  straight  lines  and  constant  tensions  are  low  order 
approximations  to  the  actual  behavior  of  flexible  lines,  it  is  important 
that  the  analyst  carefully  consider  how  the  structure  is  modeled. 

Region.s  where  large  curvature  and/or  highly  variable  tensions  are  expected 
should  be  approximated  with  more  elements  than  in  those  regions  where 
tensions  vary  slightly  and  lines  are  more  nearly  straight.  In  general, 
more  elements  mean  more  accuracy  and  more  cost.  Good  modeling  practice 
leads  to  a  rational  compromise  between  accuracy  and  cost.  Achieving 
this  compromise  requires  engineering  judgment  that  comes  from  experience. 

The  bottom-limited  catenary  element  uses  classical  catenary  equations 
in  an  iterative  procedure  to  develop  stiffness  and  force  relationships 
for  a  line  which  must  interact  with  the  bottom.  Two  nodes  are  used  in 
the  element  definition,  and  adjustments  are  made  in  the  equations  to 
deal  with  line  stretch  and  the  amount  of  line  lying  along  the  bottom. 

Only  gravity  loads  are  applied  to  this  element  (current  loads  are  neglected). 
Only  approximate  mass  relations  are  provided  for  treating  the  dynamics 
of  line  pick-up  and  lay-down,  so  this  element  should  be  used  with  caution 
in  transient  dynamic  solutions. 

Since  SEADYN  deals  with  large  deflection  effects,  the  position  and 
velocity  of  all  nodes  and  the  unstretched  lengths  of  all  elements  roust 
be  considered  in  each  step  of  the  solution.  This  is  referred  to  as  a 
geometric  nonlinearity,  and  it  poses  some  problems  not  encountered  in 
small  displacement  analyses.  In  static  analyses,  where  inertia  effects 
are  ignored,  it  is  possible  to  have  a  set  of  unstretched  lengths  and 
nodal  positions  that  represent  an  unstable  structural  configuration. 

Unless  the  specified  loads  and  solution  procedurc.s  produce  appropriate 
movements  to  modify  the  position  of  the  nodes  and  reorient  the  elements, 
the  structure  represents  a  mt'chanisra.  This  means  that  it  is  not  capable 
of  providing  a  .static  load  path  between  the  points  where  loads  are 
applied  and  where  the/  are  supported.  Another  form  of  instability  can 
also  occur.  Tliis  is  the  more  familiar  buckling  instability  in  which  an 
apparently  stable  structural  configuration  will  suddenly  deform  to  a 
radically  different  shape  with  only  minor  changes  in  loading. 

The  majority  of  numerical  solution  techniques  used  in  SEADYN  can  be 
cla.ssified  as  initial  value  methods.  This  means  that  a  solution  step 
proceeds  from  a  state  where  all  pertinent  data  are  presumed  to  be  known 
to  a  state  where  estimates  are  made  of  the  effects  of  loading  changes 
using  some  sort  of  a  predictor.  Wien  this  estimate  is  within  certain 
reasonable  bounds  of  accuracy  it  can  be  improved  by  iterative  corrections. 

In  those  situations  where  the  initial  state  is  not  accurately  described 
or  where  it  represents  an  unstable  state,  the  predictor  .is  usually  very 
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inaccurati? ,  if  not  undofined.  SKADYN  provid«!;  various  moans  of  dealing 
with  this  problem,  but  uni ortimately  there  is  no  all  pow<^rful  method 
that  works  every  time  with  little  or  no  special  input.  IL  is  in  dealing 
with  this  proble'v  that  the  new  user  (and  often  the  experienced 
will  encounter  the  most  frustration. 

The  various  solution  methods  are  discussed  in  the  following  para¬ 
graphs.  The  user  will  readily  find  the  methods  that  work  effectively  on 
well-posed  problems .  In  these  cases  more  than  one  approacii  may  be 
available,  and  the  choice  of  which  one  to  use  can  he  based  on  economics, 
accuracy,  or  user  experience.  The  discussion  focuses  on  difficulties 
that  can  be  encountered  in  highly  nonlinear  or  poorly  posed  problems  in 
the  hope  that  the  user  can  identify  the  difficulty  and  see  how  to  adjust. 

The  first  step  in  any  analysis  is  to  describe  the  isiitial  geometric 
state.  As  noted  above,  thi.s  requires  the  specification  of  the  length 
(stretched  or  unstretclied)  of  every  element  and  the  spatial  position  of 
each  node.  The  ideal  is  to  have  a  consistent  set  of  nodal  positions  and 
element  lengths  that  describe  an  equilibrium  state  under  a  known  loading. 
The  natural  way  of  describing  such  a  state  is  to  give  the  spatial  coordi¬ 
nates  of  each  node  and  the  element  tensions.  Unstretched  element  lengths 
are  then  easily  calculated  from  the  strains  (obtained  from  the  tensions 
and  material  properties)  and  the  stretched  lengths  (obtained  from  the 
nodal  positions).  When  such  a  configuration  description  is  available, 
the  solution  process  can  start  simply  by  specifying  a  loading  sequence 
that  begins  with  the  loads  represented  by  the  initial  equilibrium  state. 

Most  often,  the  analyst  does  not  have  an  exact  equilibrium  state 
from  which  to  start.  Initial  unstretched  element  lengths  are  known  and 
the  general  form  of  the  structural  layout  can  be  defined,  but  the  correct 
nodal  positions  and  tensions  for  a  particular  loading  are  not  known. 

For  simple  structural  forms  it  may  be  possible  to  guess  a  set  of  nodal 
positions  that  approximates  the  equilibrium  state  with  the  distance 
between  the  nodes  representing  the  unstretched  rather  than  the  stretched 
lengths  (catenary  response,  for  example).  Generally,  oven  this  approximate 
shape  will  be  awkward  to  specify.  It  may  be  more  convenient  to  list  the 
unstretched  element  lengths  independent  of  the  positions  of  the  nodes  in 
the  guessed  configuration.  In  either  of  these  situations  the  element 
tension  data  are  extraneous  because  the  correct  tensions  can  be  obtained 
only  after  the  correct  nodal  positions  are  obtained.  In  the  first  case, 
the  nodal  positions  imply  unstretched  lengths,  which  in  turn  means  all 
the  element  tensions  are  zero.  In  the  second  case,  the  nodal  positions 
imply  guesses  of  the  stretched  lengths  which,  when  taken  with  the  speci¬ 
fied  unstretched  lengths,  imply  a  uonequi 1 ibrium  set  of  tensions. 

Neither  of  these  approaches  will  give  a  stable  state  from  which  to 
proceed  directly  through  an  analysis  sequence.  An  equilibrium  state  for 
the  structural  configuration  must  be  found  for  some  initial  loading 
state.  When  the  starting  point  is  not  an  equilibrium  state,  the  initial 
value  solution  methods  fail;  the  predictor  becomes  inaccurate  or  undefined. 
This  is  true  in  static  analyses  because  the  predictor  is  expressed  in 
the  form  of  a  tangent  stiffness  matrix  that  is  singular  unless  the 
structure  is  stable.  However,  it  is  possible  to  construct  a  very  approxi¬ 
mate  predictor  for  unstable  systems  using  a  set  of  artifical  teiisions. 
SKADYN  provides  for  input  of  artificial  tensions  and/or  a  psiMido-tensioii 
generator  called  numerirnl  damping.  The  predictors  constructed  in  this 
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manner  can  be  very  inaccurate,  and  the  appropriate  choices  for  the  set 
of  tensions  and  numerical  damping  parameters  are  highly  problem  dependent. 
Some  trial  solutions  may  be  required  to  find  the  appropriate  choices  for 
a  specific  structure. 

Various  numerical  solution  schemes  are  provided  that  attempt  to 
converge  on  the  correct  equilibrium  state  when  the  structure  is  initially 
unstable  or  nearly  so.  All  of  them  rely  on  some  combination  of  iteration 
and  incrementing  (applying  the  load  in  steps).  The  multiplicity  of 
methods  is  maintained  in  SEADYN  because  no  universally  effective  method 
has  been  found.  Most  problems  can  be  handled  quite  well  once  the  analyst 
has  experimented  with  the  solution  parameters  and  the  specific  problem 
to  identify  and  adjust  to  any  sensitivities.  One  should  plan  to  spend  a 
little  time  becoming  familiar  with  the  initial  responses  in  highly  non¬ 
linear  problems.  The  discussions  in  Sections  8  and  9  will  be  helpful. 

Static  Solution  Methods 

A  brief  description  of  the  static  solution  methods  follows.  More 
detailed  descriptions  are  given  in  References  8  and  9.* 

Residual  Feedback  Method  (RFB).  This  is  an  incremental,  self-correct¬ 
ing  procedure  that  presumes  loads  or  imposed  displacements  are  applied 
in  a  sequence  of  steps.  The  first  step  is  a  simple  linear  one.  Each 
step  after  the  first  adjusts  the  load  increment  by  an  estimate  of  the 
equilibrium  error  from  the  previous  step.  The  incremental  stiffness 
matrix  is  recalculated  at  each  step  and  reflects  the  nonlinearities 
apparent  at  the  last  completed  step.  This  method  is  the  least  expensive 
to  apply,  the  most  numerically  stable,  and  the  least  accurate.  It 
requires  explicit  listing  of  the  steps  and  performs  only  that  number  of 
steps  with  no  iteration.  It  requires  a  stable  state  to  start  from, 
which  can  be  artifically  produced  by  fictitious  tensions  or  initial 
numerical  damping.  The  RFB  method  can  fail  if  the  incremental  stiffness 
matrix  becomes  singular  or  very  ill-conditioned  after  the  first  step. 

The  final  state  can  have  large  equilibrium  errors  since  a  full  iteration 
to  convergence  is  never  done.  More  accurate  solutions  can  be  obtained 
by  increasing  the  number  of  steps. 

Modified  Newton-Raphson  Method  (MNR) .  This  method  can  be  used  in  a 
fully  iterative  (single-step)  form  or  an  incremental-iterative  form. 

The  method  iteratively  evaluates  the  difference  between  the  external 
nodal  forces  and  the  internal  reactions  to  search  for  the  displaced 
state  that  saLisifies  equilibrium.  The  tangential  stiffness  matrix  is 
used  to  estimate  the  displacement  change.*:  for  each  iteration.  This 
matrix  can  be  recalculated  at  each  step  or  at  user-specified  intervals. 

When  the  incremental  form  is  used,  the  user  must  estimate  the  number  of 
steps  to  be  used.  Depending  on  the  rate  of  convergence  at  each  step, 
the  step  size  can  be  increased  or  decreased  automatically  by  the  program. 
Various  schemes  are  used  on  the  iterations  to  accelerate  convergence. 


*Another  solution  form  is  also  in  SEADYN.  It  is  called  a  sequence  of 
linear  increments  (SLI),  This  is  essentially  Euler's  method  and  it 
requires  very  small  steps  to  get  accuracy.  It  is  of  no  value  in 
getting  initial  configurations,  and  its  use  is  not  recommended. 
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The  first  iteration  at  each  load  step  can  extrapolate  from  the  converged 
state  in  the  previous  step.  Fictitious  tensions  can  be  used  in  the 
starting  state,  and  numerical  damping  can  be  employed  on  each  step.  The 
MNR  method  is  conditionally  stable.  It  can  diverge  if  the  predictor  and 
convegence  accelerators  are  inefficient  or  inappropriate.  The  convergent 
behavior  can  be  obtained  very  slowly  with  excessive  numerical  damping. 
Convergence  of  the  iterations  is  determined  from  simultaneous  occurrences 
of  small  displacement  changes  and  small  values  for  the  force  residual 
(difference  between  external  and  internal  forces).  The  user  can  specify 
the  convergence  tolerance,  the  frequency  of  tangent  stiffness  matrix 
calculations,  the  limits  on  iterations,  and  the  characteristics  of  the 
extrapolations  and  convergence  accelerators.  Default  values  are  taken 
when  no  selection  is  made. 

Viscous  Relaxation  Method  (VRR) .  This  is  a  generalized  form  of  the 
Newton-Raphson  method,  which  can  automatically  adapt  the  characteristics 
of  the  solution  steps  to  the  behavior  being  sensed.  An  artificial  time 
parameter  is  used  to  produce  load  steps  or  iterations.  When  applying 
the  load  in  steps,  a  series  of  numerically  damped  load  steps  are  produced 
that  have  a  form  similar  to  the  RFB  method.*  When  the  full-load  level 
is  reached,  the  VRR  method  iterates  while  adjusting  the  damping  level 
and  pseudo-time  step  to  move  to  the  equilibrium  state.  As  convergence 
is  achieved,  the  method  degenerates  to  the  Newton-Raphson  method  with  a 
stiffness  matrix  evaluation  at  each  iteration.  This  is  the  most  robust 
and  the  most  expensive  method.  It  is  often  capable  of  getting  convergent 
solutions  when  the  others  fail.  The  VRR  method  can  fail  if  there  are 
excessive  numbers  of  iterations  or  an  inappropriate  selection  of  control 
parameters.  The  methods  used  to  adjust  the  solution  characteristics  are 
heuristically  defined  and  can  work  well  on  some  problems  and  be  ineffective 
on  others.  The  user  can  select  the  initial  damping  level,  the  initial 
step  size,  the  number  of  starting  increments,  and  the  convergence  tolerances. 
Default  values  are  provided,  but  the  appropriate  choice  of  parameters  is 
problem  dependent.  Inappropriate  selections  of  damping  can  lead  to  very 
slow  convergence  or  wildly  oscillating  behavior  (cyclic  calculations 
without  convergence).  Convergence  is  signaled  by  a  very  low  value  of 
the  force  residual  or  by  the  simultaneous  occurrence  of  low  force  residual 
and  small  pseudo-velocities  at  all  of  the  nodes. 

Even  when  a  stable  structural  state  has  been  obtained,  it  is  possible 
to  develop  solution  instabilities  while  subsequently  applying  additional 
loads  and/or  movement.  These  problems  result  from  physical  instabilities 
in  the  structure  (e.g.,  buckling)  or  from  numerical  instabilities  in  the 
solution.  Often  both  situations  occur  together.  In  general  one  does 
not  encounter  the  classical  bifurcation  buckling  behavior  exemplified  by 
an  Euler  column.  An  instantaneous  shift  in  deformation  behavior  is  a 
small  difiplacement  ab.straction  of  a  large  displacement  nonlinear  phenomenon. 
Incicmental/iterative  solutions  would  show  a  buckling  instability  as  a 
large  change  in  deflection  for  a  small  change  in  load.  Graphically  this 
would  appear  as  a  strongly  curved  load/deflection  plot  with  a  slope 
tending  to  zero.  Analytically  this  means  a  portion  of  the  stiffness 
matrix  contains  small  coefficients  relative  to  the  rest  of  the  matrix. 
Numerically  this  means  the  stiffness  matrix  is  ill-conditioned  (nearly 
singular) . 


*A  parallel  form  which  employs  the  SLI  format  is  referred  to  as  the  VRS  method. 
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In  modeling  cable  structures  (moorings  in  particular),  one  often 
encounters  portions  of  the  structure  with  stiffnesses  much  lower  than 
the  rest  of  the  system.  Lines  carrying  very  low  tension  have  very 
little  stiffness  in  lateral  motion  even  though  their  stiffness  in  stretching 
can  be  high.  The  angular  motion  components  on  moored  vessels  and  buoys 
can  have  much  lower  stiffnesses  than  the  linear  motion  components.  This 
disparity  in  stiffness  may  not  always  signal  a  buckling-type  problem, 
blit  it  does  produce  numerical  problems.  The  low  initial  tension  problem 
is  not  a  buckling  problem  unless  compressive  loads  are  applied.  Loads 
that  produce  tension  will  stabilize  the  structure.  (This  is  the  initial 
state  problem  discussed  earlier.)  On  the  other  hand,  the  buoy  angular 
response  problem  can  well  be  a  buckling  problem.  There  can  be  loading 
conditions  that  drive  the  buoy  angular  response  to  give  large  angle 
changes  for  small  changes  in  load  magnitude  and/or  direction.  Or  there 
can  be  limit  states  (which  represent  a  sort  of  buckling)  where  the  buoy 
position  is  unstable  and  tends  to  move  into  a  state  representing  a 
twisting  of  the  mooring  lines.  SEADYN  makes  no  attempt  to  detect  line 
crossings  or  related  phenomena. 

Attempts  to  converge  on  an  equilibrium  position  when  very  low 
stiffness  components  are  in  the  system  can  be  costly  (requiring  many 
iterations  and  small  steps).  The  buoy  angular  response  and  low  tension 
systems  are  particularly  troublesome.  In  both  cases,  the  iterations  can 
produce  large  changes  in  position,  and  convergence  will  be  hard  to 
obtain.  The  problem  is  compounded  by  numerical  errors  associated  with 
solving  the  stiffness  equations.  A  mixture  of  stiff  and  soft  components 
leads  to  numerical  ill-conditioning  errors,  which  can  be  compounded  by 
the  sequence  of  the  equation  processing.  This  sequence  is  determined  by 
the  node  numbering  scheme  in  SEADYN.  The  ideal  sequencing  causes  the 
soft  components  to  be  processed  before  the  stiff  ones.  Unfortunately, 
one  cannot  afford  the  luxury  of  optimum  ordering  in  this  sense  because 
it  can  greatly  increase  the  equation  bandwidth,  thereby  increasing  the 
solution  costs.  Consult  Reference  10  for  more  on  equation  ordering  and 
solution  errors. 

Rigid  bodies  (mooring  buoys,  ships,  platforms,  etc.)  are  interfaced 
to  cable  elements  using  multiple  attachment  points.  Nodes,  which  are 
used  on  the  cable  elements,  have  three  degrees-of-freedom  (one  for  each 
displacement  component).  The  rigid  bodies  require  six  degrees-of-freedom 
(three  displacements  and  three  angular  components).  Slave/master  constraints 
are  imposed  that  imply  the  rigid  body  kinematics.  A  pair  of  nodes  is 
used  to  define  the  six  degrees-of-freedom  for  a  rigid  body  and  all  other 
nodes  on  the  body  are  required  to  move  as  though  rigid  links  exist 
between  them  and  the  body  node  pair.  The  body  node  pair  is  the  master, 
and  all  other  nodes  on  the  body  are  rigid  link  slaves. 

Lumped  bodies  are  assumed  to  have  insignificant  spatial  dimensions 
so  that  their  effect  on  the  kinematics  of  the  line  is  negligible.  They 
are  simply  lumped  at  nodes  where  they  produce  mass  anl  drag  load  effects. 

Dynamic  Solution  Methods 

The  dynamic  solution  procedures  in  SEADYN  fall  into  two  categories: 
time  domain  and  frequency  domain.  Time  domain  analyses  are  fully  nonlinear. 
Two  basic  solution  methods  are  available  for  numerical  integration  of 
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the  nonlinear  time  domain  equations.  Both  solutions  qre  based  on  a 
generalized  form  of  the  Newmark  difference  equations  (Ref  8).  At  present 
there  are  no  time  domain  rigid  body  equations  available  in  SEADYN  so 
that  nonlinear  Lime  domain  solutions  for  ships,  platforms,  and/or  buoys 
are  not  possible.  Master/slave  relations  are  not  functional  in  the  time 
domain.  A  brief  description  of  time  domain  solution  mettiods  follows. 

Implicit  RFB  Method.  This  is  an  implicit  integration  scheme  that 
follows  the  more  traditional  Newmark  format  of  solving  a  set  of  simul¬ 
taneous  algebraic  equations  at  each  time  step.  Payout/reel-in  and 
moving  boundary  options  have  not  been  implemented  in  this  method. 
Specification  of  three  integration  parameters  and  the  time-  step  size  is 
required.  The  method  is  strongly  stable  but  can  be  inaccurate  for  large 
time  steps. 

Direct  Integration  Method  (DIM).  This  is  a  predictor-corrector 
technique  that  does  not  require  the  formation  of  a  stiffness  matrix. 
Specification  of  three  integration  parameters  and  iteration  controls  is 
required.  Time  step  size  can  be  specified  or  calculated  by  the  solution 
routine.  The  iterative  corrector  is  conditionally  convergent  and  requires 
strict  upper  bounds  on  the  time  step. 

There  are  two  types  of  frequency  domain  analyses  in  SEADYN.  One  of 
the  frequency  domain  options  is  simply  the  estimation  of  vibration  mode 
shapes  and  natural  frequencies.  It  uses  the  well  known  Jacobi  method. 

This  option  gives  information  about  the  dynamic  character  of  the  equations. 
The  information  is  made  available  on  a  data  file  but  is  not  used  elsewhere 
by  SEADYN.  In  the  present  version  of  the  program,  only  a  diagonal  mass 
matrix  is  used,  and  no  correction  is  made  for  the  lack  of  a  tangential 
added  mass  on  the  cable  elements.  All  other  dynamic  options  in  SEADYN 
make  this  adjustment. 

The  other  frequency  domain  option  is  a  quasi-linear  procedure  that 
solves  a  set  of  frequency-dependent  linearized  equations  for  steady-state 
harmonic  responses.  These  responses  represent  the  fully  coupled  behav'  r 
of  mooring  components  and  a  single-moored  vessel  subjected  to  wave 
forces  defined  through  a  wave-height  spectrum.  Spectral  superposition 
techniques  are  used  to  estimate  response  spectrum  data.  The  methods 
•ised  assume  that  the  motions  involved  are  small  amplitude  perturbations 
about  a  nonlinear  static  reference  state.  The  steady  effects  of  wave- 
induced  drift  forces  can  be  approximated  by  an  iterative  procedure  that 
automatically  generates  a  LIVE  subanalysis.  The  LIVE  subanalysis  adds 
the  drift  forces  onto  the  static  reference  state  after  a  frequency 
spectrum  pass  has  been  made.  The  user  has  the  option  of  limiting  this 
procedure  to  a  single  pass  or  repeating  it  until  no  significant  changes 
occur  in  drift  forces  and  wave  responses.  See  References  11  and  12  for 
theoretical  details.  SEADYN  has  frequency  domain  equations  for  the 
response  of  spherical  mooring  buoys.  These  equations  assume  that  the 
w.-itcr  line  is  at  the  equator.  It  is  assumed  that  the  frequency  domain 
equations  for  the  ship  or  platform  are  available  in  a  user-provided  file 
that  gives  motion  coefficients  and  loading  functions  versus  wave  frequency 
and  heading.  The  format  of  this  file  is  described  in  Appendix  A. 
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4.0  INPUT  DATA  STRUCTURE 


The  SEADYN  program  uses  a  free-field  input  format.  Information 
describing  the  problem  and  subanalyses  is  organized  into  blocks  of  data 
called  records ,  terminated  by  column  80  or  a  terminating  character  (see 
Section  5).  Examples  of  records  are  punched  cards  or  lines  separated  by 
a  return.  General  rules  for  generating  SEADYN  data  records  are  given  in 
Section  5,  followed  by  specific  descriptions.  The  data  deck  always 
begins  with  at  least  one  title  record.  An  unlimited  number  of  title 
records  can  be  used  at  the  beginning  of  the  deck.  Later  in  the  deck 
structure  when  the  logic  calls  for  a  new  title  record,  only  one  is  to  be 
given.  Immediately  following  the  title  record  set,  SEADYN  expects  to 
find  either  a  RESTART  or  PROBLEM  data  set.  These  provide  a  basic 
description  of  the  problem  either  from  a  saved  file  or  new  data  input. 

The  PROBLEM  data  set  consists  of  a  variable  set  of  records  defining  the 
geometry,  materials,  and  various  tables  to  be  used  in  the  problem  solution 
stages.  The  data  sets  are  identified  by  flag  records  or  keywords. 

These  keywords  are  limited  to  10  characters  in  length,  but  only  the 
first  four  are  required. 

The  various  solution  stages  are  described  in  SAO  data  sets  under 
the  headings  described  previously.  A  termination  of  a  string  of  solution 
stages  is  indicated  by  an  END  or  NEW  record.  END  terminates  the  run. 

NEW  returns  to  the  title  record  set  (only  one  record  this  time)  and 
expects  to  find  RESTART  or  new  PROBLEM  directives.  The  basic  program 
flow  is  represented  in  Figure  4-1.  Representative  input  decks  are 
included  later  in  the  examples.  These  decks  provide  useful  references 
for  the  input  instructions  section. 

The  free-form  input  allows  the  data  records  for  a  ship  loading  file 
to  be  read  in  the  fixed-field  format.  When  required,  this  data  can  be 
included  anywhere  in  the  data  deck  after  the  first  title  record  set. 

The  input  format  is  described  in  Appendix  B. 
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Figure  4-1.  Macro  flow  chart  of  SEADYN, 


5.0  FREE-FIELD  INPUT  RULES 


The  following  special  characters  are  recognized  by  the  subroutine 

FREINP: 

$  Record  Terminator  Flag 

Signals  no  more  data  to  be  read  for  the  record  image  being 
processed.  Multiple  record  images  (records)  can  appear  on 
a  single  record  separated  by  record  terminators.  Double¬ 
record  terminators  signal  the  end  of  a  record,  and  any  data 
following  this  are  treated  as  a  comment.  Comments  will  be 
listed  as  part  of  the  record  but  will  not  be  transmitted  to 
the  data  file. 

;  Alternate  Record  Terminator  Flag 

Performs  same  function  as  $ 

C0LUMN80  Default  Record  Terminator 

Unless  a  prior  termination  or  continuation  is  signaled, 
the  end  of  record  (C0LUMN80)  is  taken  as  a  record 
termination. 

,  Word  Delimitor  (Separator) 

Separates  sequences  of  data  entries  in  a  record  image. 
Repeated  delimiters  produce  zeros  in  the  words.  An  initial 
comma  produces  a  zero  in  the  first  word  of  the  record.  All 
words  not  explicitly  defined  are  assumed  to  be  0.0.  A 
comma  can  be  used  to  signal  multiple  records  (continuation) 
when  only  blanks  occur  between  the  last  comma  and  the  end  of 
record. 

BLANK  Separator/Delimitor 

Leading  blanks  are  ignored.  Once  the  beginning  of  a  word 
is  detected,  a  blank  will  terminate  the  word.  Any  blanks 
following  a  delimitor  are  treated  as  leading  blanks  for 
the  next  word.  The  following  are  equivalent: 


XX  yy 
XX  ,  yy 
XX, yy 

Continuation  Flag 

Signals  a  word  termination  with  the  next  word  to  be  read 
from  the  next  record.  See  for  alternate  continuation. 

W  Word  Position  Flag 

Used  to  override  the  word  sequencing  and  shift  to  a  new 
word  in  the  record.  Input  then  follows  in  sequence  from 
the  new  word  location.  The  new  word  number  is  given 
immediately  following  the  W  and  before  the  next  or  blank. 
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The  W  can  be  used  as  a  delimitor  of  the  previous;  word  on 
alJ  but  the  first  word  in  a  record.  The  combination  ",W" 
is  the  .same  as  W  alone.  The  first  word  of  a  record  is  not 
che^lccd  for  the  W  flag  so  ",W"  mu.st  be  used  to  skip  to  a 
ni  setpience  from  the  first  word  po.sition.  Any  W  after  the 
first  word  and  before  the  record  terminator  will  be 
interpreted  as  a  position  flag. 

Comment  Record  Flag 

This  ch.aracter  anywhere  on  a  record  will  terminate  it 
and  the  remainder  of  the  record  is  treated  as  a  comment. 

*  in  column  1  produces  no  record. 

'  Hollerith  (Alpha-Numeric)  Field  Delimitor 

This  character  is  the  apostrophe  (11-8-5  on  the  026 
keypunch  and  8-5  on  the  029  keypunch).  It  signals  the 
st.irt  and  end  of  a  character  string.  The  string  can  be 
any  length  up  to  a  maximum  of  ten-character  "words" 
allowed  by  the  application.  This  is  usually  8,  but  it 
can  be  more  or  less  as  the  use  dictates.  Any  legal 
character  can  be  used  in  the  string  except  the  apostrophe. 
Character  strings  can  continue  past  column  80  to  the  next 
record  since  the  string  automatically  signals  continuation 
until  the  terminal  apostrophe  is  found.  All  blanks  in  the 
string  are  counted  as  characters. 

(  Fired  Format  Initiator  Flag 

This  in  column  1  of  any  record  after  the  title  record  signals 
that  the  records  up  to  the  next  ")"  record  are  in  fixed  form 
These  records  are  written  on  a  special  data  file  in  BCD 
format . 

)  Fixed  Forma t  Te iminator  Flag 

This  in  column  1  signals  the  end  of  a  sequence  of  rigid 
format  records. 

Any  record  with  a  or  in  column  1  will  be  treated  as  a  comment 

record.  It  will  be  listed  but  will  not  produce  a  data  record. 

Each  free-field  input  deck  is  presumed  to  begin  with  one  or  more 
title  records.  Title  records  are  read  and  listed  until  a  specific  record 
terminator  is  diJtected  ($  or  ;  but  no_t  COLUMN80) .  The  record  on  which  the 
terminator  is  detected  will  be  used  as  a  page  heading  for  the  run. 
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The  KPEINP  subroutine  processes  the  entire  input  deck  and  translates 
it  into  a  series  of  data  records.  As  noted  above,  a  data  record  can 
span  more  than  one  line  or  there  can  be  one  or  more  records  on  a  line. 

After  the  title  line,  the  data  are  assumed  to  be  arranged  in  blocks 
headed  by  a  flag  record.  Each  flag  record  is  limited  to  ten  alpha-numeric 
characters.  Only  the  first  four  characters  of  the  flag  record  have 
meaning.  For  example,  the  flag  record  ELEMENTS  could  be  shortened  to 
ELEM  to  produce  the  same  result.  In  those  instances  where  the  flag  has 
only  three  characters,  there  must  be  a  blank  or  word  terminator  after 
the  third  character.  The  specific  data  order  applicable  to  the  flag 
record  is  assumed  until  the  next  fl.ag  record  is  detected.  Flag  records 
must  have  the  flag  word  in  the  first  word  position. 

Data  records  are  assumed  to  be  in  floating  point  form  unless  a 
character  is  detected  that  is  inconsistent  for  a  floating  point  number. 

In  this  case,  the  word  will  be  treated  as  a  Hollerith  word.  All  floating 
point  words  assume  a  decimal  at  the  end  of  the  word  if  none  is  given. 

Words  actually  intended  to  be  integers  are  converted  to  a  fixed-point  form 
at  the  time  they  are  used  by  the  program.  The  maximum  length  of  a  data 
record  is  determined  by  the  program  using  the  free-field  subroutine. 

The  beginning  and  end  delimiters  for  fixed  format  records  must 
appear  in  column  1  of  a  distinct  data  record  (individual  record).  This 
specifically  requires  that  the  previous  record  must  have  been  appropriately 
terminated  (no  continuation).  The  (  )  delimitors  are  the  only  things 
read  on  that  record,  and  the  next  data  processed  are  assumed  to  be  on  the 
next  record.  Only  one  rigid  format  data  set  is  allowed  in  any  run.  A 
rigid  format  data  set  cannot  be  input  before  the  initial  title  record 
set  is  completed. 
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6.0  PROBLEM  DEFINITION  DATA 


Thv.L^  are  two  w.ays  of  providing  problem  definition  data:  RESTART 
or  a  new  PROBLEM  data  set. 

RESTART  presumes  a  previously  generated  restart  file  is  available 
that  contains  all  of  the  problem  description  data  and  the  results  of 
the  SAO  calculations.  RESTART  can  be  use^’  to  continue  the  previous  SAO 
with  some  of  the  options  and  parameters  t‘  nged,  if  desired.  At  the 
completion  of  this  SAO,  any  other  appropriate  SAO  can  be  requested. 
Alternatively,  RESTART  can  be  used  to  establish  a  starting  state  for  a 
new  SAO.  The  appropriate  form  for  RESTART  input  is: 


Title  Record  Set 
RESTART 
(Restart  data) 

SAO  flag 

(SAO  data  as  needed) 
END  or  NEW 


The  PROBLEM  data  set  provides  a  complete  description  of  all  of  the 
nodes,  elements,  bodies,  etc.,  that  describe  the  problem  at  hand.  The 
order  of  the  data  records  is: 


Title  Record  Set 
PROBLEM 

(Problem  data  sets) 
SAO  Flag 

(SAO  data  as  needed) 
END  or  NEW 


The  order  of  the  problem  data  sets  is  not  rigidly  specified.  Use  the 
data  sets  as  needed  and  follow  logical  sequences.  For  example:  LINE 
requires  start  and  end  node  data  that  must  be  defined  by  NODE  and/or 
other  LINE  data.  The  data  set  descriptions  are  listed  alphabetically  by 
the  flag  record  following  the  title,  REST  and  PROB  descriptions.  A  summary 
of  flag  records  (four  characters)  follows; 
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Ti_t_l_e 

REST 

Restart  data 

PROB 

New  problem  definition 

BLOC 

Body  locations 

BODY 

Define  lumped  body  table 

ELEM 

Line  or  cable  element  definitions 

FLOW 

Flow-field  library  definitions 

FLU  I 

Fluid  media  definitions 

INVE 

Modify  component  inventory 

LIMI 

Limit  set  definitions 

LINE 

Generate  lines  of  nodes 

LLOC 

Limit  locations 

MATE 

Material  table  definitions 

NODE 

Node  point  definitions 

SHIP 

Ship  data  definitions 

STRU 

Strum  string  definitions 

TENS 

Initial  tension  input 

TFUN 

Time  function  library  definitions 

The  next  sections  present  detailed  instructions  for  inputting  this 
information  to  SEADYN,  including  expected  units  where  applicable 
(F  =  force,  L  =  length,  T  =  time). 
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6.1  Title  Record  Set 


Any  number  of  title  records  can  be  used  to  Desin  the  data  set.  At 
leas'-  one  is  required.  The  last  title  record  must  be  terminated  with  a 
record  terminator  (,  or  $).  These  record  terminator  characters  cannot 
appear  anywhere  else  in  the  title  set. 

Subsequent  requirements  for  a  title  record  set  (e.g.,  following  a  NJ^ 
flag)  must  be  limited  to  a  single  title  record  using  the  character  string 
deiimitOL's . 


21 


6.2  RESTART  -  Restart  data  record  -  Must  itnmediately  follow  title  records 
when  it  is  used. 


Word 


1 


2 

3 

4 

5 


Variable 

Name 


TYPE 


NTAPE 

NFILE 

IRST 

IDCHK 


Description 

Solution  type  being  restarted  (DEAD,  LIVE,  DYN,  NEW). 
(Job  will  abort  if  DEAD,  LIVE  or  DYN  is  given  and 
does  not  match  restart  data) 

File  code  number  for  the  data  file  (1,  2,  3  or  4) 

The  number  of  the  restart  record  to  be  read  from  the 
data  file.  (Default  is  the  last  one  written) 

Restart  file  flag.  (Same  function  as  in  SAO  data,  see 
Section  7.1.12  SAVE) 

Identification  check  flag 

1  -  Read  the  title  record  from  the  file  and 
compare  the  first  10  characters  with  the 
10  characters  given  in  CHKWRD  (Word  15); 
abort  if  they  do  not  match 

0  -  No  identification  check  is  made 


(The  data  words  beyond  this  point  are  not  used  if  Word  (1)  is  NEW) 


6 

NIPR 

New  value  for  number  of  steps  between  printing. 

If  this  is  >0  and  not  equal  to  the  value  given 
previously,  the  output  interval  will  be  changed 
by  the  ratio  NIPR  (new)/NlPR  (old) 

7 

IBG 

New  value  for  optional  debug  output  flag 

8 

DTT 

New  step  size 

9 

THAX 

New  maximum  time  (default  is  old  one) 

10 

GAMNEW  1 

11 

BETNEW 

New  integrator  parameters  (default  is  old  values) 

12 

ALPNEW  ] 

13 

DTOUT 

New  output  time  interval 
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Variable 


Word 

Name 

Description 

14 

DTRSRT 

New  restart  save 

time 

interval 

15 

NOTES; 

CHKWRD 

Label  check  word 

(see 

Word  5),  10  characters 

1.  The  data  file  referred  to  by  NTAPE  and  the  restart  record  referred 
to  by  NFILE  must  have  been  created  on  a  previous  SAO  using  the  SAVE 
option. 

2.  TYPE  (Word  (1))  is  used  to  signal  the  type  of  restart  to  be  executed. 
NEW  means  the  data  on  the  file  describes  the  beginning  point  of  a  new 
SAO  set,  and  an  SAO  flag  record  is  expected  on  the  next  input  record. 

The  DEAD,  LIVE,  and  DYN  types  signal  that  the  SAO  data  on  the  file  is 
the  same  as  the  SAO  form  indicated,  and  that  the  previous  SAO  is  to  be 
continued  using  the  indicated  parameter  changes.  Subsequent  SAO  sets 
can  follow  after  the  completion  of  the  restarted  SAO.  A  restart  to  con¬ 
tinue  a  TSSS  SAO  should  have  TYPE  =  "LIVE”. 


6.3  PROBLEM  -  Ne?w  problem  definition  -  Must  immediately  follow  title 
records  when  it  is  used. 

Variable 

Woi^  Name  Description 

1  NN  Number  of  nodes  in  model 

2  NE  Number  of  line  elements  in  model 

3  IDIR  Gravity  direction  code  (see  Note  1) 

0  -  no  gravity  loading 

±1  -  for  X  direction 

±2  -  for  Y  direction 

±3  -  for  'I.  direction 

4  ^JDYN  Dynamic  solution  option  flag 

^0  -  mass  and  acceleration  storage  will  be 
allocated 

-1  -  no  dynamics  terms  calculated  and  storage 
truncated  (DYN,  MODE,  FREQ  not  allowed) 

5  IPRO  Input  print  option  flag 

0  -  echo  all  input  data  with  interpretation 
-1  -  suppress  input  echo 

1  -  echo  input  plus  print  slopes  and  approximate 
transit  times  for  each  element 

6  INDRAG  Drag  model  override  flag  (see  Note  2) 

0  -  use  default  <lrag  model  regardless  of  flags 
on  individual  components 
^0  -  user-defined  drag  models  will  be  used  .is 
called  for 

7  GRAV  Gravitational  acceleration  (LT 

Defaults  to  32.17  ft/scc'  if  IDlR  #  0  and  GRAV 
not  given 

8  NSFILF,  Ghip  load  data  file  flag 

0  -  no  data  file 

>0  -  ship  load  data  for  /NSFII.E/  ships  are 

provided  as  rigid  format  input  in  this  run 
<0  -  data  file  available  with  /NSFILE/  ships  on  it 


NOTES: 

1.  It  is  assumed  that  the  global  axes  are  selected  such  that  i)iie  of' 
them  coincides  with  the  direction  of  gravity.  Tims,  IDIR  -  -2  moans  that 
gravity  acts  iu  the  negative  Y  direction  (i.e.,  +Y  is  up). 

2.  INDRAG  can  be  used  to  select  user-specified  drag  coefficients,  since 
it  is  passed  to  the  DRAGCO  subroutine  whenever  it  is  called.  Tf  it  is 
zero  it  overrides  all  other  drag  coefficient  specifications  and  the  pro¬ 
gram  will  use  the  built-in  or  default  values  described  in  Appendix  C. 
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0  . /*  BJiOC  -  Lumpotl/ I'igi  (i  l)0(ly  location  sp«?ci  f'i  cat  ions 


Word 

1 

2 

3 

4 
3 


Variable 

Name 


LBOUN  (I) 


l.BEG 

LEND 

LINC 

LIMSET  (3) 


Description 

Body  number  (see  BODY  table) 

Beginning  node  number  for  body  location(s) 
Ending  node  number  for  body  loc<ition(s) 
Node  number  increment 

Limit  set  number  (same  as  in  LLOC  record) 


NOTES: 


1.  Cylindrical  buoys  use  the  attached  elements  to  define  the  orientation 
of  the  cylindric.nl  axis.  The  first  two  connecting  elements  are  used. 

The  slopes  of  these  two  elements  will  be  averaged  to  get  the  orientation 
of  the  cylinder.  If  only  one  element  i:?  connected  to  the  body,  its 
direction  v^ill  be  used. 

2.  Mooring  buoys  (rigid  bodies)  require  a  minimum  of  four  jiodes  in  their 
definition.  Tlte  buoy  location  must  bo  a  twc.-node  pair  to  give  posLtioi) 
and  angle  (six  degree.S'o}'- I'rocdom) .  Two  additional  slave  nodes  mu.sl  be 
lie<i  to  the  buoy.  These  two  nodes  are  automatically  defined  from  the 
first  two  nodes  (in  numerical  order)  slaved  to  the  buoy  location  node. 
Recall  that  slave  nodes  need  not  have  linos  attached  to  them.  The 
location  node  for  a  rigid  body  is  the  first  of  tbe  two  node  pair  (the 
one  that  gives  the  position  coordinates).  When  less  than  three  lines 
are  attached  to  a  mooring  buoy,  the  roll  motion  about  the  line  between 
the  first  two  slave  nodivs  is  assumed  to  he  fi.xed  to  avoid  eipi.ilion 
singularity  problems. 

3.  Input  Cor  Words  (3)  and  (4)  is  not  required  if  only  one  body  is 
present . 

4.  No  more  than  10  lines  can  bo  connected  to  any  body  location  with 
limits  imposed  on  it  with  an  Ll.OC  or  BLOC  data  set. 

5.  The  present  version  of  SEADYN  does  not  allow  rigid  bodies  (using 
slave  and  angle  nodes)  in  the  time  domain  solutions  (l)YN). 
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6.!> 


YjitA 

1 

2 


4 


6 


7 


8 


9 

10 


BODY  '  Defines 
t  he  pretilem. 

Viiriiible 

Name 


I 

IDRB  (1) 


ADM  (I) 
DBU  (I) 
BLEN  (I) 


BAMC  (1) 

BWND  (T) 

BSCD  (I) 

BNOM  (1) 
MEDMB  (I) 


lup'ped  body  table  (A  catalog  of  bodies  available  to 
See  BLOC  record  to  identify  where  tli<^  bodies  are  used.) 


Description 

Body  table  number 

Drag  coefficient  number 

0  -  use  default  drag  (see  Appendix  C). 

“  use  DRAGCO  subroutine  unless  INDR  AO  =  0 

Body  buoyancy  (negative  means  weight)  (F) 

Body  diameter  (L) 

Body  length  (L) 

0  -  sphere  or  lump 
X)  *  in-line  cylindrical  body 

Added  mass  coefficient 

(Default  for  sphere  =  0.5,  for  cylinder  -  1.0) 

*> 

Wind  drag  coefficient  foi  .surface  buoy  (L‘') 


2 

Surface  current  drag  coefficient  (1,  ) 


«'i)  ‘s' 

Mass  moment  of  inertia  (Fl.T^) 

Fluid  medium  in  which  ADM(I)  is  defined. 
See  FLUID  record  (Default  =  1). 


NOTES: 

1.  lumped  bodies  can  be  spherical  or  au  in-line  cylinder.  These  bodies 
provide  only  loads/mass  to  the  nodes  where  they  are  located.  I'hey  <lo 
not  have  stiffness  terms,  nor  do  they  have  wave-induced  loads  in  the 
FREQ  subanalysis. 


2.  Bodies  used  as  mooring  buoys  cannot  be  cylindrical.  They  are 
treated  as  rigid  bodies  in  static  analyses  and  are  assumed  to  be  spherical 
bodies  with  the  water  line  at  the  e<|Mator  in  FREQ  subaiialyses  (see  Note  2 
on  BLOC  record). 


3.  Environmental  loads  on  the  body  are  determined  from  the  relative 
velocity  between  the  fluid  and  node,  the  dimensions  and  form  of  tlu' 
body,  and  the  drag  coefficients  (default  or  user  defined). 
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6.6  I'.LKMKNT  -  l.iuo  or  ili'fiiitt  ions 

Var i abJo 
Word  Naaio 

1  I  KlrnuMit  iiiiinhrr 

2  IT  (1,  I)  Nodo  nninbcr  lor  first  ond  of  elemc'iit 

3  IT  (2,  1)  Node  number  for  second  end  of  element 

U  Not  used  at  present  but  the  word  position  must  be 

accounted  for. 

5  MAT  (I)  Materiai  number 

6  KOMP  (I)  1'ypc  code 

0  -  element  can  not  support  compression 
1  -  element  can  support  compression 
-1  -  element  is  a  catenary  (see  Note  *.) 

7  KSKP  Generation  code  (see  Note  1) 

8  ISIGO  Initial  state  flag  (give  on  first  element  record 

only;  see  Note  2) 

0  -  equilibrium  initial  state  with  valid  strotclied 
lengths  and  tensions.  Word  (10)  not  used. 

1  -  compatible  initial  state  defining  unstretched 

lengths  with  dummy  tensions  given  in  Word  (9). 

2  •  gucs.sed  deformed  state.  DnsLrelched  lengths 

given  in  Word  (10).  Dummy  tensions  given  in 
Word  (9). 

9  SIG  (1)  Initial  tension  (see  Note'  3) 

10  nso  (I)  lliistretched  length  (retquired  for  catenary  element, 

see  Note  A) 

11  MEDIUM  (I)  Fluid  mediun^  flag  (the  elome>nl  is  assume'd  to  leiiiain 

in  the  original  medium  at  all  times).  See  FLUID 
record  (default  is  flui«l  1). 

NOTES: 

1.  Elements  must  be  in  ascending  element  number  order,  and  the  first  and 
last  ones  must  be  input.  Any  omitted  elements  will  be  generated  using 
the  KSKP  value  ou  the  element  card  omitt«>d  clniieuls. 

Elements  are  incrementeil  by  1;  node  nvimbers  are  incremented  by  KSKP.  ITu' 
element  preceding  the  omitted  elements  is  used  to  begin  the  increments. 
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For  example:  if  element  1,  defined  by  nodes  6  and  12,  preceded  the 
omitted  elements  and  KSKP  =  4,  then  the  next  element  generated  would  be 
element  2  defined  by  nodes  10  and  16.  The  generated  elements  will  be 
assumed  to  have  the  same  properties  as  given  by  the  element  previous  to 
the  omitted  ones. 

2.  The  initial  state  flag  (ISIGO)  has  three  basic  options: 

ISIGO  =  0:  The  preloads  are  presumed  to  represent  an  equilibrium 

state.  The  element  lengths  obtained  from  the  nodal 
positions  are  used  with  these  preloads  and  the 
constitutive  relations  to  calculate  the  unstretched 
lengths.  DSO(I),  Word  (10),  is  not  required. 

ISIGO  =  1:  The  configuration  represented  by  the  nodal  positions 

is  an  unstretched  compatible  arrangement  of  cables. 

The  preloads  (SIG(I))  are  merely  estimates.  The 
unstretched  lengths  are  obtained  from  the  nodal 
positions . 

ISIGO  =  2:  The  configuration  represented  by  the  nodal  positions 

is  a  compatible  estimate  of  the  deformed  configuration. 
The  values  of  the  unstretched  lengths  are  given  on  the 
element  records.  Element  loads  are  calculated  from 
the  unstretched  lengths,  stretched  lengths  represented 
by  the  nodal  positions,  and  the  material  constitutive 
relations . 

3.  When  a  string  of  elements  lies  along  a  catenary  curve  that  has  been 
generated  by  a  node  generation  record  (LINE  record),  the  initial  tension 
for  each  element  can  be  obtained  internally  from  the  catenary  generation 
using  the  TENS  record.  The  initial  tension  will  be  taken  to  be  the 
catenary  tension  half-way  between  the  two  nodes  defining  the  element. 

SIG  (I)  is  not  input  if  the  TENS  record  is  used. 

4.  A  catenary  element  is  signaled  by  KOMP(I)  <  0.  The  unstretched 
length  input  is  required  in  this  case.  The  catenary  element  uses  the 
positions  of  the  two  nodes  and  the  unstretched  length  to  estimate  stiff¬ 
nesses  and  internal  reactions.  This  is  done  by  iteratively  solving  the 
classical  catenary  equations  with  adjustments  for  line  stretch.  Only 
the  weight  of  the  catenary  is  used  in  this  iteration  (current  loading  is 
neglected).  The  lowest  node  on  the  element  is  assumed  to  be  at  a 
bottom-limited  condition.  The  line  is  assumed  to  lie  on  the  bottom 
rather  than  to  sag  when  the  line  length  is  sufficient.  Note  that  the 
catenary  line  generation  of  the  LINE  record  (Section  6.11)  is  a  node 
generator  and  has  no  function  relative  to  defining  catenary  elements. 

The  unstretched  length  may  be  obtained  from  the  LINE  input  data  by 
entering  a  minus  times  the  catenary  number  from  the  LINE  record  (see 
Note  4  of  Section  6.11)  for  Word  10.  Note:  The  present  algorithm  fails 
when  the  top  node  of  a  catenary  element  is  on  the  bottom. 
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6.7  FLOW  -  Flow-field  library  definitions 


Word 

1 

2 

3 

4 


Variable 
Name  _ 

I 

I FLOOD  (I) 
FLPAR  (1,  I) 
FLPAR  (2,  I) 


Description 

Flow-field  number  (10  max) 

Flow-field  code  (see  Note  1) 

Optional  flow-field  parameters  (see  Notes) 


12  FLPAR  (10,  I) 


NOTES: 

1.  Flow-field  codes  select  a  particular  form  of  built-in  function 
or  signal  a  call  to  the  subroutine  CURUSR.  Positive  field  codes 
refer  to  built-in  functions,  while  negative  codes  refer  to  user 
subroutines.  A  request  for  a  built-in  function  beyond  those  currently 
defined  will  cause  an  abort  of  the  run.  The  parameters  are  used  as 
appropriate  for  the  built-in  functions,  and  they  are  provided  in  the 
calling  sequence  to  the  user  subroutines. 


2.  The  user  subroutine  calling  sequence  is: 

CALL  CURUSR  (T,  N,  NN,  X,  V,  FLPAR  (1,  I)) 


where  T 
N 
NN 
X 
V 


current  time 
ABS  (IFLCOD  (D) 
number  of  nodes 
nodal  position  vector 
nodal  flow  velocity  vector 


The  nodal  vector  storage  format  is  X  (I,  J)  and  V  (I,  J),  where 
I  =  1,  2,  3  for  X,  Y,  and  Z  components,  respectively;  and  J  =  node  number. 


3.  The  only  built-in  flow  field  at  present  is: 
IFLCOD 

1  -  uniform  velocity  field 
FLPAR 

1  -  X  component  of  flow  velocity  (LT_J) 

2  -  Y  component  of  flow  velocity  (LT^f) 

3  -  Z  component  of  flow  velocity  (LT”^) 
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6  8  FI.UID  -  Fluid  media  definitions  (used  only  if  IDIR  ^  0  on  PROB  record). 
Must  precede  MATE  and  BODY  records. 


Word 

Vari able 
Name 

Description 

1 

FDEPTH  (I) 

Coordinate  at  fluid  surface 

2 

-- 

Property  code 

0  -  read  properties  from 

Words  (3)  and  (4) 

1  -  1.30  X  10‘^  ft^/soc. 

64.0  Ib/ft^  (seawater) 

2  -  1.68  X  lO'^  ft^/sec, 

0.0765  Ib/ft^  (airl 

3 

FVISC  (I) 

2  -1 

Kinematic  viscosity  (L  T  ) 

4 

FGAM  (I) 

specific  weight  (FL 

NOTES; 

1.  Provide  one  record  for  each  fluid  in  ascending  order.  The  first 
fluid  given  is  assumed  to  extend  infinitely  below  its  surface. 

2.  The  FDEPTH  for  the  highest  f.»uid  need  not  be  given.  It  is  assumed 
to  extend  infinitely  above  the  previous  fluid. 

3.  Pressure  effects  from  the  accumulated  depths  are  ignored. 

4.  The  present  version  of  the  program  only  allows  two  fluids. 
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6.9  INVKNTORY  -  Units  conversion  for  component  inventory  (see  Appendix  D 
for  the  inventory  contents). 

Variable 

Word  Name  Descript  i_o  n 

1  FRCVY  Conversion  factor  for  weights  and  strengths 

(Default  =  1.0) 

2  FLNVY  Conversion  factor  for  buoy  diameters  and  lengths 

(Default  =  1.0) 

3  FDIVY  Convcfrsion  factor  for  line  diameters  (Default  =  1.0) 

A  INVPNT  Inventory  Print  Flag 

0  -  do  not  print  inventory 
1  -  print  all  inventory  entries 

NOTES : 

1.  The  value  obtained  from  the  inventory  is  multiplied  by  the  corresponding 
factor  to  get  it  into  the  units  implied  by  the  problem  input  data. 

2.  Inventory  data  will  not  be  used  if  this  record  is  not  input. 


6.10  LIMIT  -  Limit  set  definitions  (used  only  if  IDIR  t  0  on  PROB  record) 
Var i ab  le 

Word  Name  Description 


1  LIMIT  (1)  Limit  set  number 

2  CORDLM  (I)  Vertical  coordinate  of  limit  (L) 

3  TOLIM  (I)  Limit  tolerance  (L) 

A  RELFAC  (I)  Release  factor  (default  =  1.001) 

5  JANCR  (1)  Fixity  code 

0  -  buoyant  limit  (surface  bound)  with  limit 
imposed  only  on  vertical  component 
1  -  weight  limit  (bottom  bound)  with  limit 
imposed  only  on  vertical  component 
3  -  weight  limit  with  all  three  components  held 
when  limit  is  imposed 

NOTES : 

1,  Repeat  records  for  as  many  limits  as  needed  up  to  50. 

2.  This  data  set  defines  a  table  of  limit  conditions.  These  are  assigned 
to  specific  nodes  using  LLOC,  BLOC,  or  possibly  LINE  records  (see  Note  3 
of  LINE  record). 


3.  RELFAC  gives  the  ratio  of  the  sum  of  the  vertical  components  of  the 
line  tensions  to  the  external  loads  at  the  limited  node  that  must  be  ex¬ 
ceeded  to  release  a  limit  in  a  transient  (DYD)  analysis.  External  loads 
include  weight/buoyancy  of  the  lines,  bodies  at  the  node,  and  any  current 
and  point  loads.  The  release  factor  is  always  1,0  in  any  static  analysis. 


6.11  LINE  -  Generates  lines  of  nodes 
Variable 

Woj^  Name _  Description 

1  NUMN  Number  of  nodes  to  be  generated 

2  NOLD  Beginning  reference  node  (lowest  node  on  a  catenary) 

3  I  Ending  reference  node  (highest  node  on  a  catenary) 

4  KSKP  Mode  numbering  increment  (default  =1) 

5  NBC  Boundary  condition  flag 

0  -  do  not  copy  constraint  codes.  Set  all 
constraints  to  zero. 

1  -  copy  constraint  codes  from  node  MOLD 

2  -  copy  constraint  codes  from  node  I 

6  NFIRST  Number  of  first  generated  node  (default  =  NOLD  +  KSKP) 

7  NCAT  Generation  code 

0  -  straight  line 

1  -  catenary  with  sag  laid  flat  at  limit 

(default  =  vertical  coordinate  of  NOLD) 

2  -  catenary  with  unconstrained  sag 

3  -  catenary  as  in  1  except  no  nodes  placed 

on  bottom 

8  GCATW  Catenary  weight  per  unit  length  (FL  ^) 

9  GCATH  Catenary  horizontal  component  of  tension  (F) 

10  LIMCAT  Limit  set  number  for  NCAT  =  1  (see  Word  (1)  of 

LIMIT  record)  may  be  used  in  place  of  the  LLOC 
data  for  nodes  generated  on  the  bottom  by  the 
NCAT  =  1  option  (see  Note  3) 

NOTES: 

1.  Both  the  beginning  and  ending  nodes  must  have  been  previously  defined 
by  NODE  or  LINE  records.  The  beginning  and  ending  nodes  need  have  no  num¬ 
erical  relation  to  the  generated  node  numbers;  only  the  spatial  relation 
of  the  generated  nodes  is  significant. 

2.  The  value  of  KSKP  can  be  positive  or  negative.  Incrementing  starts 
from  NFIRST. 
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3.  All  generated  nodes  are  evenly  spaced  on  the  generated  straight 
line.  Nodes  generated  on  a  catenary  are  placed  to  give  uniform  angle 
change,  except  for  NCAT  =  1  with  line  on  bottom.  If  sag  below  the 
bottom  is  detected  with  NCAT=1,  the  tangent  point  is  calculated.  If  the 
distance  between  NOLD  and  the  tangent  point  is  greater  than  0.75  times 
the  uniform  spacing  distance,  a  node  is  placed  at  the  tangent  point.  If 
the  distance  is  greater  than  1.50  times  the  uniform  spacing  distance, 
NFIRST  is  placed  half-way  between  the  tangent  and  NOLD.  The  node  at  the 
tangent  is  then  NFIRST  +  KSKP,  otherwise  the  tangent  is  NFIRST.  The 
remaining  generated  nodes  are  placed  on  a  new  uniform  angle  spacing  from 
the  tangent  point  to  the  last  node.  All  nodes  placed  on  the  bottom  are 
given  the  limit  set  number  LIMCAT.  If  LIMCAT  is  not  given,  the  LLOC  (or 
BLOC)  records  should  be  used  to  set  limits. 

A.  The  catenary  options  provided  here  are  simply  for  the  purpose  of 
generating  nodes.  They  have  no  direct  function  relative  to  the  generation 
of  catenary  line  elements  (see  Section  6.6).  Information  needed  to 
calculate  the  unstretched  length  of  a  catenary  element  can  be  passed 
through  this  node  generation  input.  This  is  done  by  using  NCAT  =  3  and 
following  Note  A  of  Section  6.6.  This  information  pass-through  can  be 
done  even  when  no  nodes  are  to  be  generated  by  setting  NUMN  =  0.  This 
particular  situation  arises  when  a  single  element  leg  is  desired  and  the 
horizontal  component  of  pretension  is  known  instead  of  the  line  length. 
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6.12  LLOC  -  Limit  location  specification 
Variable 


Word 

Name 

Description 

1 

LIHSET 

Limit  set  number  (see 

Word  (1)  on  LIMIT  record) 

2 

LBEG 

First  node  where  limit 

is  applied 

3 

LEND 

Last  node  where  limit 

is  applied 

4 

LINC 

Node  number  increment 

NOTES: 

1.  The  maximum  number  of  limit  locations  is  50. 

2.  No  more  than  10  line  elements  can  be  connected  to  a  limited  node. 

3.  Generation  of  nodes  on  a  catenary  can  cause  generation  of  up  to  two 
limited  nodes  per  catenary.  This  situation  is  described  in  Note  3  of 
the  LINE  record . 

6.  Body  location  data  can  also  be  used  to  locate  node  limits  (see 
BLOC  record). 


6.13 

MATERIAL  -  Material  table  definitions 

Word 

Variable 

Name 

Description 

I 

I 

Material  number 

2 

IDKG  (I) 

Drag  coefficient  code 

0  -  use  default  drag  coefficients  (see 

Appendix  C) 

4  0  and  INDRAG  ^  0  -  use  DRAGCO  subroutine 

3 

DIAM  (I) 

Cable  or  line  diameter  (L).  Used  only  for  fluid 
load  computation. 

4 

G3  (1) 

Weight  per  unit  length  (negative  means  buovanl) 

(FL 

5 

HF,D  (I) 

Reference  medium  code.  See  FLUID  record 
(default  is  fluid  1). 

6 

CAMC  (I) 

Added  mass  coefficient  (default  =  1.0) 

V 

TENULT  (I) 

Ultimate  tension  capability  (F)  (see  Note  3  for  use) 

8 

ilE  (I) 

Option  flag  for  tension  (T)/strain  (c)  properties 

0  -  use  exponent  form  T  =  af. 
n  -  n  points  in  tabular  form  (maximum  n  is  20) 
(see  Notes  1  and  2) 


9 

n 

(1. 

I) 

First  tension  in  table  or  a 

(F) 

10 

STR 

(1, 

I) 

First  strain  in  table  or  b 

11 

TTD 

(I) 

Damping  parameter  CAj  (FT) 

(see  Note  A) 

12 

HK 

(I) 

Damping  parameter  EA^  (F) 

(see  Note  A) 

13 

TT 

(2, 

I) 

Second  tension  in  table  (F) 

lA 

STR 

(2, 

I) 

Second  strain  in  table 

(repeat  pairs  for  all  table 

points) 
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NOTES : 


1.  When  point  input  is  used  and  data  are  required  outside  of  the  given 
values,  the  slope  of  the  last  (first  if  only  one)  segment  is  assumed  to 
be  continued. 

2.  Input  must  be  in  increasing  tension  order.  Thus,  materials  with 
compression  stiffness  must  begin  by  listing  the  largest  compressive  load 
first  and  progressing  to  the  largest  tensile  load. 

3.  The  ultimate  tension  input  is  used  only  by  the  frequency  domain 
dynamic  solution.  When  TENULT  (I)  is  greater  than  zero,  the  random  esti¬ 
mates  for  tensions  will  be  factored  and  compared  to  the  ultimate  tension. 

U,  There  are  two  material  damping  models  provided: 

KELVIN  Model 


EA  is  provided  in  the 
tension/strain  data 
CAj  is  the  dashpot  coefficient 

NOAA-REID  Model  CAj 

I— Ih-AA/V 


KA 

- - ' 

EAj  is  an  additional  stiffness  parameter.  Additional  information 
regarding  the  applicability  and  use  of  these  models  can  be  found  in 
Reference  13. 


_ ,c;/ 

HD- 

AAAr 
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6.14  NODE  -  Node  point  definitions 


Variable 
Word  Name 

1  1 

2  KSKP 

3  XO  (3*1-2) 

4  XO  (3*1-1) 

5  XO  (3*1) 

6  NODFIX  (3*1-2) 

7  NODFIX  (3*1-1) 

8  NODFIX  (3*1) 

9  DFLAG 


Description 
Node  number 
Node  generation 

•v 

X 

Y  ^Global  coordinates  of  node  I  (L  for  position, 
dimensionless  or  degrees  for  angles) 


'Displacement  constraint  codes  (see  Note  2) 


Angular  coordinate  flag  --  "DEGREE"  -  means  the 
nodal  coordinates  of  angle  nodes  are  input  in 
degrees.  Otherwise  angles  arc  in  radians 
(see  Note  1). 


NOTES : 


1.  Each  node  is  assumed  to  require  three  components  to  describe  its 
position  relative  to  a  global  right-handed  cartesian  system.  Rigid 
bodies  requiring  six  degrees-of -freedom  (ships  and  mooring  buoys)  use 
two  consecutive  nodes.  The  first  gives  position  and  the  second  gives 
angular  coordinates.  All  angular  computations  are  doitc  in  radians,  and 
angular  output  is  in  radians.  Input  of  initial  angular  positions  can  be 
either  in  radians  or  degrees  as  indicated  by  DFLAG.  No  equations  are 
available  for  rigid  bodies  in  the  time  domain  solutions  (DYN). 

2.  Constraint  codes  define  the  active  and  constrained  components  of 
movement . 

0  -  no  constraint  (free) 

1  -  no  movement  (fixed),  subject  to  release  by  limit  set  conditions 
-N  -  means  the  node  is  slaved  to  node  N 

2  -  identifies  a  component  reserved  for  payout  that  will  be 

unconstrained  (free)  when  activated.  (A  flag  of  1  remains  in 
effect  on  a  node  activated  by  payout.)  Cannot  be  freed  by  a 
limit  condition. 

3  -  identifies  a  component  that  is  fixed  and  cannot  be  freed  by  a 

limit  set  condition. 
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3.  Slave  nodes  must  be  numbered  last.  This  means  the  node  numbers  of 
slave  nodes  must  be  larger  than  any  of  the  active  and  roaster  nodes.  The 
input  routine  counts  the  slave  nodes  and  deducts  that  number  from  the 
total  nodes  (Nh)  to  get  the  number  of  active  nodes.  Slave  node.s  need 
not  have  lines  connected  to  them.  Slave  nodes  cannot  be  used  in  time 
domain  solutions  (DYN) .  Master  nodes  are  assumed  to  be  a  node  pair 
(see  Note  1)  since  the  slave  nodes  implies  angular  responses. 

4.  A  node  with  any  or  all  constraint  codes  set  to  1 ,  2,  or  3  is  still 
considered  an  active  node  even  though  it  cannot  move.  It  is  assumed 
that  any  of  the  constrained  components  can  be  modified  by  FREE,  MOVE,  or 
PAYOUT.  LIMIT  checks  will  affect  only  components  with  a  constraint  code 
of  1  or  0. 

5.  Input  of  nodes  need  not  be  in  numerical  order.  Nodes  can  be  entered 
more  than  once,  and  the  last  entry  will  be  the  one  used  in  the  solution. 

Omitted  nodes  can  be  generated  using  a  straight  line  with  uniform  spacing 
with  the  aid  of  the  KSKP  parameter.  If  KSKP  ^  0  on  a  node  record,  then  that 
node  is  designated  as  the  last  node  on  the  line,  and  the  one  input  preceding 
it  is  designated  as  the  first  node  on  the  line.  Nodes  are  generated 
evenly  along  the  line  between  these  two  points  with  node  numbers  incre¬ 
mented  by  KSKP  from  the  first  to  the  last  node.  The  difference  between 

the  node  numbers  at  the  ends  of  the  line  must  be  an  integer  multiple  of 
KSKP,  and  KSKP  cannot  be  negative.  The  generated  nodes  will  have  the 
same  constraint  codes  as  the  first  node  on  the  line.  More  general  node 
generation  schemes  are  available  using  LINE. 

6.  All  of  the  NN  nodes  must  be  accounted  for  in  the  combined  specifica¬ 
tions  for  nodes  and  generation  schemes. 

7.  Slave  nodes  need  not  have  elements  connected  to  them.  This  allows 

the  investigation  of  the  response  of  specific  points  on  a  ship  at  locations 
other  than  the  attachments  of  mooring  and  working  lines. 

8.  All  points  of  attachment  of  lines  to  rigid  bodies,  such  as  ships  or  mooring 
buoys,  must  be  slaved  to  the  primary  node  defining  the  body.  The  program 
imposes  no  limit  to  the  number  or  relative  locations  of  these  attachments, 
except  that  the  frequency  domain  solution  for  mooring  buoys  uses  the 

first  two  attachments  to  define  the  local  coordinate  system  for  the 
motion  equations.  These  two  attachments  are  also  used  for  reference  in 
the  CHEK  option  (see  Note  2  on  BLOC  record). 

9.  The  initial  orientation  of  a  ship  relative  to  the  global  coordinate 
system  must  be  input.  The  angle  required  in  this  definition  can  be  seen 
in  the  following  sketch.  (The  global  axes  are  X,  Y,  Z  and  the  local 
ships  axes  are  1,  2,  3,  representing  aft,  starboard,  and  up,  respectively.) 


6.n 

yilP  Sliip  (l.'ita 

d('f  iiiit  inup 

Word 

Va  riablo 

Name 

Description 

] 

ISHIP  (1) 

Node  where  ship  is  located 

2 

LSHP  (I) 

Load  function  option 

-1  -  Search  file  for  equivalent  ship 

0  -  use  app^lytical  load  functions 
n>0  “  use  n  ship  from  loading  file 

3 

CPROP  (1) 

Propeller  resistance  coefficient  (default  =  l.'Kl) 

4 

OR  (1) 

Longitudinal  resistance  coefficient  for  hull  to 
be  used  only  for  analytical  functions 
(default  =  program  calculates  one,  see  Note  5) 

5 

cs  (I) 

Hull  wetted  surface  coefficient  to  be  used 
only  for  analytical  functions  (default  =  2.70) 

6 

CMS 

Amidships  coefficient  to  be  used  to  calculate 
for  analytical  functions  (default  =  0.98) 

7 

SLT  (1) 

Total  length  of  ship  (L) 

8 

SAH  (I) 

End  projected  wind  area  (L^) 

9 

SAS  (1) 

Side  projected  wind  area  (1,^) 

10 

SLWL  (I) 

Watei—iine  leagtli  (L) 

11 

SBEAH  (I) 

Beam  at  midships  (1.) 

12 

SDRFT  (1) 

Draft  at  midships  (L) 

13 

SUSP  (I) 

Volume  displacement  (1.^) 

U 

APROP  (I) 

Propeller  projected  area  (L^) 

15 

FSFRW  (I) 

Load  table  wind  force  conversion  factor 
(default  =  1.0) 

16 

FSFRC  (I) 

Load  table  current  force  conversion  factor 
(default  =  1.0) 

17 

FSLEN  (1) 

Load  tabic  length  conversion  factor 
(default  =  1.0) 

18 

FSVEL  (1) 

Load  table  velocity  conversion  factor 
(default  =  1.0) 
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Variable 


Word 

N;imc 

Description 

19 

SHIPK  (1,1) 

2  2 

Heave  restoring  coefficient  (default  “  10  ) 

20 

SlllPK  (2,1) 

7‘> 

Koll  rostorittg  cooflicirnl  ■•=  10  ‘  ) 

2) 

SHTPK  (3,1) 

Pitch  restoring  coefficient  (default  =  10^^) 

22 

SHIPK  (4,1) 

Heave/pitch  restoring  coefficient  (default  -  0) 

23 

WDEPTH  (T) 

Water  depth  at  .ship  location  (default  =  10  times 
draft)  (see  Note  6) 

NOTES; 

1.  Linearized  ship 

restoring  coefficients  can  be  input  if  desired. 

Otherwise  the  ship  will  be  ass>m»od  fixed  in  heave,  roll,  and  pitch 
during  static  analyses.  During  frequency  domain  analyses,  the  restoring 
matrix  is  obtained  from  the  ship  motion  file. 

2.  Words  (7)  through  (14)  may  be  omitted  if  LSHP  (I)  =  n  and  the  moored 
ship  is  the  same  as  on  file  n  of  the  ship  loading  file  (i.e.,  no  similarity 
scaling  required). 

The  conversions  factors.  Words  (15)  through  (18),  multiply  the  values 
of  force,  length,  and  velocity  from  the  ship  loading  file  to  get  values 
consistent  with  the  units  implied  by  the  rest  of  the  input  data. 

4,  The  local  coordinate  system  for  ships  asramies  the  1,  2,  and  direc¬ 
tions  are  .aft,  starboard,  and  up,  respectively.  This  is  consistent  with 
the  loading  conventions  of  Appendixes  B,  K,  and  0. 

5.  When  analytical  loading  functions  are  used,  the  longitudinal  resistance 
coefficient  can  he  input  or  calculated  using  a  table  look  up.  The  total 
coefficient,  CR  (I),  consists  of  three  parts: 


=  residuary  coefficient 
Cj.  =  frictional  resistance  coefficient 
=  fouling  resistance  coefficient 
The  value  for  C^.  is  always  calculated  in  the  program  using: 


0.4.S6 


(loflloKe^ 


05 


1700 

R 


R  >  5  X  10 
e 


where  is  Reynold's  Number  based  on  the  longitudinal  coaiponent  of  the 
velocity;  or 


-  0.002  R  <  5  X  10^ 
f  e 
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6.16  STRUM  -  Strum  string  definitions 


Variable 

Word  Name _  Description 

1  ISTRNG  (I)  Number  of  elements  in  string  (max  20) 


2  KSTRNG  (1,  I)  List  of  string  elements 

3  KSTRNG  (2,  I)  Must  be  listed  in  sequence  from  one  end  of  the 

string  to  the  other.  It  does  not  matter  which 
end  is  listed  first. 


31  KSTRNS  (30,1) 

NOTES: 

1.  A  string  consists  of  adjacent  (contiguous)  elements  to  be  considered 
in  estimating  drag  amplification  factors  for  strumming  analyses.  An 
element  can  be  contained  in  more  than  one  string.  The  drag  amplification 
factor  for  each  element  listed  in  a  string  is  the  largest  value  it  has 

in  any  string. 

2.  The  maximum  number  of  strings  is  30.  The  maximum  number  of  elements 
per  string  is  20. 

3.  Drag  amplification  estimates  are  made  on  LIVE,  DYN,  and  TSSS  subanalyses 
with  nonzero  relative  fluid  velocities  whenever  these  strum  strings  are 
defined.  The  procedure  assumes  that  the  strings  are  supported  at  each  end 
and  an  auxiliary  vibration  mode  analysis  is  conducted  on  each  string  to 
estimate  its  possible  involvement  with  vortex  shedding-induced  strumming. 

The  general  approach  of  Skop,  Griffin  and  Remberg  (Ref  lA)  is  used.  The 
occurrence  of  the  drag  amplification  estimates  is  controlled  by  the 
changes  in  relative  velocity  and  the  CEPS  parameter  input  on  the  SOLU 
record  (Section  7.1.1A)  associated  with  the  SAO  data  set. 
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6,17  TENSION  -  Initial  tension  input 


Variable 


Word 

Name 

Description 

1 

LBEG 

Beginning  element  number 

2 

LEND 

Ending  element  number 

3 

LINC 

Element  increment 

4 

LCODK 

Option  code 

0  -  assigned  initial  tension,  the  value  in 
Word  (5) 

N  -  estimate  initial  tension  from  the  N 
catenary  defined  on  the  LINE  records 
(see  Note  2) 

.SiC  Tension  to  be  assigned  to  eJeinerits  (F). 

NOTES: 

1,  Catenary  line  tensions  (LCODE  =  N)  for  the  catenary  are  calculated 
for  eacli  e j.eaient  based  oii  the  position  on  the  catenary  at  the  midpoint  of 
rich  element  (see  ITNF,  record  and  ELEMENT  record,  Note  3). 

2.  The  TENSION  record  must  not  precede  the  LINE  record  which  generates 
the  catenary  this  data  record  refers  to. 
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6,18  TFUNCTION  -  Time-function  library  definitions 


Variable 

Word  Name 

1  1 

2  ITFCOD  (1) 

3  TPARM  (1,  I) 

6  TPARM  (2,  I) 


Description 

Time-function  number  (20  max) 
Time-function-type  code  (see  Notes) 
Optional  time-function  parameters 


22  TPARM  (20,  I) 


NOTES: 

1  Tirae-function-type  codes  select  a  particular  form  of  built-in  function 
or  sienal  a  call  to  the  subroutine  TFNUSR.  Positive-function  cooes 
refer  to  built-in  functions,  while  negative  codes  refer  to  the  i^ser 
subroutine.  A  request  for  a  built-in  function  beyond  those  currently 
defined  will  cause  an  abort  of  the  run.  The  parameters  are  used  as 
appropriate  for  the  built-in  functions,  and  they  are  provided  in  the 
calling  sequence  to  the  user  subroutine. 

2.  The  user  subroutine  calling  sequence  is: 

CALL  TFNUSR  (T,  F,  N,  TPARM  (1,  T)) 


where  T 
F 
N 


current  time 

the  returned  value  of  the  time  function 
.ABS  (ITFCOD  (1)) 


(See  Appendix  H.) 
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The  built-in  functions  presently  available: 
ITFCOD 

1  Ramp  Bui Id -up /Decay 


F 


2  Shifted  Sine 

F  =  A  Sin  (u)  (T-T^))  +  B 


3  Ramped  Sine 


Tr 


4  Decayed  Sine 


TPARM  INPUT 
12  3  4 


T  F  F 

1  ‘2  '2 


A  T„  T  B 
0  s 


A  T-  T  B 
0  s 


ITFCOD 


Tabular  Input  (max  of  9  points  plus  end  values) 


F 


Function  value  for  times  less  than 

First  time  point 

First  time  function  value 

etc. 


18  Tg 


20  Function  value  for  times  greater  than 

last  time 

ITFCOD 

A  Random  File  Input 
TPARM 


1 

NOTAPES 

Number  of  excitations  (5  max) 

2 

TTZERO 

Starting  time  on  excitation  tapes 

3 

SCAFACT  (1) 

Multiplier 

4 

(2) 

5 

(3) 

6 

(A) 

7 

(5) 

8 

SHIPVEL  (1) 

Velocities  of  excitation  points 

9 

(2) 

10 

(3) 

11 

(4) 

12 

(.^) 
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7.0  SUBANALYSIS  OPTION  (SAO)  DATA 


Each  subanalysis  is  headed  by  an  option  flag  record; 

DEAD  static  analysis  with  gravity  and  point  loads 

LIVE  static  analysis  with  arbitrary  combined  loads 

TSSS  time-sequenced  static  solutions  (approximate  dynamic 

analy.sis  using  LIVE  that  allows  moved  nodes,  payout, 
and  other  time-varying  loads  but  neglects  inertia 
effects) 

DYN  transient  nonlinear  dynamic  analysis  (time  domain) 

FREQ  frequency  domain  dynamic  solution  (response  spectrum 

analysis)  for  wave  excitation 

MODE  determination  of  natural  frequencies  and  mode  shapes 

for  current  position 

CHEK  component  adequacy  checks 

END  run  termination 

NF.W  define  new  problem  --  title  card  read  next 

An  SAO  flag  record  has  the  flag  name  as  Word  (1)  with  no  other  data 
in  the  record.  SAO  data  are  grouped  following  the  SAO  flag  to  identify 
solution  characteristics,  loading,  boundary  conditionf; ,  and  outjuit 
requests.  Appropriate  default  values  are  assumed  when  no  data  arc  given 
Unless  explicity  .stated,  the  data  records  have  no  required  order,  since 
each  data  set  is  idontified  by  a  keyword. 
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7.1  DEAD,  LIVE,  TSSS,  DYN  Data  Set 


The  data  set  description  for  DEAD,  LIVE,  TSSS,  and  DYN  are  given  in 
alphabetical  order  of  the  keyword: 


Keywo  rd 


CURR  flow-field  specification 

FIX  applies  temporary  fixed  conditions  to  nodes 

FREE  releases  previously  imposed  fixity 

TMPA  impacting  body  input 

INIT  describes  dynamic  initial  conditions 

KEEP  retains  the  data  defined  in  the  preceding  SAG 

for  this  SAG 

LOAD  specifies  load  conditions 

LVAR  specifies  load  variation  codes 

MOVE  specifies  displacement/velocity/acceleration  at  nodes 

OUTP  selects  output  data 

PAYG  defines  payout/reel-in  of  lines 

SAVE  defines  restart  file  save  intervals 

SBUO  specifies  solution  option  characteristics 

SOLD  solution  option  characteristics 

STEP  solution  step  size  data 

SURF  surface  current  data 

TIME  time  step  data 

WIND  surface  wind  data 
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7.1.1 

CURRENT  -  Flow- 

field  specification 

Word 

Variable 

Name 

Description 

1 

OPTION 

"CURR" 

2 

NFLUID 

Flow-field  number  from  library  defined  in 
FLOW  record  (see  Notes) 

3 

CURMUL 

Flow-field  scale  factor  (multiplier  for  the 
values  obtained  from  the  flow  library; 
default  is  1.0) 

4 

NFLVRY 

Flow-field  variation  code 

For  LIVE  SAP 

1  -  increase  flow  field  incrementally 
0  -  hold  flow  field  at  full  amplitude 
“1  -  decrease  flow  field  incrementally 

For  DYN  and  TSSS  SAP 

0  -  hold  flow  field  amplitude  constant 
or  use  the  time  variation  implied 
by  th|j^ClJRUSR  subroutine 
>0  -  use  n  "  time  variation  function 
from  TFUN  record 

NOTES: 

1.  When  NFLUID  ^  0,  the  flow  field  will  be  evaluated  each  time  the 

fluid  loads  are  calculated.  The  flow  velocity  obtained  from  the  requested 
flow  field  (see  FLOW  library,  Section  6.7)  will  be  multiplied  by  CURMUL 
and  the  appropriate  time  variation  code  for  DYN  and  TSSS.  Note  that  for 
LIVE  subanalyses,  the  variation  code  multiplies  the  fluid  loads  instead 
of  the  fluid  velocity. 

2.  Wlien  NFLUID  <  0,  the  flow  field  will  be  evaluated  only  once  during 
the  SAO. 
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7.1.2  FIX  -  Applies  fixed  conditions  to  nodes 


Variable 
Word  Name 

1  OPTION 

2  I CODE 

3 


4 

(etc.) 


NOTES: 


Description 

"FIX” 

Constraint  code  (default  =  1,  see  Note  2) 

List  of  node  component  codes 

CODE  =  NODE*10+I 

where  I  =  1,  2,  or  3 

for  X.  Y  or  Z  direct  ion >  respecti  vely 

Maximum  list  length  is  98  (use  additional  FIX 
records  to  get  more) 


1.  FIX  causes  the  selected  code  to  be  placed  in  the  constraint  code  array 
(NODFIX)  for  the  indicated  component.  The  component  will  remain  fixed 
until  released  by  a  FREE  record  or  a  predefined  LIMIT  record. 

2.  'Ilowable  constraint  codes: 

1  -  component  fixed  unless  limit  condition  override.s. 

2  -  component  fixed  and  reserved  for  payout.  It  will  become  Iroe 

when  activated  by  a  payout  mitosis'’'  (see  PAYO  record, 

Section  7.1.11).  It  cannot  be  released  by  a  limit  condition. 

3  -  component  fixed  and  cannot  be  released  by  a  limit  condition. 

3.  The  MOVE  record  will  override  any  constraint  code. 

CAUTION:  Do  not  FIX  any  components  of  a  slave  node. 


*Mitosis  refers  to  the  halving  of  an  element  into  two  new  separate 
elements;  the  original  element  length  is  replaced  by  two  elements 
(see  Reference  1  for  further  details). 
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7.1.3  FREE  -  Removes  fixed  conditions  for  nodes 


Variable 
Word  Name 


Description 


1  OPTION  "FREE" 


2 


3 

(etc.  ) 


List  of  node  component  codes 

CODE  =  N0DE*10+I 

where  I  =  1,  2  or  3 

for  X,  Y,  or  Z  direction,  respectively 

Maximum  list  length  is  99  (use  additional  FREE 
records  to  get  more). 


NOTES; 

1.  FREE  causes  a  0  to  be  placed  in  the  constraint  rode  array  (NODFIX) 
for  the  indicated  component.  The  component  will  remain  free  until  reset 
by  a  FIX  record  or  a  predefined  LIMIT  record. 

CAUTION:  Do  not  FREE  any  components  of  a  slave  node. 
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7.1.4  IMPACT  -  Impacting  body 
Variable 

Word  Name _  Description 


1  OPTION 

2  IMPNOO 

3  IMPBOD 

4  VIB  (1) 

5  VIB  (2) 

6  VIB  (3) 

7  I  OPT 


"IMPA" 


Node  where  impact  occurs 


Body  number  (see  BODY  data  records) 


Components  of  body  velocity  at  impact 


2  f 

Weight  option 


0  -  floating  body,  no  weight  added,  only 
adds  inertias  (e.g.,  iceberg  or  snagged 
vessel) 

1  -  new  lumped  body  added  to  the  system 
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7.1.5  INITIAL  *  ll«si:rib«s  dynamic  initial  conditious 


Variable 


Word 

Name 

Description 

1 

OPTION 

"INIT" 

2 

VI  (1) 

V  \ 

X  1 

3 

VI  (2) 

Vy  \  Components  of  nodal 

4 

VI  (3) 

V  1 

z  f 

5 

NBEG 

Beginning  node  number 

6 

NEND 

Ending  node  nuinber 

7 

NINC 

Node  number  increment 

NOTES : 

1.  If  Word  (5)  is  zero  or  is  not  given,  the  velocity  components  are 
assumed  to  be  on  all  nodes  in  the  system.  In  this  case,  the  last.  INIT 
record  encountered  is  the  one  used. 

2.  Wlieu  Word  (5)  througlt  Word  (7)  are  not  zero,  the  velocity  components 
are  assigned  to  the  individual  nodes.  Repeat  the  INIT  records  as  many 
times  as  needed  to  define  all  nodal  velocities. 
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7.1.6  KEEP  -  Retains  the  data  from  the  pioeeeding  SAO 
Variable 

Word  _ Name  Description 


1  OPTION  “KEEP" 

MOTES; 

1.  This  flag  causes  the  data  initialization  to  he  bypassed,  thereby 
allowing  the  preceding  SAO  data  to  remain  in  effect.  Only  those  data 
to  be  changed  need  be  entered. 

2.  KEEP  must  be  the  first  flag  encounteretl  in  the  SAO  data  set  it  it  is 
used. 
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7.1.7 

LOAD  -  Specifies 

load  coudi lions 

Wo  rd 

Var  J  ,ible 

l)es*'ript  ion 

1 

OlTHiN 

"LOAD’' 

2 

l.SKT 

Load  set  number  (3  max,  default  -  1) 

3 

FF  (1,SKT,  I) 

F 

X 

FP  (LSKT,  1+1) 

F 

y 

i  Pojnt  load  components  (F  or  Fl.) 

1 

,5 

FP  (LSET,  1+2) 

1 

6 

NIIFG 

Peg in  node  number 

7 

MEND 

End 

node  number 

ft 

NINO 

Nod«'  nnmb<M'  increment  Idef.iulL  :=  1) 

NOTO: 

1.  I.oads  arc  assigned  as  inembors  of  up  lo  lliroo  load  sots  according  lo 
LSKT.  The  IjVAKY  record  file  is  llieii  used  lo  ideutily  the  characteristics 
of  carl)  load  set. 

2.  I.,oad  components  arc  placed  in  the  point  load  array  Ki'  l>y  load  set 
number  and  node  numbers.  Word  ((>)  must  .ilways  he  given.  Tlio  saim'  io.ad.s 
can  be  applie»i  to  a  sequence  of  nodes  by  giving  appropriate  v.ilues  Jin 
Word  (7)  and  Word  (ft). 

d.  If  the  node  is  used  for  the  angular  position  of  a  r>gid  body,  the 
load  is  assumed  to  be  a  inomenl  about  the  specified  axis  (units  Fl.) 
Otlierv'ise  the  loads  are  point  forces  (units  F) . 
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7.1.8 

LVAHY  “  Specifies 

load  vari.'ition  cod«‘s 

Word 

Variable 

Name 

l)<'scr  ipt  ion 

1 

OFT  ION 

"l.VAK" 

2 

ILF  (1) 

Variation  code  for 

load 

set 

1 

;i 

ILF  (2) 

Va  r i a  t i on  code  f o  r 

loa  d 

set 

2 

4 

TLF  (3) 

Variation  code  for 

load 

set 

3 

5 

.JDLl) 

Or.avity  load  vari.at 

i  on 

code 

f  o  r 

siiliaiu»lyK<?s  (see  Noi<>  2) 

NOTKS: 

1.  Variat  ion  (loile : 

For  nKAl)  and  1.1  Vf: 

1  '  increment  the  load  set  factor  from  0  to  1.0  (apply  load) 

0  •  hold  the  load  sot  factor  at  1.0  (.steady  load) 

-1  -  ifioiemenl  the  lo.'id  .sot  factor  from  1.0  to  0  (remove  load) 

iMai-DYN 

0  "  hold  the  sot  factor  at  1.0  (steady  load) 

1>{)  *  use  the  r‘'time  variation  function  from  the  TFIIN  record 
sot  to  get  the  load  set  factor 

2.  The  only  loads  tised  in  a  DKAl)  SAO  are  the  internally  calculated 
gravity  loads  and  point  loads  defined  in  l.OAl)  record  data  sots.  The 
gravity  loads  are  assnnie<i  to  have  a  +1  variation  code  in  the  DKAl)  SAO. 
The  MVK  SAO  assnim's  a  variation  code  of  0  for  gravity  load.s.  When  a 
LIVE  SAO  has  not  l)een  preceded  by  a  DKAl)  SAO  it  may  be  advantageous  to 
inerenK'ut  tin?  gravity  load.s  along  with  the  LIVE  load.s.  Setting  .)1)1,D  to 

will  accomplish  this. 
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7.1.9 

MOVE  -  Specifies 

node  point  motion 

Word 

Variable 

Name 

Description 

1 

OPTION 

"MOVE" 

2 

INOD 

Node  number  to  be  moved 

3 

MC  (1) 

X  motion  code 

4 

MVC  (1) 

X  motion  variation  code 

5 

AMP  (1) 

X  motion  amplitude 

6 

MC  (2) 

Y  motion  code 

7 

MVC  (2) 

Y  motion  variation  code 

8 

AMP  (2) 

Y  motion  amplitude 

9 

MC  (3) 

Z  motion  code 

10 

MVC  (3) 

Z  motion  variation  code 

11 

AMP  (3) 

Z  motion  amplitude 

NOTES; 

1 . 

'  on  codes: 

0  -  no  motion  specified  (see  Note  4) 

1  -  displacement 

2  -  velocity  (DYN  and  TSSS  only) 

3  -  acceleration  (DYN  and  TSSS  only) 

2.  Motion  variation  codes: 

For  DEAD  and  LIVE 

1  -  increment  motion  factor  from  0  to  1.0  (apply  deflection) 
0  -  hold  the  motion  factor  at  1.0  (hold  deflected  position) 

The  motion  factor  is  the  multiplier  of  the  motion  amplitude  for  the 
increment  level  (load  factor)  at  the  current  step. 

For  DYN  and  TSSS 

0  -  hold  th>  motion  factor  at  1.0 
I>0  -  use  the  time  variation  function  from  the  TFUN  record 

set  to  get  the  motion  factor 
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3.  The  motion  amplitude  will  be  an  angflar  displacnneut .  velocity,  or 
acceleration  on  nodes  used  to  define  angular  motion  of  rigid  bodies. 

U.  DYN  .and  TSSS  require  all  three  motion  components  at  a  node  to  be 
defined.  Individual  components  can  be  defined  in  DEAD  ana  LIVE.  A  zero 
motion  variation  code  in  DYN  and  TSSS  means  a  fixed  component  that 
reinaif\s  at  its  initial  position. 

5.  The  maximum  niimber  of  components  (X,  Y,  Z)  of  imposed  displacements 
in  DEAD,  LIVE,  and  MODE  is  30.  The  maximum  number  of  moved  nodes  in  DYN 
and  TSSS  is  5. 

6.  The  constraint  code  array  (NODFIX)  entry  i.s  superseded  by  a  MOVE 
iiistri'ction .  The  constraint  code  is  restored  at  the  completion  of  the 
SAO. 
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7.1.10  OUTPUT  -  Specifies  output  data 


Variable 
Word  Name _ 

3  OPTION 

2  NIPR 

3  DTOUT 


4  IBG 

5  TFXOUT 

NOTES: 


Description 

"OUTP” 

Number  of  solution  steps  between  output  prints 
Output  time  interval  (T)  (DYN  and  TSSS  only) 
or 

Output  angle  interval  (degrees)  (LIVE  with 
heading  change) 

Debug  output  parameters  (see  Appendix  I) 

Output  current  list  of  node  component  fixity 
codes  if  this  is  nonzero. 


1.  NIPR  is  used  to  calculate  the  output  interval  on  DEAD  and  LIVE  by 
multiplying  the  first  step  size  by  NIPR.  The  output  interval  remains 
the  same  even  though  the  step  size  is  changed  unless  NIPR  =  1.  When 
NIPR  =  1,  tlie  output  is  at  every  step,  regardless  of  step  size. 

2.  DYN  and  TSSS  can  use  NIPR  or  DTOUT  to  determine  the  output  interval. 

If  both  are  given,  DTOUT  has  precedence. 

3.  Output  is  given  at  the  beginning  and  end  of  each  subanalysis. 

Output  is  also  given  when  the  load  or  time  parameter  exceeds  the  value 
of  the  parameter  at  the  last  output  plus  the  output  interval.  Extra 
output  records  can  be  produced  when  a  restart  save  request  does  not 
correspond  with  an  output  interval  (see  SAVE  record). 

4.  Output  for  LIVE  SAO  with  heading  changes  will  occur  at  angle  intervals 
.selected  by  DTOUT.  If  DTOUT  is  not  given,  the  results  will  be  output  at 
every  angle  calculated.  NIPR  has  no  effect  on  heading  change  output. 
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7.1.11 

Word 

1 

2 

3 

4 

5 

6 

7 

8 
9 

NOTES; 


PAYOUT  -  Payout/reel-in  data 
Variable 

Name _  Description 

OPTION  "PAYO” 

JOP  (1)  Initial  node  where  payout/reel-in  occurs 

JPELT  (I)  Initial  element  number  for  payout 

PAYV  (I)  Payout  rate  (LT”^) 

+  for  payout 
for  reel-in 

AMAXL  (I)  Mitosis  length  (see  Notes)  (L) 

NGROW  (I)  Number  of  elements  available  for  growth 

NSHRINK  (I)  Number  of  elements  available  for  reel-in 

NELPOI  (I)  Element  number  increment  (default  =1) 

NPOVRY  (I)  Number  of  time  variation  function 

0  -  constant  rate 

n>0  -  use  n^  function  defined  in  TFUN 
record 


1.  Payout/reel-in  is  approximated  by  incrementally  changing  the  un¬ 
strained  element  length.  When  the  unstrained  length  exceeds  AMAXL  (1) 
plus  a  reference  length  during  payout  and  NGROW  (I)  >  0,  then  a  mitosis 
operation  will  be  performed.  The  payout  element  will  be  divided  into 
two  elements.  The  original  payout  node  will  be  assigned  to  the  division 
point,  .'nd  the  new  node  introduced  in  the  new  element  will  become  the 
payout  node.  The  new  element  is  obtained  by  adding  NELPOI  (1)  and  JPELT  (1). 
During  reel-in  the  reverse  process  will  occur  when  the  unstrained  length 
is  less  than  AI'AXL  (I)  and  NSHRINK  (I)  >  0.  In  this  case  the  new  element 
is  identified  by  subtracting  NELPOI  (I)  from  JPELT  (1).  AMAXL  (I)  =  0.0 
causes  the  mitosis  check  to  be  ignored.  The  reference  length  for  payout 
mitosis  is  one  of  the  following: 

A  -  The  initial  upstretxhed  length  of  the  payout  element  (JPELT  (I)) 
if  the  next  element  to  be  activated  has  a  different  material 
number  or  if  this  is  the  first  mitosis  for  that  payout  end. 

B  -  AMAXL  (I)  for  all  other  situations. 
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2.  The  elements  available  for  payout  are  inactive  until  a  mitosis  acti¬ 
vates  them  by  assigning  a  length  to  them.  These  elements  must  be  included 
in  the  total  number  of  elements,  and  their  nodes,  material,  etc.,  must  be 
defined  in  the  usual  manner.  The  nodes  for  these  inactive  elements  must 
also  be  defined  in  the  usual  manner  with  their  coordinates  defined  to  be 
the  same  as  the  initial  position  of  the  payout  node  (causes  the  element 
length  to  be  zero).  These  nodes  must  be  given  constraint  codes  for  all 
degrees-of-f reedom.  The  appropriate  codes  are: 

1  -  if  the  component  is  to  remain  fixed  and  subject  to  limit  checks 

(if  any)  after  mitosis 

2  -  if  the  component  is  to  be  free  after  mitosis  but  possibly 

subject  to  subsequent  limit  checks 

3  -  if  the  component  is  to  remain  fixed  and  not  subject  to  limit 

checks  after  mitosis 

3.  The  elements  available  for  reel-in  must  be  active  portions  of  the 
structure  and  sequentially  connected  to  the  reel-in  point  (no  branches). 
The  nodal  constraint  codes  on  the  reel-in  node  will  be  assigned  to  the 
new  reel-in  node,  and  the  deactivated  nodes  will  be  appropriately  con¬ 
strained  automatically. 

A.  Payout/reel-in  can  occur  only  at  fixed  nodes  or  nodes  defined  in 
MOVE  data  sets.  In  either  case,  the  displacements  of  all  three  components 
must  be  defined. 
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SAVE  -  Defines  restart  tile  save  intervals 


Variable 

Word  Name 

1  OPTION 

2  IRST 


3  DTRSRT 

A  NIXPRN 

NOTES: 

1.  Restart  files  are: 


Description 

"SAVE" 

Restart  file  flag 

0  -  do  not  save  records 
-n  -  rewind  file  and  save  every  n  output 
record 

nil  -  e^^end  present  file  and  save  every 
n  output  record 

Restart  save  time  interval  (T)  (DYN  and  TSSS 
SAO  only) 

Output  suppress  flag  (see  Note  3) 


01  -  DEAD 

02  -  LIVE  (TSSS) 

03  -  DYN 

2.  Vfhen  DTRSRT  is  given  for  DYN  or  TSSS,  IRST  is  used  only  to  signal 
file  rewinds.  The  restart  data  is  then  written  when  the  current  time  is 
greater  than  or  equal  to  the  time  of  the  last  save  plus  the  save  time 
interval . 


3.  If  the  step/time  on  the  SAVE  record  does  not  coincide  with  the  OUTPUT 
record  request,  then  an  extra  output  record  will  be  produced  unless 
NIXPRN  f  0. 
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7.1.13 

SBUOY  -  Specifies 

vertical  motion  for  a  body  on  the  surface 

Word 

Variable 

Name 

Description 

1 

OPTION 

"SBUO" 

2 

NODSUR 

Node  where  body  is  located 

3 

NMOTN  (I) 

Motion  code 

1  *  displacement 

2  -  velocity 

3  -  acceleration 

A 

IBS  (I) 

Time  function  number  (see  TFUN  record, 
Section  6.18) 

5 

SBAMP  (I) 

Surface  motion  amplitude 

NOTES: 


1.  The  motion  indicated  by  these  parameters  will  be  imposed  on  the 
vertical  component  whenever  the  body  is  held  to  a  limit.  This  means  the 
body  node  must  appear  on  the  LLOC  record.  When  the  limit  restraint  is 
exceeded,  the  body  is  released  from  the  motion  and  the  limit  condition. 
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7.1.14 

SOLUTION  - 

Specifies  solution  option  characteristics 

Word 

Variable 

Name 

Descript  Lon 

1 

OPTION 

"SOLU" 

2 

SOPTN 

Solution  method: 

MNR  (default  for  DEAD,  LIVE,  TSSS,  FREQ) 

SLl 

RFB 

VRsJ  ,  , 

}  static  only 

DIM  (default  for  DYN)  dynamic  only 

3 

DMU 

Numerical  damping  factor  (see  Note  1) 

4 

RERR 

Residual  norm  error  bound  (default  =  0.001) 
(see  Note  2) 

5 

DERR 

Displacement  and  pseudo-velocity  norm  error 
bound  (default  =  0.001)  (see  Note  2) 

6 

DALPHA 

Proportional  damping  multiplier  of  mass  matrix 

(DYN) 


or 

Alpha  integration  parameter  (VR  solutions  - 
detauH  =  1.0) 

or 

Ship  angle  damping  for  MNR 

7  DBETA  Proportional  damping  multipl ier  of  stiffness 

matrix  (DYN) 

or 

Initial  pseudo-time  step  (VR  solution.s  - 
default  =  1 .0) 

or 

Buoy  angle  damping  for  MNR 

8  SRCHFC  ID  search  factor  for  MNR  method 

=0.0  -  no  ID  search  of  alternating  estimates 
<0.0  -  no  alternating  estimate  checks 
>0,0  -  the  ID  search  initial  guess  factor 
(see  Note  4) 
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Variable 

Word 

Name 

Description 

9 

PARMT 

MNR  method  extrapolation  parameter  (default  = 
0.5)  (see  Note  5) 

10 

CEPS 

Strum  update  parameter  (default  =  0.0001) 

0.0  -  no  drag  amplifications  for  this  SAO 
<0.0  *  calculate  drag  amplifications  only 
once  per  SAO 

>0.0  -  calculate  drag  amplification  at  the 
beginning  of  this  SAO,  and  each  time 
the  relative  velocity  norm  changes 
more  than  CEPS  (sec  Note  6). 


11 

NUP 

Update  option  flag  (let  default) 

12 

JMPDT 

Step-size  control  number  (default  =  2) 

13 

LMITER 

Iteration  limit  (default  =  200) 

14 

KONVRT 

Number  of  trials  before  divergence  abort 
(default  =  3) 

15 

NRUP 

Newton-Raphson  update  interval.  Used  on  MNR 
method  to  signal  how  often  to  recalculate 

the  tangent  stiffness  matrix,  K.J..  This  is 
a  packed  word; 

•  Units  and  tens  digits  give  number  of  steps 
before  new  K.p. 

•  Hundreds  and  above  digits  give  the  number 
of  alternating  sign  trials  on  a  given 
step  before  new  K.j,  (default  =  205,  i.e., 
two  alternating  tries  and  five  steps) 

NOTES : 


1.  The  numerical  damping  factor,  DHll,  is  used  to  avoid  problems  with 
singular  or  ill-conditioned  stiffness  matrices.  It  has  no  influence-  on 
the  DIM  method.  It  should  he  used  with  caution  with  the  SLl  and  KFB 
methods  since  it  alters  the  stiffness  matrix,  hence  the  equilibrium 


equations , 

in  these  two  cases 

Possible 

DMU  i?  1.0 

very  heavy 

1.0 

^  DMU  ^0.1 

heavy 

0.1 

^  DMU  ^  0.01 

moderate 

0.01 

^  DMU  i  0.001 

light 

0.001 

g  DMU 

very  light 
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If  DHU  <  10  and  a  singularity  is  oiicountorpd ,  thon  DMu  is  sol  to  0.001, 
and  the  step  is  tried  again.  On  a  repeated  singularity,  DMU  is  increased 
to  0.1,  and  the  step  is  tried  once  more.  If  the  .singularity  pcrsi.sts, 
the  calculation  is  aborted. 

2.  The  error  bounds  are  used  to  test  for  convergence  of  tlie  iterations. 
RERR  is  used  only  in  the  MNR  and  VRR  methods.  DERR  is  used  in  MNR,  VRK, 
VRS  and  DIM  methods, 

The  convergence  criteria  are: 

RNORM  ^  /RERR/ 

DNORM  ^  /DERR/ 

where  RNORM  and  DNORM  represent  norms  of  the  nodal  force  residual  and 
displacement  increments,  respectively.  Both  of  these  inequalities  must 
be  satisfied  for  convergence  of  the  MNR  iterations,  while  only  the  dis¬ 
placements  are  checked  in  the  DIM  method.  The  VRR  solution  checks  the 
residual  (RERR)  as  well  as  displacements  and  velocities  (DERR).  The 
residual  is  not  used  in  the  VRS  solution. 

Some  flexibility  in  the  form  of  these  norms  is  available.  If  the 
value  of  RERR  is  greater  than  zero,  then: 


RNORM  = 


where:  ••  the  i  compoutMJt  of  the  nodal  residual 

N  =  the  total  number  of  nodal  degrees-of-freedom 
T 

max  =  the  maximum  element  tension 
If  the  value  of  RERR  is  negative,  then: 


RNORM  = 


N 

I  R 
i=l  ^ 

”  2 
I  R 

i=l  ^ 


where:  M  =  the  total  number  of  nodal  components  including  the  fixed 

nodes;  i.e.,  the  reactions  arc  included  in  the  denominator 
(MNR  method  only) . 

If  the  value  of  DERR  is  greater  than  zero  on  static  analyses,  then: 


DNORM 


/  ^  2 
=  f  I  A(Aq,)VNl 


where: 


the  i*"*'  component  of  the  change  in  the  displacement 
increment . 


If  DERR  is  negative  on  static  analyses  (except  VR  methods),  then: 


DNORM 


In  the  DYN  snbanalysis,  the  norm  selection  is  the  opposite  of  that 
selection  used  in  static  analyses.  See  the  discussion  of  norms  in 
Reference  8. 


3.  The  terms  DALPHA  and  DBETA  can  be  used  to  specify  internal  damping 
proportional  to  the  mass  and/or  stiffness  matrix.  The  assumed  form  of 
the  d.imping  matrix  is: 

(C]  =  a  (M)  +  p  IK] 

This  procedure  can  be  used  in  DYN  subanalysis. 

A.  The  ID  sc'arch  factor,  SRClfFC,  can  bo  used  with  the  MNR  method  in  an 
attempt  to  enhance  a  poor  initial  configuration.  See  the  discussion  of 
the  MNR  method  in  Reference  8. 


5.  PARMT  is  an  extrapolation  parameter  used  in  the  incremental  MNR  method. 
The  starting  estimate  for  the*  displacement  on  all  but.  tin*  first  :>tcp  is 
the  accumulated  displaccaiients  freun  the  previous  step  plus  PARMT  Limes 
the  change  in  displacement  calculated  from  the  preceding  step. 


6.  When  strum  strings  are  defined  in  the  STRUM  recor'd  (Section  6.16) 
and  CEl'S  is  not  zero  for  LIVE,  DYN  or  TSSS  SAO's,  drag  amplification 
factors  will  be  coniput.ed  for  the  strings.  This  computation  will  be  done 
at  the  beginning  of  the  subanalysis  and  whenever  the  relative  velocity 
between  the  fluid  and  cable  elements  changes  significantJy .  A  significant 
change  is  indicated  by: 


VNORM.  -  VNORMo 
_ r _ R 

VNORM„ 


^  CEPS 


where  VNORM 


IVF(I)  -  UD(I))' 


il/2 


R  -  time  when  drag  amplifications  were  last  calculated 
i  =  present  time 

VF  =  fluid  velocity  components  at  each  node 
1)1)  =  nodal  velocity  components  (DYN  only) 
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7.  1  Ab 

STEP  -  Solution 

step-size  data  (static  solntions) 

Wo  rd 

V.i  riable 

Name 

Descript  ion 

1 

OPTION 

"STEP" 

2 

JSTEPR 

Number  of  steps  for  load  incrementing 
(default  =  1) 

3 

NSTUP 

Start  up  option 

n>l  -  divide  first  step  into  n  substeps 
in  an  arithmetic  progression 

k 

IIEDINC 

Surface  load  heading  increment  (degrees) 

3 

HKDEND 

Surface  load  total  heading  change  (degree) 

NOTES: 

1.  JSTEl'R  indioal«'s  the  number  of  loacl/niovoment  steps  requested. 

During  this  phase  of  the  SAG,  ttw'  load  and  movement  magnitudes  are  varied. 
When  HEUINC  is  not  zero,  a  sequence  of  static  configurations  will  be 
generated  following  the  convergence  or  completion  of  the  magnitude 
incrementing.  This  additional  loading  sequence  moves  the  surface  wind 
and  current  headings  through  the  excursion  defined  by  HKDENl)  and  the 
heading  variation  codes. 

2.  Se<;  Note  9  of  Section  6,1A  (NOl)K  record)  for  global  heading  conventions. 

3.  See  Section  7.1.16  (SURF  record)  and  Section  7.1.18  (WIN!)  record) 
for  initial  lieadings  and  variation  codes. 
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7.1.16  SDRFACK 


Surf:ii:<.'  cinreut  daUi 


Word 

Variabl e 

Name 

Description 

1 

OPTION 

".SURF" 

2 

CURNT 

Surface  current  velocity  (I.T  ^) 

3 

GAD 

Initial  current  heading  (degrees) 

4 

KODEC  (i) 

Variation  code  for  fixed  heading  (load  varies) 

5 

KODEC  (2) 

Variation  code  for  variable  heading  (load 
constant) 

NOTES: 


1.  See  Note  9  for  Section  6.14  (NODE  record)  for  globnl  heading  convention. 

2.  See  Section  7.1.15  (STEP  record). 

3.  Variation  codes: 

0  -  hold  constant  at  value  given 
-I  -  increment  down  (not  for  heading  variation) 

1  -  increment  up 

4.  Shipr.  and  buoys  located  at  the  surface  limit  the  use  this  of  those 
data  and  the  auxiliary  loading  data  to  calculate  static  loads.  Ships 
will  obtain  loading  coefficients  from  a  data  file  or  built-in  functions 
(see  SHIP  record,  Section  6.15,  and  Appendixes  B,  F,  and  G) .  Buoys  will 
use  the  BWND  and  BSCD  parameters  (see  BODY  record,  Section  6.5).  If  the 
buoy  is  pxilled  away  from  the  limit,  its  loads  will  he  calculated  from 
the  GURR  record  (Section  7.1.1)  defined  flow  field  and  the  body  drag 
functions  (built-in  or  user  defined).  See  Word  (2)  on  BODY  record. 

Section  6.5. 
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/ .  1 . )  7 

TIME  -  Time 

step  data 

Word 

Variabl  e 
Name 

Descripl ion 

1 

OPTION 

"TTME" 

2 

DT 

Time  step 

3 

TMAX 

Final  time 

A 

T 

Initial  time  (defaults  to  0) 

5 

DTU 

Update  time  (defaults  to  DT);  always  let 
default  for  DIM  solutions 

6 

ALPNEW 

a  integration  parameter  (default  =  0) 

7 

BETNEW 

P  integration  parameter  (default  =  1/12) 

GAMNEW 

Y  integration  parameter  (default  =  1/2) 

NOTES: 

1.  If  DT  i  0.0,  the  time  step  will  be  internally  estimated.  It  DT  =  -A, 
then  the  estimated  value  will  be  multiplied  by  A.  If  DTU  =  “B,  then  DTO 
is  set  to  B  times  the  DT  estimate. 

2.  The  integration  parameters  are  those  of  the  generalized  Newmark 
difference  equations  (Ref  A). 

3.  A  nonzero  initial  time  can  be  used  tor  TSSS  or  DVN  SAG  to  udjust  to 
time  in  the  TFUH  records. 
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7.1.18 

WIND  -  Surface 

wind  data 

Word 

Variable 

Name 

Description 

1 

OPTION 

"WIND" 

O 

4- 

WIND 

Wind  velocity  (LT  ^) 

3 

WAD 

Initial  wind  heading  (degrees) 

4 

KODEW  (I) 

Variation  code  for  iixed  heading  (load  varies) 

5 

KODEW  (2) 

Variation  code  for  variable  lieading  (load 

constant) 


NOTES: 

1.  See  Note  9  of  Section  6.14  (NODE  record)  for  global  heading  convention. 

2.  See  Section  7.1.15  (STEP  record) 

3.  Variation  codes; 

0  -  hold  constant  at  value  given 
~1  -  increirent  down  (not  for  heading  variation) 

1  -  increment  up 


7,J 


7.2  MODE  -  Mode  shape  calculation 


tfSOL  selects  solution  options.  An  MSOL  recor  is  no  n 
defaults  are  acceptable. 


7.2.1  MSOl  -  Specifies  solution  format  parameters 
Variable 

Word  _ Name  Description 

1  OPTION  "MSOL" 

2  MODEIl  Mode  shape  order  flag 

0  -  list  mode  shapes  in  order  of  increasing 
f  •'equency 

1  -  list  in  order  calculated 

3  M0DEI2  Mode  shape  output  flag 

0  -  print  all  mode  shapes 
n  -  print  n  mode  shapes 
-n  -  print  n  mode  shapes  and  write  them  on 
logical  unit  20 


7!i 


7.2.2  FIX  -  Applies  fixed  conditions  to  nodes 

This  record  is  Identical  to  FIX  in  Section  7.1.2. 
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7.3  FREQ  -  Response  spectrum  calculation 


The  frequency  domain  subanalysis  data  set  is  composed  of  three 
levels  of  data  following  the  FREQ  flag  record.  The  first  level  describes 
the  setup  of  the  implied  LIVE  solution  for  drift  force  updates.  The 
second  level  specifies  the  setup  of  the  frequency  domain  solution.  The 
third  level  gives  the  data  to  be  processed  for  each  wave  heading  considered. 

The  characteristics  of  the  implied  LIVE  step  to  provide  drift  force 
updates  can  be  specified  by  the  appropriate  records: 


OUTP 

SAVE 

SOLU 

STEP 


See  Section  7.1 


These  must  appear  between  the  FREQ  keyword  and  any  of  the  FREQ  setup  cards. 

The  data  records  for  the  FREQ  solution  setup  are: 

FSOL  -  selects  frequency  domain  options 

SPECTRUM  -  specifies  spectral  characteristics  (required) 

EXTERNAL  -  ship  motion  file  conversion  factors 

RESULTS  -  optional  output  requests 

Each  of  these  is  optional  except  for  SPEC,  and  can  be  given  in  any 
sequence  in  advance  of  the  first  HEAD  record. 

The  data  records  for  the  wave  headings  are  grouped  by  HEAD  flag 
records.  The  HEAD  record  sets  contain  one  or  more  RAND  or  RECU  record 
sets . 
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7.3.1  FSOL  -  Specifies  frequency  domain  options 


Variable 
Word  Name 

1  OPTION 

2  IDRITR 


3  IFRQUP 

U  ICMCHF 

5  ICONMS 

6  IROLIT 

7  ALPHA 

8  BETA 


Description 

"FSOL” 

Drift  force  iteration  flag  (see  Note  1) 

0  -  no  drift  force  iterations 
1  -  perform  drift  forces  iterations  to 
convergence 

-1  -  do  a  single  drift  force  update  with  no 
repeat  of  frequency  solution 
±2  -  same  as  il  except  use  only  surge  force 
(for  single-point  moorings) 

Iteration  option  flag  for  drift  force  updates 

0  “  start  from  original  static  reference 
1  -  start  from  last  configuration  found 

Component  check  file  flag 

0  -  save  response  file  OA  for  use  by  CHEK 
record  (see  Note  2) 

Mass  matrix  format 

0  -  lumped  mass 
1  -  consistent  mass 

Roll  damping  iteration  flag 

0  -  do  not  iterate 

n  -  iterate  no  more  tlian  n  times  to  estimati 
ship's  roll  damping 

Proportional  damping  multiplier  of  mass  matrix 

Proportional  damping  multiplier  of  stiffntiss 
matrix 


NOTES: 

1.  If  the  static  reference  st.ite  is  to  be  modified  for  the  effects  of 
wave  drift  forces,  then  the  accumulated  drift  forces  for  a  pass  througli 
the  wave  spectrum  must  he  statically  applied.  A  single  up<iate  simply 
assumes  the  dynamic  solution  is  unchanged  by  the  new  static  reference 
state.  The  iterative  solution  repeats  t  lie  frequency  pass  until  negligii>lc 
changes  in  position  due  to  drift  forces  are  ca I < ul ated .  fhis  convergence 
is  checked  using  the  following  measure; 
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1/2 


DN 


N 

I 

i=l 


/DN'  -  DN/  ^ 
DN 


DERR 


where;  X  =  nodal  point  coordinates 

N  =  total  number  of  coordinates  (3  times  the  number  of  active 
nodes) 

DERR  =  displacement  error  tolerance  (see  SOLD  record,  Section  7.1.1^) 
The  prime  indicates  the  next  value  of  DN. 

2.  Dynamic  data  for  the  CHEK  record  must  be  passed  on  file  04  or  the 
CHEK  results  will  include  only  the  static  reference  state. 

3.  Roil  iterations  at  each  frequency  seek  to  converge  on  a  roll  angle 
estimate  by  adjusting  the  roll  damping.  Convergence  is  assumed  if  a 
change  of  less  than  1  degree  is  found  on  two  successive  cycles  or  when 
the  number  of  cycles  equal  IROLIT. 

A.  Internal  damping  proportional  to  the  mass  and/or  stiffness  matrix  is 
used  when  ALPHA  and/or  BETA  are  non-zero.  The  form  is: 

(Cl  =  0(  (Ml  +  p  IK] 


81 


7.3.2  SPKCTRUM  -  Specifies  spectral  characteristics 


Word 

Variable 

Name 

Description 

1 

OPTION 

"SPEC" 

2 

SPECA 

A  coefficient  (L^  T  ^) 

3 

SPECB 

B  coefficient  (T  ^) 

♦ 

DOMG 

Am,  frequency  increment  (T  S 

5 

OMGMN 

u)  .  ,  lower  bound  on  frequency 
min’  ^  ^ 

(T‘^) 

6 

OMGMX 

uppi^r  bound  on  frequency 

(t‘^) 

7 

AMPMN 

Cut-off  amplitude  for  waves  (L) 
(default  =  0.0001,  see  Note  2) 

NOTES: 

1 .  The 

wave  spectrum 

is  assumed  to  have  the  form: 

S(U)) 


A/u)^ 


where  S(ui)  is  based  on  twice  tl>e  square  of  the  wave  height.  Any  spectrum 
having  this  general  form  can  be  input.  Values  for  common  spectra  arc: 


SPECTRUM 


A 


Pierson-Moskowi  t7.  135.0 


where 


S 


Bretschneider  A200  II  ^/T  ^ 

s  s 

I.S.S.C.  2760  II  2/T  ^ 

s 

V|^  =  wind  speed  (knots) 

11^  -  significant  wave  height  (ft) 

=  significant  wave  period  (sec) 
u)  “  circular  frequency  ( radians/scc) 
(u))  =  spectral  ordinate  (ft^-sec) 


B 

4  A 

9.7  X  10  /V|^ 

1050/T^^ 

690/T^^ 


2.  Regular  wave  re.sponse.s  are  calculated  for  waves  with  frequencies  be¬ 
tween  U)  .  and  ui  ,  which  have  wave  amplitudes  (based  on  the  specLiimt) 
greater'  ^ilian  Increment  ing  will  not  procee<l  beyond  u)  regard  .-ss 

of  the  wave  amplitude. 
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7.3.3 

RXTEjW^  -  Specifies 

external  ship  motion  file  conversion  factors 

Word 

Variable 

Name 

Description 

1 

OPTION 

"KXTE" 

2 

FRCFAC 

Force  conversion  factor 

3 

FACLEN 

Length  conversion  factor 

4 

TIMFAC 

Time  conversion  factor 

5 

NFILEF 

File  format  code  (default  =  0) 

J  -  NSRDC  ship  motion  file 
0  -  NCKL  ship  motion  I'ilo 


NOTES : 

1.  The  conversion  factors  are  used  as  multipliers  of  the  data  on  the 
external  coefficient  file  (ship  motion  file).  Eacli  of  them  has  a  default 
value  of  1.0. 

2.  The  file  formats  are  described  in  .Appendi.x 


7.3.4  RESULTS  •  Speoit'i«?s  optional  output  requests 
Variable 

\^rd  _ Name  Description 

1  OPTION  "RESU" 

2  lUNRKS  Unrestrained  buoy  and  ship  motion  outputs 

0  -  no 
1  -  yes 

3  IBGF  Debug  output  flag  (see  Appendix  I) 


NOTES ; 

1.  The  IBGF  flag  produces  extra  output  for  the  frequency  response 
solution  only.  The  IBG  flag  (Section  7.1.10,  OUTPUT  record)  will 
produce  extra  output  only  for  the  implied  I.IVK  soliition  used  for  diift 
force  updates. 
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7.3.5  HDADING  ■  Wavt;  iicadiug  record 
Variable 

DescHi>tion 

1  OPTION  "HEAD" 

2  GHED  Wave  heading  in  global  system  (degrees) 

NOTES: 

1.  HEAD  data  sets  consist  of  the  record  followed  by  a  series  of  records 
of  RAND  and/or  HEGU  re.sponse  requests.  Each  HEAD  data  set  can  be  terminated 
by  a  DONE  record,  another  HEAD  record,  or  a  valid  SAO  flag. 


7.3.6 

RANDOM  -  Random 

Variable 

response  data  requests 

Word 

Name 

Description 

1 

OPTION 

"RAND" 

This  flag  record  is 
following  form: 

Variable 

followed  by  a  string  of  records  having 

the 

Word 

Name 

Description 

1 

VTYPE 

"NODE",  "SHIP",  "TENS",  "DONE"  (see  Notes) 

2 

NUMC 

Node  or  element  number 

3 

NDRT 

If  VTYPE  =  "NODE"  gives  the  global 
direction 

1  =  X 

2  =  Y 

3  =  Z 

component 

A 

SFSTAT 

Static  load  factor  ("TENS"  only) 

5 

SF3 

(1/3)  load  factor  ("TENS"  only) 

6 

SFIO 

(1/10)  load  factor  ("TENS"  only) 

7 

NOTES: 

STIOO 

(1/100)  load  factor  ("TENS"  only) 

1.  The  spectral  response  of  the  ship,  element  tensions,  and  any  of  the 
nodal  displacement  components  can  be  calculated.  VTYPE  signals  the  one 
that  is  desired.  Any  number  of  these  cards  can  be  provided.  The  input 
is  terminated  by  a  card  with  VTYPE  =  "DONE".  Termination  can  also  be 
accomplished  with  a  REGU,  HEAD,  or  any  valid  SAO  flag.  The  response 
data  calculated  represent  the  average  of  the  1/3,  1/10,  1/100  highest 
responses . 

2.  VTYPE  =  "SHIP"  requires  no  other  entries  on  the  card. 

3.  VTYPE  =  "TENS"  requires  the  specification  of  the  element  number  in 
NUMC.  If  values  are  given  for  SFSTAT  through  SFIOO  and  if  the  element 
material  has  the  ultimate  tension  specified,  then  estimates  of  the 
factored  maximum  loads  will  be  printed. 
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4.  VTYPE  =  "NODE"  requires  the  specification  of  the  node  number  in  NUMC 
and  the  component  direction  in  NDKT. 

5.  A  TENS  record  must  be  included  for  each  element  that  will  be  involved 
in  a  subsequent  coinponeat  adequacy  check.  For  example,  if  a  check  of  an 
anchor  capacity  is  to  be  made  including  dynamic  effects,  every  element 
that  connects  to  the  anchor  (or  fixed  node  where  the  anchor  is)  must  be 
called  out  on  a  TENS  record. 
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7.3.7 

REGUUR  -  Regular 

response  data  requests 

Word 

Variable 

Name 

Description 

1 

OPTION 

"RECU" 

This  flag  record  is 
following  form: 

followed  by  a  string  of  records  having  the 

Word 

Variable 

Name 

Description 

1 

NODE 

Node  for  which  response  is  requested. 

2 

LOCAL 

0  -  produces  output  in  global  system 

1  '  produces  output  in  ship's  local  system 

3 

WRFQ 

Circular  fro4uency  of  response  (T 

4 

WAMP 

Wave  amplitude  for  response  (L)  (default  =  1.0) 

5 

lOCODE 

0  -  output  is  displacements 

1  -  output  is  position 

6 

.NllM 

Number  of  divisioa.s  per  cycle  (default  -  30) 

NOTE^: 

1.  These  records  will  generate  the  displacement  or  position  versus  time 

for  all  three  nodal  components  through  one  cycle  of  inotinn. 

2.  The  local  coordinate  system  will  not  be  used  if  position  output  is 

requested , 

3.  There  is  no  limit  to  the  number  of  requests  that  Cf.n  be  made. 

4.  A  DOlfE,  RAND,  HEAD  or  any  valid  SAG  flag  record  will  terminate  this 
option. 
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7.3.8  DONF.  -  Optional  torminator  flag 
Variable 

Done  ___Name  Description 

1  OPTION  "DONE" 

NOTES: 

1.  Can  be  used  to  terminate  REGU  and  RAND  record  strings  and  HEAD  data 
sets.  Optional  terminations  are  taken  from  the  flag  records  if  a  DONE 
record  is  not  given.  Acceptable  terminators  are  RAND,  REGU,  HEAD,  or  any 
of  the  SAG  record  flags. 
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7.4.1  CONF  -  Wave  heading  configuration 

Use  only  if  previous  SAO  is  FREQ  with  more  than  one  wave  heading. 


Word 

Variable 

Name 

Description 

1 

Type 

"CONF" 

2  --  Wave  heading  set  number  from  previous  FREQ. 

(default  =  1) 

NOTES : 

1.  These  data  are  required  to  get  solution  data  for  wave  headings  other 
than  the  first  one  given  in  the  foregoing  FREQ  solution  (see  Section  7.3). 
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7.4. 2 

ANCH^  BUOY,  LINE 

-  Component  requests 

Word 

Variable 

Name 

Description 

1 

CTYPE 

"ANCH",  "BUOY",  "LINE",  ' 

’DONE" 

2 

NELD 

Node  or  element  number 

3 

ICODE 

Component  code 

4 

JCODE 

Bottom  factor  code 

.S 

CAPY 

Component  capacity 
(default  =  use  inventory 

with  SAFAC) 

6 

CMPID 

Component  Identifier  for 

Anchor  weight 

Chain  size 

Line  diameter 

Buoy  O.D. 

Inventory 

(input  in  inventory  units) 


7  SAFAC  Safety  factor  (default  =  1.0) 

NOTES; 

1.  CTYPE  identifies  the  type  of  component  to  be  checked.  Only  anchors, 
buoys,  or  lines  are  allowed.  Any  number  of  cards  can  be  given.  Input 
is  terminated  by  CTYPE  =  "DONE". 

2.  Negative  ICODE  with  CTYPE  =  ANCH  means  the  vertical  load  will  be 
used  in  the  capacity  check. 

3.  The  CHEK  SAO  can  be  initiated  following  any  other  subanalysis.  If 
the  previous  one  was  not  FREQ,  the  present  state  is  evaluated.  If  it 

was  FREQ  and  the  dynamic  response  file  was  generated,  the  dynamic  tensions 
will  be  included  in  the  check  for  those  elements  of  the  file.  (See  Note 
2  of  Section  7.3.1,  FSOL  record) 

4.  If  CTYPE  is  "ANCH"  or  "BUOY",  then  NEED  is  the  node  where  the  anchor 
or  buoy  is  located.  Anchor  checks  can  be  made  for  fixed  nodes  as  well 
as  for  active  nodes  where  limits  have  been  defined. 

5.  No  more  than  ten  lines  can  be  attached  to  any  anchor  or  buoy  at  an 

active  node  with  limits  (see  Note  3,  Section  6.12,  l.LOC  record).  Anchor  checks 
at  fixed  nodes  allow  up  to  twenty  lines  to  be  attached  to  the  fixed  node, 
which  has  no  limit  conditions  specified. 
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CompoiitMil  codes; 


ANCHORS 

1.  Navy  standard  stookless  with  stabilizers 

2.  NAVSHJHS  lightweight 

3.  NAVFAC  STATO 

A.  Imbedment  (no  inventory) 
b.  Stock  (Admiralty)  (no  inventory) 

6.  Mushroom  (no  inventory) 

B^£YS 

1.  Bar  riser  chain 

2.  Spherical  or  other  (no  inventory) 

LINES 

0.  Chain 

1.  Samson  2-in-l”  braided  nylon 

2.  Samson  2-in-lg  power  braid 

3.  Samson  2-in-K  stable  braid 

4.  Samson  Bine  Streak 
f).  Other  (no  inventory) 

Bottom  factor  codes  for  anchors: 

1 .  Compacted  sand 

2.  Stiff  dense  clay 

3.  Sticky  clay  of  medium  density  (cohesive) 

A.  Soft  mud  (fluid),  loose  coarse  sand,  gravel 

5.  Hard  bottom  (rock,  shale,  boulders) 

Contents  of  the  inventory  are  presented  in  Appendix  D. 
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7.5  ^JD>  NKW  '  Problem  termination 


The  fc!ND  flag  record  signals  the  completion  of  the  SAC  sequence  with 
no  more  data  to  be  read.  The  run  is  terminated. 

The  NEW  flag  record  signals  the  completion  of  the  SAO  se<iucnce. 

The  next  data  record  is  a  title  card  to  begin  a  new  problem.  This  title 
card  must  use  the  character  string  delimiters,  and  only  one  title  record 
is  allowed. 


8.0  KXAMPLE  PROBLEMS 


This  section  coiitaLiis  four  example  problenis  and  the  description  oi 
some  optional  approaclies  to  solving  them.  Althoiigli  not  exhaustive, 
these  problems  demonstrate  many  of  the  features  of  SEADYN.  The  models 
have  purposely  been  kept  small  and  geometerical ly  simple  to  make  them 
convenient  to  handle.  The  user  should  be  aware  that  SEARYN  is  capable 
of  modeling  much  more  complicated  geometries. 

® 1  Towed  Body  Example 

A  spherical  body  towed  at  the  end  of  a  280-foot  line  will  be  used  to 
demonstrate  various  methods  for  starting  the  problem,  the  effects  of  grid 
coarseness,  and  changing  a  quasi-static  solution  to  a  dynamic  solution. 
The  pertinent  problem  data  are; 

Tow  Cable  Data 

Unstretched  length .  280  ft  =  3,360  in. 

Drag  diameter . 0.35  in. 

Cable  weight  (in  seawater)  .  0.169  Ib/ft  ~  0.014  Ib/in. 

EA . 1.92  X  10^  lb 

-1/2 

Normal  drag  coefficient,  ....  1.5  +  4.0  (R^)  ' 

Tangential  Drag  Coefficient,  C.,.  .  .  0.2  C,, 

I  w 

Body  Data 

Weight  (in  seawater)  .  580.9  lb 

Effective  diameter . 1.0  ft  =  12  in. 

The  first  part  of  the  problem  is  to  determine  the  position  of  the 
cable  and  body  when  to'.<;d  in  a  straight  line  at  a  constant  velocity  of 
10.5  knots.  A  convenient  model  for  this  is  to  assume  the  tow  point  is 
fixed  and  the  system  is  subjected  to  a  uniform  current  equal  to  the  tow 
speed.  This  allows  the  static  solution  form  to  he  used. 

Although  the  general  shape  of  the  system  at  the  tow  speed  can  be 
guessed,  the  exact  positions  of  the  nodes  and  the  tensions  in  the  various 
portions  of  the  cable  are  not  easy  to  compute.  Three  approaches  to 
calculating  the  steady-state  towing  shape  are  outlined. 
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API’I\0A('I1  i:  i>y  lsMiiiiill>>  t  !l<'  lint-  i  •:  .Slr.ll.ullt 

down.  Api'ly  tho  ^;r.ivity  lo;id  v>  i  Ui  i«o  \<iii'eiu  to  i;i'i  t  ho  siispondcd 
lonisioii  (list  I'ihiil  ion.  I  no  roiiicnl  :i  1  i  y  apply  tin'  v'lii'ionL  to  do|)luy  I  lio 
cabU'  in  an  .i|ipio\iin,it o  yolution  (lUb  iiiot.honi.  lloialo  on  tho  dopioyo.l 
ytato  to  yatiBly  oqu  I  I  i  hi  i  tini  iiiolhod).  'I'ht'  in|int  loi-  aoooiiip  1  i  oh  i  no 

this  is  shown  in  Table  8--1.  It  r.hoiild  be  not  «'d  ili.ii  lin-  first  step 
requires  artiTio.il  tensions  to  j:et  st.irtoo.  Nuimiii'.il  damp  i  iip,  and  the 
MNR  solution  work  well  here.  Kither  VK  nioHiod  could  ai.se  In'  list'd. 

.M’l’ROACII  2-.  Ih'gin  hy  .issiiiii i iif,  the  ],iMe  is  deployed  horizontally 
in  an  unstretched  state.  .Apply  the  gravity  and  ciuit-'iit  simultaneously 
ill  Lncrenients  to  get  an  afiproximate  solution  (RKIl  method).  Dummy  tensions 
are  assumed.  Iterate  ou  the  deployed  .state  to  satisfy  etiui  I  i bri urn  (MNR 
method).  The  input  is  shown  iu  Table  8-2. 


Table  8-1.  TOWKl)  HODY  AlTKOACll  1  ..JN^UT 

TOWED  liODY  DEMONSTRATION  PRODLEIl 

APPROACH  1  ---  START  FROM  VERTICAL; 

PROB 

11, 10,  2,, 1,1,386 

FLUl 


0,0,  .002.'i2,  .037 
BODY 

1. , -580.9, 12 

MATE 

1.1.. 35..014112.1,W8,0,1,92E5.1 


NODE 

1..  .3360 
11,1,W6,1,1 ,1 
El, EM 

1.1.2., ],W8,1 
10  10  11, , 1 
B],0C ;  1 , 1 
FLOW 


*»OnY  NODE 

'■TOWF.D  NODE  IMPLIED  CENEKATION 
•'■ASSUME  IINSTRETCHED  LENOTHS 
•"LOCATE  BODY  AT  NODE  1 


1,1,  212.8 
DEAD 

SOLD ,,.001 

LIVE 

CDRR , 1 , 1 , ,  1 
S0I,U,RFB 
STEP, 10 

LIVE 

CIIRR,1 

END 


*DF.PL0Y  LINK  STRAIGHT  DOWN  WITH  MNR  t  DAMPING 
•>M0VE  TO  APPROX  TOWF.D  SHAPE  WITH  10  STEP  RFB 

*ITERATF,  TO  CORRECT  STATE  WITH  MNR 
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Tolilo  8-;;.  TOWKU  KOUY  Al’lMWAai  >’ 

TOWl'l)  nODY  OKNONSTKAIMON  i'ROIU.KN 

APl’KOACH  SiAKi  KKOM  HORIZONTA!,; 

I'ROl);  i  1 , 10.:;, ,  1 , 1 
FIAM  ;0,0,  .00252,  .0;)7 
1)ODY;1,0,-;>80.9,  12 
13LOC ;  1 , 1 
MATK 

1  , 1 , .  35 ,  .  0 1 A 1 1 2 ,  1 , Wft ,  0 , 1 .  ')2K5 , 1 
NODK 

1 . .  3360 

11, W6, 1,1,1 
LINE  ;9, 1,11 
ELKM 

1,  l,2,,l,W8,l 

10.10.11 . . 1 
TENS;1, 10,1, ,500 
FLOW; 1,1,  212.8 
I.IVK 

CIIRR,1  ,1 
S0Ll),RFB 
STEP,  100, 10 
LVAR,W.5,1 

LIVE 

S0LII,VRR 
CORK , 1 

END 


■^AS.SIJMK  UNSTRETCHEl)  LENGTHS 

■VAPPLY  DUMMY  TENSIONS 

•'MOVE  TO  APPROX  STATE  Wl'l'll  RFB 

'■SET  GRAVITY  LOAD  FLAG 


APPROACH  3:  Hogiii  l>y  guessing  an  angle  Lluil  approximates  the 
deployed  stale.  Deploy  tlie  oust  retched  cable  .at  that  angle.  Iterate  on 
this  guessed  slate  lo  satisfy  ecpii  1  ibrium.  Input  I'oi'  this  approach  is 
given  in  Table  8-3.  The  VRR  method  is  used  in  ihi.s  case  since  the 
quality  of  the  initial  guess  is  very  low  and  divergent  MNR  itecations 
can  he  expected.  'I'he  VRR  melltod  automatically  cointHuisat  i*s  for  t  uc  lack 
of  initial  tension. 


The.se  three  approactu’s  are  not  the  only  possihli'  c>pi  ioti;  .  In  (-.(ch 
approach,  the  beginning  state  was  descril)ed  by  the  positioii:;  of  ih;'  nodes 
in  an  unstrelcht'd  state  (ISIGO  '  1).  Other  approaches  c  ui  i'>  dev  i. o' 
quite  easily. 

The  luimher  of  elements  for  t  liese  tests  was  chosen  somew’liat  e.rhi  t  rari  ly . 
The  results  of  changing  the  niimbei'  of  eleimnits  ar«'  sntnnia  r  i  zed  in  I'ignrc  8-1. 
The  figure  sliows  quil<>  accurate  results  even  for  very  crude  model.s. 

These  models  have  assumed  uniform  spacing  ol  tlu'  nodes  foi  input  convoniein  o. 
3'h<-'  equilibrium  state  suggests  it  would  he  more  appropr  t  .it  «*  lo  put 
shorter  elements  neai'  the  lowed  body  wlu'ie  the  curv.it  lire  is  the  gre.iti'sl. 
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Th*'  St  o;i(ty-s(  )U  lowing*  roni  i  gttral  ion  fill)  l'(t  tisfMi  ;is  the  start  iitj; 
point  lor  it  (lynamir  analyais  in  whirli  the  toW  velocity  and  diioction  are 
eltanged.  Tal'le  8-A  shows  the  nddit  ioii:il  ininil  retpiired  to  contiiitte 
t  rom  any  el'  the  previous  eotiveixtod  static  :;tates  to  a  4tynaiiiie  soltilioit 
wltere  the  tow  point  is  movtMl .  Tlie  k«‘y  here  is  the  cfi:inp,e  in  flow-lield 
vt'loeity  to  /’.oio  and  assip.ning  the  initial  towing  velocity  to  all  oi  the 
nodes.  The  movonvt'iit  Itinction  given  del'itiest  a  sviddeii  cliange  of  dirttetion 
of  the  towitig  at  the  same  veloeity.  Arbitrary  movement  can  b«^  defined. 

Ta  1)  1  e  8  -  .1 .  TOWKD  BODY  APPROACH  .T 

TOWED  fiODY  DEMONSTRATION  PKOHI.RM 

APPROACH  3  START  KKOM  SCOPED  LINE; 

PROhjll.lO,  2,,  1,1, 386 
FLUlj 0,0, .00252,  .037 
1U)DY;1 ,0,-.580.9,12 
fiL0C;l,l 
MATE 

1 . 1 . .  35  ,  014112 , 1 ,W8 ,0 , 1 . y2E5 , 1 
NODE 

1..  31 95.. 55. 1038. 297 

11, W6, 1,1,1 

LINE; 9, 1,11 
El, EM 

1 . 1 . 2 . .  1 , W8 , I  ^ASSUME  UNSTRKTCHED  LENGTHS 

10.10.11., ! 

FLOW;!,],  212.8 

LIVE  A  iterate  on  GUESS  TO  CONVERGENCE  WITH  VR 

S0H),VRR 
CURR.l 
END 


Table  8-4.  TOWED  BODY  DYNAMIC  INPUT 

TFUN 

1,3, 1,1, 1,1,. 866, 90,. 866 
2, 3, 0,1, 0,1, .5, 90, .5 

DYN 

INIT,-212.8 

MOVE.ll ,2,1,-212.8,W9,2,2,212.8 

TIME,, 30  AEiNl)  OWN  TIME  STEP 

OUTP,W3,l 
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One  general  comment  is  needed  on  tiris  problem.  Care  mu.st  !)e  taken 
to  assure  all  data  are  provided  in  cciisistenl  units.  The  clioire  oi  units 
is  arbitrary  since  SKADYN  does  not  presume  .any  units  (except  for  certain 
identified  defaults).  The  clioice  made  here  is: 

Length  -  inches 
Force  -  pounds 
Time  -  seconds 

Note  that  lengths,  diameters,  forces,  velocities,  fluid  properties,  and 
gravitational  acceleration  all  use  these  units.  In  all  of  the  approaches, 
the  tow  cable  drag  coefficients  are  given  by  the  DRAGCO  subroutine  shown 
in  Table  8-5. 


Table  8-5.  DRAGCO  SUBROUTINE 

SUBROUT  I NE  DRAGCO  (I  ND''  ,  I  TYPE ,  NUM ,  REN ,  RET ,  CN ,  CT ) 
CN=0. 

CT=0. 

IF  (REN.LT.IOO. )  GO  TO  500 
GO  TO  (100, 200, 300). I TYPE 
C  SPHERICAL  BUOY 

100  CN=.A7 

GO  TO  500 

C  CYLIUDRICAL  BUOY 

200  GO  TO  300 
C  CABLE 

300  CN=1.5+4./SQRT(REN) 

CT=0.02*CN 
500  RETURN 
END 


Buoj'  Rc'l. 1X0 lion  Kxainplo 

ConsicitM  ,!  huoy  rrstjoined  liy  a  singio  line  tliat.  i  b  ;<miio  !'(.■(!  to 
t.lie  1  toor .  Af;siini<’  Ihal  Ihe  bitoy  ix  snagged  and  lioiiieii  ov-  r  to  a 
;('rto;in  posiljcfi  and  held  tbere.  Then,  alter  seme  time,  the  inioy  and 
I  itie  are  releasfHi  to  move  dyiniinica  1.1  y  to  the  ortginoi  gravity  lo.-ded 
state.  A  SF.ADYN  analysis  of  this  scenario  is  presented  n  three  .stages: 

STAGE  J  “  Establish  the  gravi ty- 1  o.aded  initial  .state. 

STAGE  2  -  Move  the  buoy  over  to  the  restrained  (.snagged) 
position . 

STAGE  3  -  Kelea.se  the  restraint  .and  follow  the  dynamics  back 
to  the  initial  state. 

The  actual  execution  of  these  stages  is  dependent  on  how  the  snagged 
position  is  defined.  The  problem  statement  implies  the  dynaiiucs  of  the 
snagging  operation  are  not.  iiriportant  and  only  the  calcul.ition  of  the 
firical  static  state  in  the  restrained  position  is  needed.  Two  d.ifferc-nt 
definitien.s  of  the  snagged  state  are  pursued  here.  The  first  presumes 
the  snagged  position  of  the  buoy  is  known.  The  second  presumes  the 
magnitude  and  direction  of  the  snagging  load  is  specified. 

Recalling  the  previous  example  problem,  there  are  various  approaches 
that  can  be  used  to  get  the  initial  static  state.  In  both  cases  a 
gravity-load  vertical  state  is  used  as  a  starting  point.  Loads  or 
displace ments  are  then  applied,  and  iterations  are  carried  to  convergence 
at  the  final  .state.  The  dynamic  phase  is  initiated  simply  by  releasing 
tile  displacement  constraint  or  force.  Input  for  these  examples  is 
presented  in  Table  8-6.  The  second  case  is  assumed  as  a  restart  from 
tire  fir.st  to  demonstrate  the  RESTART  section.  The  vi.itcous  relaxation 
method  was  chosen  for  both  of  these  since  it  is  most  likely  to  remain 
stable  and  obtain  convergence. 

Tlie  dynamic  solution  used  the  direct  iteration  method  since  this 
has  proven  to  be  the  most  reliable.  The  algorithm  was  left  to  choose 
the  time  step  needed  since  no  external  loa<ling  with  its  own  time 
characteristics  wa.s  present. 
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Table  8-6.  BUOY  RELAX-VriON  INPUT 


BUOY  RELAXATION  DEMONSTRATION  PROBLEM  —  INITIAL  POSITION  GIVEN; 


PROB;U,10,-2,  ,1 
FLUI;0,I 
BODYjl,  ,.129,2 
BLOC ; 1 , 1 

MATE ; 1 , , . 15 , . 00055 ,W9 , 2 . 36 , 1 
NODE 


1..  .5.551 

11, 1, 0,. 104, W6, 1,1,1 
ELEM 

1.1.2. . 1 
10,10,11,  ,1 
DEAD 

SOLD, , .01 
SAVE,-1 

LIVE 

SOLU.VRR 
MOVE, 1,1, 1,4, 1,1 


*1SIGO=0  WITH  NO  loads  SAME  AS  ISlGO^l 

-GET  VERTICAL  STATIC  STATE 
*SET  NUMERICAL  DAMPING 
-SAVE  DEAD  RESULTS 
-'MOVE  TO  SNAGGED  POSITION 


'2 

*RELEASE  FOR  RELAXATION 


TIME,, 10 
OUTP,  ,1 

NEW 

"BUOY  RELAXATION  DEMONSTRATION  PROBLEM  -  FORCES  SPECIFIED" 

REST ; NEW, 1  *GET  THE  LAST  DATA  SAVED  FROM  PREVIOUS  DEAD 

LIVE  *APPLY  SNAG  LOAD 

SOLU.VRR 

LOAD, 1,. 115,-  0265,,!  -'LVAR  DEFA” 

DYN  *RELEASE  FOR  REL' ■  ON 

TIME, ,10 
OUTP ,  ,  1 
END 


8.3  Anchor  Last  Example 

An  anchor  Iasi  deployment  scheme  for  tethered  buoys  is  readily 
modeled  In  SEADYN .  The  starting  state  before  the  release  of  the  anchor 
requires  some  careful  con.sideration ,  however,  when  special  situations 
are  modeled.  Static  analysis  problems  can  occur  when  floating  lines 
with  low  initial  tensions  are  involved.  This  presents  a  problem  bccrause 
slight  changes  in  tension  distribution  can  cause  large  position  changes 
In  such  situations  it  is  best  to  make  a  reasonable  guess  on  the  initial 
state  with  no  tensions  and  go  immediately  to  the  dynamic  analysis  following 
the  release  of  the  anchor. 

This  is  a  relatively  simple  problem  intended  to  demonstrate  the 
mechanics  of  getting  the  static  .ind  dynamic  solutions  accomplished.  The 
problem  presented  here  is  sufficiently  well-conditioned  to  allow  the 
initial  static  state  to  be  calculated.  The  problem  characteristics  are 
given  in  Figure  8-2.  The  input  data  are  given  in  Table  8-7.  The  input 
assumes  the  anchor  is  released  at  time  zero,  and  the  top  buoy  is  held 
for  10  seconds  before  it  is  released.  The  line  between  the  two  buoys  is 
assumed  to  be  neutrally  buoyant.  The  initial  state  holds  node  1  fixed 
and  applies  1,000  pounds  horizontal  force  to  tiode  10.  This  stretches 
the  neutrally  buoyant  line  out  straight  and  produces  a  catenary  configu¬ 
ration  in  the  heavy  line.  The  guessed  initial  input  uses  uns*^  retched 
lengths  (ISICO  =  1)  and  the  catenary  line  generator.  The  initial  distance 
between  nodes  4  and  10  is  selected  to  give  the  desired  line  length  for 
the  catenary. 


Table  8-7.  ANCHOR  LAST  DEPLOYMENT  INPUT 


DOUBLE  BUOY  DEPLOYMENT  USING  ANCHOR  J.AST  METHOD; 

PROB;10,9,1,1,1 

FLUI;0,1 

LIMI 

1 . ,  1  *SURFACE 

2.5000.1., !  *BOTTOM 

BLOC 

1,1,W5,1  ••BUOY  I  AT  NODE  I 

2,4,W5,1  *ni]OY  2  AT  NODE  4 

3,10,W5,2  * ANCHOR  AT  NODE  10 

NODE 

1,W6,3,3,3  'WILD  TOP  BUOY 

4  ,  ,  ,  1200  ,  ,3  ^VERTICAL  HOLD  ON  BUOY 

10.. . 2400.. 3.,  1  '■^HOLD  ANCHOR,  LEAVE  HORIZONTAL  FREE 

LINE;5,4,10,W7,2,.r),  1000  *GENERATE  €.4 TENARY  .AND  NODES  2,3 

2,1,4 

ELEM 

1.1.2.. 1,W8,1,1000 
^,A„5,,2 

9.9.10. . 2 
MATE 

1. . .3,W9,1.K4,1 

2.. .2..5,W9,1.E4,1 
BODY 

1 ,  ,2000,8, W8,20  'BUOY  1  WITH  SURFACE  DRAG 

2,  ,1000, 6, W8, 15  -BUOY  2  WITH  SURFACE  DRAG 

3.,  -5000, 5  ^ANCHOR 

DEAD  *S INGLE  STEP  VRR  SOLUTION 

SOLU.VRR 

LOAD, 1,,  1000,, 10 

v-DROP  ANCHOR 

FIX, 1,41  ^CHANGE  CODE  TO  ALLOW  BUOY  PULL-UNDER 

FREE, 101 
TIME,  .10. 

OUTP,,2 

DYN  '''KELE.ASE  BUOY  1 

FREE  ,11,12 
TIME,, 1000 
OlJTP,,iOO 

END 


8.4 


^ '‘2DL.^  Moving  Ship  Example 

The  deployment  of  a  cable  system  by  paying  out  line  from  a  moving 
.ship  is  demonstrated  in  this  example.  The  problem  selected  is  the  final 
stages  of  deployment  of  a  trapezoidal  array  system  (Ref  15).  Figure  8-3 
shows  the  starting  configuration.  The  system  contains  tVi;o  clump  anchors 
and  two  suspension  buoys.  In  the  sequence  modeled,  the  first  cliutip 
anchor  is  in  place  and  the  final  anchor  is  being  lowered  by  a  line 
paying  out  from  a  cable  ship.  The  payout  rate  and  ship  velocity  are 
chosen  .so  as  to  pi  ace  the  anchor  at  the  desired  point.  The  demonstration 
problem  tises  a  straight  line  ship  velocity  (away  from  the  cable  system) 
of  1  ft/sec  and  a  constant  payout  rate  of  2  ft/sec. 

The  element  model  used  is  the  simplest  possible:  only  one  element 
per  line.  The  payout  line  element  is  assumed  to  grow  to  a  maximum 
length  of  10,000  foot  and  then  be  divided  into  two  elements  (mitosis). 

Up  to  this  point  element  5  and  node  6  are  not  involved  in  the  solution. 
When  mitosis  occurs,  node  5  is  placed  at  the  dividing  point,  node  6 
takes  the  position  of  the  moving  payout  point,  and  element  5  enters  the 
system  as  the  upper  half  of  the  payout  line  with  an  unstrotched  length 
equal  to  the  mitosis  length.  Element  5  is  then  the  payout  element,  and 
element  4  maintains  a  constant  unstrotched  length  equal  to  the  unstretched 
length  just  before  mitosis  minus  the  mitosis  length.  When  the  anchor 
(node  4)  reaches  the  bottom,  it  is  held  fixed  in  three  components. 

Table  8-8  lists  the  input  for  a  DYN  analysis  using  the  DIM  method. 
Table  8-9  presents  input  for  a  TSSS  analysis  of  the  same  time  sequence. 

The  VRR  form  was  used,  but  the  MNR  method  could  also  have  been  selected. 
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Table  8-8.  PAYOUT  FROM  MOVING  SHIP  INPUT 

PAYOUT  FROM  A  MOVING  SHIP; 

PROB;6,5,2,l,l 
FLUI; 14400,1 
BODY 

1. . 3280.8 

2. . 2160.8 

3..  5120. 5 
BLOC 

1,2 

2  3 

3 1 4 , W5 , 1  *STOP  ANCHOR  AT  BOTTOM 

LIM1;1,0,5,,1 

MATE; 1,,. 1667 ,W9,1.E6,1 

NODE 

1,W6, 1,1,1  *GLOBALLY  FIX  ANCHOR  NODE 

2..  7500. 13300 

3.,  10400, 14200 
4,  ,23000,5700 

5.. 29739.14400,W6,2,2,2  *ORIGINAL  PAYOUT  NODE 

6. . 29739. 14400, W6, 1,1,1  *NEW  PAYOUT  NODE 
ELEM 

1.1.2..  1 

5.5.6. . 1 
DEAD 

SOLU.VRR 

DYN 

MOVE,5,2,,l,l,W9,l 

PAyO,5,4,2,5000,l,l 

S0LU,W12,-1  *RESTRICT  TIME  STEP  GROWTH 

TIME,. 5, 3600 
OUTP,,100 
END 


Table  8-9.  TSSS  INPUT  FOR  PAYOUT  PROBLEM 


TSSS 

SOLU,VRR 

TIME, 100,4000 

MOVE,5,2,,l,l,W9,l 

PAYO,5,4,2,5000,1,1 

OUTP,,100 

END 
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9.0  HELPFUL  HINTS  FOK  MODEL] N(5  AND  TKOllBLESHOOTl NO 


Th(!  most  oonmioii  problem  encountered  by  a  i.ew  user  of  a  complicated 
coinputer  program  is  the  mastery  of  the  injnit  syntax.  It  is  hoped  tli.il. 
the  example  problems  will  alleviate  a  large  number  of  tlicsc  problems. 

An  at.lempt  to  tormulale  SKADYN  input  for  similar  problems  sfiould  be  made 
before  going  on  to  the  user's  own  probD.'ins.  A  variety  of  small  pioblcms 
should  be  tried  until  the  syntax  is  relatively  familiar. 

Beyond  the  syntax  problems,  one  encounters  the  modeling  errors  and 
the  need  for  understanding  the  quirks  of  tlte  solution  methods.  Finite 
element  modeling  requires  that  certain  abstr.actions  and  approximations 
be  made  to  get  from  the  real  world  to  the  numorical  model.  One  must 
decide  how  many  no<les  to  use,  where  to  place  them,  how  to  number  them, 
what  to  lump  where,  etc.  Since  SEADYN  uses  straight-line  elements, 
tlu're  will  always  be  a  disparity  between  the  actual  length  and  the 
modeled  length  on  curved  lines.  Also  there  will  always  be  the  problem 
of  step  changes  in  tension  from  one  element  to  the  next.  These 
approximations  are  inherent  in  the  modeling  process,  ai>d  the  analyst 
needs  to  become  familiar  with  the  consequences  of  approximating  for  his 
particular  situation.  Answers  to  questions  of  how  many  nodes  and  where 
to  place  them  are  strongly  problem-dependent.  In  general,  enough  nodes 
should  be  provided  to  model  geometry  and  tension  variations  to  the  level 
that  is  important  to  the  problem.  Each  element  assumes  a  uniform  material 
so  each  material  must  be  represented  by  at  least  one  element.  Highly 
curved  regions,  and  regions  where  tensions  can  vary  significantly  require 
shorter  elements.  How  short  is  short  enough  can  be  answered  from  grid 
convergence  studies  on  a  particular  problem.  The  result  presented  in 
Figure  8-1  is  an  exantple  of  such  a  convergence  study. 

The  solution  of  simultaneous  equations  in  SKADYN  follows  a  Gai):sian 
elimination  procedure  that  attempts  to  minimise  computer  storage  by 
taking  advantage  of  equation  symmetry  and  handedness.  This  means  that 
the  amount  of  storage  required  and  the  number  of  matrix  operations  arc 
sensitive  to  the  way  the  nodes  are  numbered.  As  a  general  rule,  the 
largest  difference  between  node  numbers  defining  any  element  determines 
the  bandwidth.  The  smaller  the  difference  the  smaller  the  bandwidth  and 
the  smaller  the  cost.  In  calculating  the  maximum  node  number  dilfereiue, 
the  master  node  number  or  the  master  node  number  pins  one  (since  node 
pairs  are  used  to  get  rotation)  must  be  used  in  pi. ice  of  the  slave 
nodes.  Another  general  rule  for  node  numbering  '  ■'ulioned  in  Section 
3.0)  is  to  number  such  that  the  softer  (lower  si  Incss)  compoiUMils  are 
encountered  first.  This  reduces  the  numerical  error  propagation  in  the 
analysis.  Often  it  is  not  obvious  what  components  are  soft  relative  to 
others  and/or  requirements  for  meeting  this  condition  are  in  conflict 
with  the  optimum  bandwidth  numbering.  Judgment  and  experience  on  the 
part  of  the  analyst  are  needed  in  making  compromises  in  this  area.  It 
should  be  noted  that  all  solution  methods  in  SEADYN,  except  the  DIM 
method,  require  the  solution  of  simultaneous  equations  and  are  bandwidth 
sensitive.  The  MODE  analysis  uses  the  Jacobi  iteration  method  and  is 
also  bandwidth  sensitive. 
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Krrors  Hint  rnii  occur  bocauso  of  node  numbering  can  show  up  n.s  a 
singular  system  of  equations.  The  message  printed  would  lie  "EQlJAiTON 
DECOMPOSITION  FAll.ED,  etc,"  (see  Section  10.3).  This  message  also 
results  from  problems  with  input  errors  and  for  unstable  structures,  so 
out’  should  not  be  hasty  in  concluding  it  is  the  result  of  accumulated 
numerical  error.  A  svHicific  situation  where  this  can  occur  is  in  a 
multi  leg  mooring  of  a  ship.  Certain  pretensions  and  positions  of  t  he 
legs  can  cause  the  ship's  responses  to  be  stiffer  than  the  lower  portions 
of  the  mooring  legs.  If  the  ship's  nodes  are  numbered  before  the  mooring 
legs,  numerical  error  propagation  can  cause  a  singular  matrix  error  at 
the  lower  ends  of  the  lines.  Renumbering  the  nodes  so  that  the  ship  is 
last  will  remove  the  problem. 

A  very  common  error  leading  to  a  singular  stiffness  matrix  is  an 
Lncoinplel.e  set  of  boundary  conditions  or  an  inconsistent  set  of  preten¬ 
sions  that  causes  a  portion  of  the  structure  to  go  slack.  Boundary 
conditions  must  be  given  to  fully  restrict  rigid  body  motion  of  the 
structure.  This  means  six  degrees-of-f reedom  must  be  constrained. 

Since  SEADYN  treats  geometric  nonlinearities,  the  rotation  degrees-o f- 
frcedoin  represented  by  elements  displacing  one  end  laterally  relative  to 
the  other  can  he  restrained  by  tension  in  the  clement.  For  example,  a 
single  element  held  fixed  at  one  end  (three  constraints),  with  a  tensiht 
load  at  tlu'  other,  is  a  stable  system  if  the  stretch  in  the  element  is 
ill  balance  with  the  load.  The  tension  in  the  element  resists  the  pivoting 
of  the  element  about  the  fixed  end.  There  is  no  twist  degree-of-freedom 
in  the  element,  so  no  constraint  is  required  for  it.  If  the  end  load  is 
not  applied  in  line  with  the  element,  a  rotation  must  occur  until  the 
('lement  aligns  wi  til  the  load. 

Section  ,3.0  disctissed  the  problem  of  obtaining  a  stable  initial 
configuration  of  the  finite  element  model  to  be  used  as  a  starting 
point.  The  analyst  should  always  keep  in  mind  that  the  finite  element 
model  is  an  approximation  of  the  actual  structure.  Even  if  exact  data 
describing  an  initial  static  state  of  the  actual  .structure  were  available, 
use  of  that  data  in  a  finite  element  model  would  not  represent  a  converged 
equilibrium  state.  Some  minor  adju.stments  in  position  of  the  nodes  and 
element  tensions  would  inevitably  be  calculated  by  the  numerical  iteiMlioii 

Two  procedures  for  getting  stable  initial  configurations  have  lieeu 
used  with  some  success.  The  examples  demonstrate  both  of  them.  The 
first  apiiroach  is  to  use  a  inimerically  damped  MNR  method  or  a  combination 
of  a  few  Ill'll  Steps  followed  by  MNK  iterations.  The  second  approach  is 
to  use  the  VRll  method  with  or  without  increments.  In  either  of  these 
approaches  It  may  be  nece.ssary  to  adjust  the  analysis  p,i raniot  ers  to  got 
convergence.  Two  general  forms  of  behavior  can  lead  to  a  lack  of  conver- 
gence ,  very  slow  convergence,  or  wildly  oscillating  behavior. 

If  d.iinping  is  too  large  for  the  particular  situation,  movement  will 
be  sluggish,  and  the  force  residual  will  not  change  significantly  in  any 
one  iteration.  This  can  occur  in  the  VRR  solution  if  the  structure  is 
stiff  (c.g.,  heavily  loadtnl,  high  KA  lines).  The  VRR  solntion  will  sig¬ 
nal  this  situation  by  repeated  out.put  of  the  SLOW  CONVEllCiKNCE  mess.ige 
with  neither  of  the  TNl'llEASING  NORM  mess.iges.  The  .ippropriate  action  in 
this  case  is  to  try  the  analysis  again  with  a  much  smaller  initial  num¬ 
erical  damping  (Word  (3)  on  the  SOLD  record).  The  default  value  is  0,001, 
but  in  some  situations  it  may  be  necessary  to  make  it  much  smaller. 
Alternatively,  the  MNR  method  could  bo  tried  with  default  ))arameters. 
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L 

I  Tile  vrt'y  slow  eoiivcrgenoc  behavior  of  the  VUH  method  raiRp.es! tliat  the 

j  algoritliiii  wauls  to  adapt  near-zero  damping,  whieh  is  the  Newt  on-Raphsoii 

i  iiK'thod.  Slii{’.f>ish  behavior  of  the  MNK  inclhod  also  i)ie<iiis  luiinerieal  damping, 

i:  is  too  ltlj>;h. 

i  The  more  eoinnion  silnatioii  is  the  divereenl  behavior  eviileiui'd  by 

i  very  large  (olteii  increasing)  oscillating  results  I  rom  tlie  iteiatioiis, 

i  This  usually  means  large  damping  should  ho  used  or  strategy  modifications 

made  vvliich  are  suggested  by  the  pattern  of  the  unstable  behavior. 
Divergence  is  most  often  encountered  using  the  MNR  method  when  a  poor 
starting  point  is  given  that  requires  a  significant  angular  movoment , 

The  VRR  approach  will  generally  work  much  hotter  than  the  MNR  approach 
in  these  (ases  since  it  can  more  readily  adapt  to  the  problem  by  damping 
and  step-size  adjustments.  The  adaptive  features  of  llu'  VRR  melbod  have 
been  dt'vclopcd  somewhat  empirically  by  following  a  particular  divergent 
pattern  and  then  constructing  a  strategy  to  detect  and  react  to  it . 

This  means  it  is  possible  that  tlie  action  taken  is  inappropriate  in  some 
new  situations  or  Vliat  a  nonconvergenl  behavior  with  a  different  pattern 
is  undetected.  These  situations  usually  produce  copious  nies.sages  with  a 
mixture  of  INCREASING  NORM  and  SLOW  CONVERGENCE  forms,  and  no  clear 
pattern  of  reduction  of  the  force  residual  is  obtain'd.  The  solution 
j  will  then  fail  by  running  to  the  limit  number  of  iterations  (Word  (11) 

ij  on  SOliU  record)  or  the  job  limits  (time  or  output  lirres).  When  the  VRR 

j  method  runs  to  the  limit  iterations  without  converging,  the  pattern  of 

!lhe  messages  and  the  values  of  the  residual  and  velocity  nortiis  should  be 
:  studied  before  attempting  a  rerun.  Sluggish  behavior  will  be  indicaled 

’,t  by  small  velocity  norms,  and  large  residual  norms  with  small  position 

I'  changes.  In  this  case,  the  actions  noted  above  should  be  consid(>red. 

i|  If  the  pattern  shows  large  velocity  norms  and  no  definite  residual  norm 

!  pattern,  an  increase  in  damping  could  be  effective.  Applying  the  loading 

in  increments  could  also  be  considered. 

If  very  erratic  norm  changes  occur  and  the  values  are  very  large, 

5'  one  should  review  the  problem  data  to  see  if  anythiug  is  wrong  there. 

1!'  Finding  nothing  wrong,  a  different  strategy  should  be  considered.  This 

I  can  involve  a  change  in  the  initial  state  input  and/or  a  different  mix 

of  analysis  parameters.  Very  abrupt  load/displacement  changes  should  be 
introduced  incrementally.  The  physical  possibility  of  unstable  states 
should  he  investigated  (it  can  be  that  the  load  level  requested  produces 
a  physical  instability). 

One  mistake  that  is  sometimes  made  in  dealing  with  the  staging  of 
subanalyses  is  the  sudden  omission  of  loads  or  constraints  imposed 
during  the  preceding  stage.  The  load  and  displacement  data  defined  at 
the  beginning  of  each  stage  nuist  bo  consistent  with  the  total  values 
I  obtained  at  the  end  of  the  preceding  stage.  For  ox.ainplc,  assume  Si  .age  A 

applied  load  A  to  the  structure  with  displacement  constraints  A.  At  the 
I  beginning  of  Stage  B,  the  loads  and  constraints  defined  must  be  load.s  A 

and  displacement  constraints  A  unless  it  is  an  initiation  of  a  dynamic 
sequence  with  suddenly  released  conditions.  If  Stage  R  is  intended  to 
begin  from  load  A  and  impose  additional  loads  B,  it  is  necessary  to 
redefine  loads  A  with  a  zero  variation  code  on  the  I.VAR  record,  and 
define  loads  B  with  the  +1  variation  code.  Stage  B  will  tlu'ii  maintain 
loads  A  and  incrementally  add  loads  B.  If  loads  A  were  omitted,  the 
initial  loading  would  not  be  in  balance  with  the  nodal  posit  ions  and 
I  element  tensions.  This  would  then  lead  to  an  abrupt  change  in  position 

I  and  possibly  a  divergent  solution  step. 
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MODK  aniilysrs  do  not  roquite  the  specif icalion  of  loads,  but  the 
d i sp.l aceiiUMit,  conslrai iiLs  in  force  at  tfie  end  of  the  |>rei;»‘d iny;  slaj'e  must 
remain  in  force.  If  these  displacements  w<*re  imposed  with  a  MOVI^  record, 
then  a  KIX  record  or  MOVK  record  would  he  re<iuired  for  these  components 
ill  the  MODK  SAO  set  . 

Dynamic  analyses  present  a  new  set  of  problems.  The  most  effective 
form  for  time  domain  dynamics  is  the  DIM  method.  Generally,  tlie  default 
parameters  will  he  cipj)ropri ate .  The  choice  of  a  fiiiie  step  is  the  main 
concern.  The  DIM  method  will  select  its  own  time  step  xf  none  is  given; 
however,  the  selection  algorithm  is  not  completely  reliable.  One  problem 
that  occurs  for  strongly  nonlinear  dynamics  is  that  the  appropriate  lime 
step  may  vary  with  time.  A  new  upper  limit  on  the  Lime  step  is  ealcnlalod 
only  when  there  is  a  signal  that  the  analysis  is  not  converging.  The 
signals  used  are: 

a.  A  displacement  change  norm  exceeding  1.F.12 

b.  A  large  ntimher  of  iterations  without  convergence  (LMITKR  on  the 
SOLD  card) 

c.  A  persistently  increasing  displacement  change  norm 

In  these  situations  the  time  step  upper  hound  will  he  recalculated  and 
the  step  .^ize  reduced.  It  is  possible  that  the  problem  will  not  be  de¬ 
tected  early  enough  to  avoid  fatal  growth  of  errors  in  the  analysis. 

Time  domain  analyses  which  make  repeated  time  step  changes,  should  be 
rerun  with  the  step  forced  to  remain  below  the  range  of  the  changes. 

This  xs  done  hy  specifying  l)T  (Word  (2)  on  the  TIMK  record)  and  setting 
JMPUT  to  -I  (Word  (12)  on  the  SOLD  record). 

Time  domain  analyses  without  material  damping  tend  to  be  sensitive 
to  stop  size  and  will  produce  spurious  oscillations  if  the  step  size 
and/or  convergence  tolerance  arc  not  appropriate.  Even  with  damping, 
erroneous  oscillations  can  be  induced  by  setting  the  convergence  tolerance 
(DKRR)  too  high.  The  correct  value  of  DKRR  is  somewhat  problem-dependent. 
The  default  value  of  0.001  is  a  middle-of-the-road  choice.  One  sho,ild 
not  set  it  much  higher  without  some  verification  that  the  results  arc 
acceptable.  It  m.xy  be  necessary  to  set  it  lower.  One  example  ot  this 
is  when  the  time  step  is  desired  to  he  much  smaller  than  the  upper 
hound,  and  strong  nonl  inear it ies  are  present  (c.g.,  payonl  or  bottom 
inter.act ion) .  There  is  a  trade-off  between  step  size,  DKRR,  and  the 
number  of  iterations  to  convergence.  Although  the  optimum  mix  of  slei>s 
and  iterations  is  hard  to  determine,  experience  has  indicated  that  DT 
and  DEKR  values,  which  give  between  three  and  .ihoiit  seven  or  eight 
iterations  per  slop,  are  reasonable  for  strongly  nonlinear  problems. 
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JO.O  SUMMARY  UK  MAJOR  DIACNOSTH'.  MKSSACKS 


Tlw;  SK.ADYN  p iot>,r;<ni  makrs  chocks  ol  the  input  liala  aiul  atti'iiipts 

to  aid  the  us«.'f  in  finding  J\is  crruis  by  i.'finling  v.irif>n;;  iiK'ssagcs .  No 
attcnipl  has  tM'cn  made  to  be  compi'chcns i v<-  in  this  fealiiit!  since  it  is 
v<?ry  difticull  to  foia-sce  .ind/or  detect  many  of  the  possilile  errors. 

The  input  routine!;  that  process  t.lie  I’ROB  and  SAG  dal.i  sets  produce 
various  diagnost  ic.s  tl'.at.  evaluate  errors  detec  ted  in  llic  in|iut  .  These 
messages  are  gent'rally  sel  f-<!xp lan.it ory.  They  d(!al  mniniy  with  program 
restri  etions ,  .such  as  the  maximum  Tuunber  of  items  aliowed,  or  the  complete¬ 
ness  and  coiisi.stency  of  the  data  provided.  The  ti.ser  should  hiiv«*  little' 
difliculty  interpret  tng  the  prohlem  detected,  .and  shoulcl  he-'  able  to  make 
appropriate  corrections  witli  the  aid  of  this  manna  I . 

Dttring  the  processing  of  lh<'  snban;»  I  yses ,  cln-iks  are  tiMdc'  on  the 
validity  of  the  requests  and  the  convergence  of  llic'  analysis  procedure.';. 
The  messages  that  e.-ut  he  printed  are  listed  below  witl;  a  brief  der.e  r  i  pt  i  on 
of  the  probable!  eanse  nnd/or  cute.  The  action  taken  after  the'  erroi 
det.ection  is  indicated  by  the  following  codes; 

(K)  Fatal,  run  aborted. 

(N)  Abort  analysis  case  and  .seek  .a  new  problem  definition 
by  searching  the  deck  for  a  NEW  .SAG  record. 

(O)  Abort  present  SAG  activity  and  go  to  the  next  SAG  record. 

(fi)  Skip  this  request  ami  go  to  liie  next  card. 

(C)  Continue  calculation  with  action  as  indicated. 
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10.  I 


SysltMii  LV'Scr  i  pt  i  on  (^hroks 

(K)  nF.rrii  cdKRKcr  1  ON  i.:ruor  --  ship  okakt  kxci-.i.ds  watkr  iii.i-rii 

(ii'oss  iiipiil  ('rear.  (Iliock  lor  HroppoH 

(N)  KLl'.MKN'l'  OKNl'.KAl'ION  KRUOK 

First  oloinont  was  not  input,  or  rnrOs  .ire  not  in 

incronsiiiR  olomeiit  tunnber  order,  or  last  elcMU'ni.  w.is  nof, 
input . 

(F)  KRKOK  IN  (WIND/CURRENT)  LOAD  TABLE  ON  SHIP  XXX 
HEADING  =  XXX 
LAST  TABLE  ENTRY  -  XXX 
SYMMETRY  FLAG  -  XXX 

Heading  re<piost<*d  exceeds  the  largest  value  iii  the  table. 

Check  ship  load  input  table, 

(N)  IMPROPERLY  DEFINED  MOORING  BUOY  .AT  NODE  XXX  NO.  OF  SL.AVES  FOUND  XXX 

Mooring  buoys  require  at  least  two  slaves. 

(N)  NODES  OMITTED  IN  GENERATION 

Message  is  (ollow<?d  by  a  list  of  ones  and  .:eros  (ti*n  per  1  use) 
corresponding  to  the  nodes  in  the  system.  The  zeios  indicate 
wbit'b  nodes  were  not  defined  eitlu'r  by  input  or  generation 
soqvienco.  All  nodes  must  be  accounted  for.  Following  the 
integers,  the  nodal  coonlinate.s  are  priiUeil  to  aid  in  finding 
the  error.  The  problem  usually  comes  from  improper  node 
geiuM'ai  ion  input  . 

(N)  SHIP  DATA  INPUT  ERROR  ON  XXX  (Ship  No.) 

NO.  OF  SHIPS  ON  FILE  =  XXX 
LOAD  FUNCTION  OPTION  =  XXX 

Blank  record  for  moored  ship  data  requested  definition  of  ship 
from  load  fill?  with  no  load  file  defiu,.vi, 

OR 

Atti'mpted  ship  scaling  with  no  load  file  defined. 

(N)  TAPE  POSI  TION  INC.  OR  FORMAT  ERROR 

unable:  to  find  ship  data 

ITEM  XXX  LAST  READ  IS  FOR  HEADING  XXXX 
AND  WAVEI.ENGTH  XXXX. 

WANTED  HEADING  XXXX  WITH  WAVE  LENGTH  XXXX. 
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The  ship  motion  file  is  not  formulated  properly  or  other 
input  or  equipment  malfunction  has  made  reaoing  the  file 
impossible . 

(C)  WARNING- -UN ITS  DO  NOT  APPEAR  TO  BE  CONSISTENT 
GRAVITATIONAL  ACCELERATION  FROM  TAPE  IS  XXX 
UNIT  LABELS  FROM  TAPE  ARE  XXX  XXX  XXX 

The  ship  motions  file  conversion  factors  do  not  properly 
convert  the  GRAV  on  the  file  to  the  GRAV  specified  in  this 
run.  Calculation  still  proceeds. 
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10.2  Subanalysis  Option  Errors 

(S)  ANCHOR  NOT  ON  BOTTOM  AT  NODE  XXX 
SKIP  REQUEST  FOR  XXXX 

Anchor  weight  is  not  sufficient  to  hold  the  lines  at  this 
node  and  it  has  been  lifted  off  the  bottom  or  has  not  reached 
the  bottom. 

(F)  BLANK  COMMON  ON  TAPE  LARGER  THAN  SPACE  AVAILABLE 

The  data  saved  on  the  file  requires  more  storage  than  is 
presently  available.  The  two  numbers  printed  are  the 
required  and  current  values  of  NCOM. 

(S)  CAPACITY  AND  COMPONENT  ID  BOTH  ZERO 
SKIP  REQUEST  FOR  XXXX  (CTYPE) 

Check  input  card. 

(C)  DATA  NOT  AVAILABLE  FOR  XXXX 
SKIP  THIS  REQUEST 

Regular  wave  response  data  were  requested  for  a  frequency 
outside  of  the  range  that  was  generated  on  the  Ship  Motion 
file. 

(F)  INCONSISTENT  ANALYSIS  REQUEST  ON  RESTART 

Attempted  a  restart  of  DEAD,  LIVE,  or  DYN,  and  the  file 
was  not  from  that  type  of  subanalysis. 

(N)  INSUFFICIENT  STORAGE  TO  PROCEED 
COMMON  SIZE  =  XXX 
STORAGE  NEEDED  =  XXX 
HALF  BANDWIDTH  =  XXX 
DEGREES  OF  FREEDOM  =  XXX 
BASE  SIZE  =  XXX 

Subanalysis  request  (DEAD,  LIVE,  DYN)  cannot  be  processed  due 
to  storage  limitations.  Problem  must  be  reformulated  or  NCOM 
increased.  (See  Appendix  J.) 

(N)  INVALID  CALCULATION  OPTION  =  XXXX 
CASE  TERMINATED 

The  frequency  domain  calculation  option  was  not  RAND,  REGU,  or 
DONE.  Records  are  out  of  sequence  or  record  mispunched. 
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(N)  INVALID  COMPONENT  TYPE  XXXX 
ASSUMED  OF  INPUT 

Occurs  when  random  response  requests  arc  being  processed  and 
a  card  is  encountered  which  does  not  have  SHIP,  NODE,  TENS  or 
DONE  keyword.  The  items  to  this  point  are  processed,  and  l.he 
case  is  aborted  with  the  additional  message. 

(S)  INVALID  COMPONENT  TYPE 
SKIP  REQUEST  FOR  XXXX 

The  component  number  is  not  one  recognized  by  the  component 
inventory. 

(C)  INVALID  NODE  NUMBER  =  XXX 
SKIP  THIS  REQUEST 

Regular  wave  response  data  were  requested  for  a  node  number, 
which  is  less  than  1  or  greater  than  the  number  of  nodes 
in  the  model . 

(S)  MORE  THAN  TWENTY  LINES  ON  ANCHOR  AT  NODE  XXX 
SKIP  REQUEST  FOR  XXXX 

The  fixed  node  where  anchor  capacity  check  is  requested  has 
too  many  elements  connected  to  it  to  make  the  check. 

(S)  NO  DYNAMIC  TENSION  PROVIDED  ON  ELEMENT  XXXX 
IGNORE  DYNAMICS 

The  random  response  data  for  this  element  was  not  requested  by 
a  TENS  record  for  this  wave  heading  in  the  FREQ  SAO. 

(S)  NO  LINES  CONNECTED  TO  NODE  XXX 
SKIP  REQUEST  FOR  XXXX 

Check  node  number. 

(N)  NOT  ENOUGH  STORAGE  FOR  FREQUENCY  SOLUTION 
NEED  XXXX,  WITH  XXXX  AVAILABLE 

Storage  inadequate  for  FREQ  analysis.  Increase  NCOM. 

(See  Appendix  J.) 

(O)  NOT  ENOUGH  STORAGE  FOR  MODE  SHAPES 
NEED  XXXX,  WITH  XXXX  AVAILABLE 

Storage  inadequate  for  MODE  analysis.  Increase  NCOM. 

(See  Appendix  J.) 
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(N)  NUMBER  OF  SHIPS  =  XXXX 

DYNAMIC  SOLUTION  PRESENTLY  LIMITED  TO  ONE  SHIP 

Frequency  lioinaiii  analysis  requested  with  more  than  one  ship. 

(N)  POSSIBLE  SEQUENCE  ERROR  --  CASE  TERMINATED 

(See  previous  explanation) 

The  dynamic  response  file  is  not  written. 

(O)  REQUESTED  DYNAMIC  EFFECTS  WITH  NO  FRF.QUENCY  DOMAIN  FILE  PROVIDED 
ABORT  ADEQUACY  CHECK 

Dynamic  response  file  was  not  saved  in  the  previous  FRFiO  SAO 
(see  Word  (3)  on  F30L  record). 

(N)  SHIP  MOTION  DATA  EXCEEDS  LIMIT 

NOB  NOH  NOK  NRV 

LIMITS  5  30  30  8 

VALUES  READ  XXX  XXX  XXX  XXX 

The  ship  motion  file  has  arrays  larger  than  the  dimensiotis  in 
SEADYN/DSSM. 

NOB  =  number  of  Froude  Numbers 
NOH  =  number  of  wave  headings 
NOK  =  number  of  wavelengths 
NRV  =  number  of  roll  angles 

(N)  SHIP  NOTION  FILE  ERROR 

WAVELENGTHS  NOT  IN  DECREASING  ORDER 

Check  format  of  ship  motion  file. 

(S)  SHIP  OUTPUT  REQUESTED  WITH  NO  SHIP  [N  THE  SYSTEM 

Random  response  request  for  ship  ignored  when  NSHIPS  <  1. 

(N)  SPECTRUM  ERROR,  NO  FREQUENCIES  FOUND  WITH  SIGNIFICANT  WAVE  HEIGHTS 

Check  spectrum  parameters  and/or  frequency  range. 

(F)  TAPE  I, ABEL  XXXXXX  DOES  NOT  AGREE  WITH  XXXXXX 

The  label  check  failed  on  restart.  The  first  si.x 
characters  on  the  RESTART  title  record  did  not  agree  with 
the  check  word  given. 
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(C)  TOO  MANY  LOAD  SETS  FOR  DRIFT  FORCE  ITERATION 
REFERENCE  STATE  NOT  UPDATED 

j'he  previous  static  analysis  used  three  load  sots.  This 
leaves  no  room  to  store  the  drift  forces. 

(N)  UNRECOGNIZED  ANALYSIS  OPTION  =  XXXX 
TRY  TO  GET  TO  NEXT  CASE 

Usually  indicates  improper  numbers  of  cards  or  cards  out  of 
sequence. 
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10.3  Solution  Option  Execution  Messages 

(N)  DIVERGENCE  ON  STEP  m  AT  XXXXX  (load  factor  or  time) 

Signals  abort  of  MNR  method  after  KNVRT  successive  step 
size  reductions. 

(C)  DIVERGING  NONLINEAR  ITERATION  AT  TIME  XXXX 
STEP  XX  WITH  A  TIME  INCREMENT  OF  XXXX 
LAST  TWO  NORMS  XXXX  XXXX  KOUNT  =  XX 

Signals  a  lack  of  convergence  in  the  Direct  Numerical  Integration 
time  domain  analysis.  The  various  situations  that  lead  to 
this  message  are: 

(a)  Acceleration  more  than  doubled  in  one  iteration  on  the 
component  w’th  the  largest  acceleration. 

12 

(b)  The  displacement  norm  exceeds  1x10 

(c)  The  displacement  norm  increased  in  three  successive 
iterations . 

(d)  The  largest  displacement  increment  exceeds  a  magnitude  of 
IxIO^*^. 

(e)  The  number  of  iterations  exceeds  LMITER  on  SOLD  record. 

This  will  be  followed  by  a  reduction  in  time  step  subject 
to  the  limits  of  KNVRT  on  the  SOLD  record. 

(C)  DIVIDED  STEP  SIZES  XXXX  XXXX 

A  constraint  overshoot  was  detected  in  the  SLI  or  RFB  incremental 
analyses.  The  step  was  divided  into  two  parts  as  indicated  by 
the  message.  The  first  part  represents  the  portion  of  the 
original  step  used  to  get  to  the  constraint.  The  remaining 
portion  of  the  step  was  then  taken  with  the  constraint  imposed. 

A  full  step  size  is  used  after  successful  completion  of  the 
divided  step.  Repeated  divided  steps  with  multiple  constraints 
can  cause  the  solution  to  fail.  In  that  case  small  step  sizes 
should  be  used. 

(N)  DYNAMIC  SOLUTION  DOES  NOT  CONVERGE  AT  TIME  XXXX  WITH  A  TIME  INCREMENT 
OF  XXXX 

LAST  TWO  NORMS  XXXX  XX.XX 

Signals  abort  of  DIM  solution  after  KNVRT  successive  step 
size  reductions. 
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(N,S)  EQUATION  DECOMPOSITION  FAILED  ON  ROW  XXX 


The  simultaneous  equations  in  the  subroutine  SLVBAN  (DEAD, 

LIVF. ,  DYN,  FREQ)  or  subroutine  COMBAN  (FREQ)  are  singular  or 
sufficiently  ill-conditioned  to  appear  singular.  The  row  of 
the  matrix  is  calculated  from  (3*{1-1)+J)  where  I  is  the  node 
number  and  J  is  the  direction  (1,  2  or  3).  Check  node  I  for 
proper  constraint.  This  also  occurs  in  poorly  tensioned 
(soft)  initial  configurations.  Check  for  zero  tensions,  etc. 

It  can  also  occur  from  wildly  divergent  VR  iterations  due  to 
improperly  posed  problem  or  poor  choice  of  parameters 
(usually  damping  is  too  small). 

If  this  occurs  with  the  SLT  or  RFB  methods  or  in  the 
steady-state  response  calculations  in  FREQ  SAG,  the  case  is 
terminated . 

If  it  occur.s  during  a  MNR  analysis,  various  attempts  are  made 
to  remove  the  singularity  and  repeat  the  step.  Failing  in 
these,  the  case  is  terminated. 

(C)  EXCESSIVE  STRAIN  ON  ELEMENT  XX  ON  STEP  YY,  XXXX 

The  VR  iterations  produced  a  strain  larger  than  1,000.0. 

The  strain  value  is  printed  at  the  end  of  the  message.  This 
is  an  indication  of  a  wildly  divergent  iteration  resulting 
from  a  poorly  posed  problem  and/or  a  poor  choice  of  analysis 
parameters.  Damping  is  increased  and  the  iteration  retried. 
Successive  errors  of  this  type  will  lead  to  an  abort. 

(N)  FAILURE  IN  VISCOUS  RELAXATION  SOLUTION 

All  attempts  to  get  convergent  VR  iterations  have  failed.  This 
occurs  when  repeated  norm  increases  or  large  strains  are  encountered. 

(N)  IMPROPERLY  DEFINED  DYNAMIC  PROBLEM 
COMPONENT  XX  HAS  NO  MASS 

A  material  property  input  error  has  been  made  or  zero  length 
element  has  been  included  in  the  model  that  is  not  one  of 
those  reserved  for  payout.  The  component  number  is  computed 
from  3*(N-1)+.J  where  N  is  the  node  number  and  J  is  the  direction 
number  (1,  2  or  3). 

(C)  INCREASING  RESIDUAL  NORM  ON  STEP  XX 

LAST  THREE  RESIDUAL  NORMS  XXXX  XXXX  XXXX 

The  VR  iterations  have  produced  successive  increases  in  the 
force  residual  norm.  Divergence  is  indicated.  Damping  is 
increased  and  the  iteration  restarted. 
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(C)  INCREASING  VELOCITY  NORM  ON  STEP  XX 
NORM  VALUES  XXXX  XXXX  XXXX  XXXX 
REPEAT  STEP 

lnd  iciites  a  strongly  increasing  velocity  beliavlor  in  the  VK 
iterations.  This  is  assumed  to  be  a  signal  of  divergence  if 
it  is  occurring  repeatedly.  Damping  is  increased  and  the  step 
is  repeated.  In  some  situations  this  message  will  occur  many 
times  without  an  abort.  This  is  because  it  gets  intermittent 
good  results  rather  than  successive  increases.  This  usually 
means  the  heuristic  scheme  for  adjusting  the  solution  parameters 
does  not  work  well  for  the  problem.  In  some  situations  this 
can  be  avoided  by  selecting  a  larger  damping  at  the  start. 

(C)  LAST  TWO  RESIDUAL  AND  DISPLACEMENT  NORMS  XXXX  XXXX  XXXX  XXXX 
KOUNT  =  XXXX 

The  four  values  printed  are  the  residual  and  displacement 
norms  for  iteration  (i-1)  and  iteration  i  in  the  MNR  method. 

This  message  signals  an  increasing  norm  indicative  of  diver¬ 
gence  or  a  lack  of  convergence  in  the  number  of  iterations 
given  by  KOUNT,  This  occurs  when  norm  values  are  very  large 
(greater  than  10,000)  or  repeated  increases  are  detected. 

This  will  be  followed  by  a  reduction  of  step  size  subject  to 
the  limits  of  KNVRT  on  the  SOLD  record. 

(C)  NEW  DAMPING  =  XXXX 

Indicates  the  action  taken. 

(C)  NEW  STEP  SIZE  =  XXXX 

Indicates  the  action  taken. 

(N)  ROW  XXX  OF  STIFFNESS  MATRIX  HAS  NEGATIVE  DIAGONAL  TERM  =  XXXX 

Message  follows  previous  message  on  SLI  or  RFB  methods. 

(C)  SINGULAR  EQUATIONS  WITH  NUMERICAL  DAMPING  FACTOR  OF  XXXX 

This  message  is  printed  from  the  MNR  method.  See  Section 
9.0  for  further  explanation.  This  message  is  printed  after 
the  three  trials  described  in  that  note.  This  will  be  followed 
by  a  reduction  in  step  size  subject  to  the  limits  of  KNVKT  on 
the  SOLU  record. 

(C)  SLOW  CONVERGENCE  ON  STEP  XX 

LAST  FOUR  VELOCITY  NORMS  XXXX  XXXX  XXXX  XXXX 
LAST  RESIDUAL  NORMS  XXXX  XXXX  XXXX 

This  is  a  progress  report  on  the  VR  iterations.  Damping 
is  reduced  or  step  size  is  increased  depending  on  the  pattern 
of  residual  norm  changes. 
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(C,N)  STEP  SIZE  REDUCED  TO  XXXX  ON  INCREMENT  XX 
RESIDUAL  AND  DISPLACEMENT  NORMS  XXXX  XXXX 

This  message  follows  the  reduction  of  step  size  in  the  MNR 
method.  It  will  follow  either  of  the  previous  two  messages 
when  the  number  step  size  reductions  is  within  the  limits 
imposed  by  KNVRT  on  the  SOLU  record. 

(C,N)  TIME  STEP  REDUCED  TO  XXXX  ON  STEP  XXXX 
NORN  =  XXXX 

This  message  follows  the  reduction  of  step  size  in  a  DIM 
solution.  It  will  follow  the  previous  message  when  the  number 
of  step  size  reductions  is  within  the  limits  imposed  by  KNVRT 
on  the  SOLU  record. 
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THE  SHIP  MOTION  FILE 


Till'  clatii  for  t  l»e  mot  ion  oqu.iL  ions  for  a  ship  di  ivon  by  harmonic 
voavos  are  providoH  to  iho  SEAOYN  progr.im  througli  the  Ship  Motion 
File.  this  sol  of  (hila  is  assumed  to  be  on  a  soqiuMUial  Oinary  file  on 
logical  unit  08,  i'lris  Appendix  describes  tl>e  format  of  ll\at  tile. 
Notations  from  ReftMonce  16  are  used. 

TIk'  (xpialions  of  motion  for  a  ship  moving  in  waves  on  a  free  .surf  are 
•ire  assumed  to  iinve  the  following  form: 


The  term.s  ot  this  equat  ion  are  assumed  to  bo  provided  on  l  he  shi|>  i»iot  ion 
file  in  a  nondimonsional  form  reflecting  the  effects  of  the  wavelength 
of  the  surface  wave  and  the  relative  heading  between  the  wave  and  the 
ship. 

The  relationships  hetweon  the  terms  of  Equation  «- 1  aud  the 
noiulimonsional  terms  on  the  file  are  given  below: 


M.  u  1  lUMUC.l.K)! 

J  ^ 

(A-2) 

(A-;i) 

Units:  FT^  l'^"' 

Units;  FTl.'^"  ‘ 


(A-  U) 


^  M(r,/L)®  -’  (l/)  iBAAStlj^^)) 


(A-b) 


Hints:  FTl, 


^  M  g  (l/'^''’‘‘^'^I|l)C(J,K)|  (A'l)) 

Units: 

...  .  IBOD(J)  +  iBOn(J  t  3)1 

¥.  -  Mg  M 

lUliV(J)  +  iHKV(.1  +  3)) 

Units:  FL*"  ^ 

where:  M  =  Ship's  mass  (TMAS)  Units:  FT^  1,  ' 

i.  =  Ship's  length  (ELI.)  Units:  1, 

g  ~  gravitational  acceleration  (GKAV)  Units;  LT  ^ 

=  the  Ith  roll  angle  index 

nCj .k)  =  m(j)  >  m(k) 

m(j)  =  0  for  j  S  3 

I  j  ^  3 

i  =  -! 

The  coefficients  are  assumed  to  be  linearized  for  tinii  motion  amplitude 
and  wave  height.  Typical  tinits  for  the  Ship  Motion  File  ar<  ; 

F  -  long  tons  (2,240  pounds) 

1.  -•  feet 
T  -  seconds 


Angles  and  angular  responses  are  assumed  to  he  in  radians. 

In  addition  to  the  ship  motion  coefficients,  the  Ship  Notion  file 
provides  coefficients  (or  estimating  the  steady-slate  approximations  of 
the  second  ordei-,  wave-induced  drift  forces.  These  dimensionless  coef¬ 
ficients  are  used  to  estimate  the  drift  forces  once  tlu'  magnitudes  of 
the  ship  responses,  il,j,  are  obtained.  The  drift  force  loinpoiu'uis  are: 

Siirge  (Units  -  F); 


(A-H) 


Sway  (Units  -  F): 


/  7  n  •  'Vl 

K,,  m(k)  i  I  }  -r-yn  ty(.j,k) 

y  j_j  1^-,  jl-in(k) 

Yaw  (Units  '  FL) : 


(A- 9) 


^  ;pixj)  -raio 


7  7  n 

+  n  R  11  - fTT  TP(J) 

1  =  1  k=l  1. 


whwre:  =  wavo  amplitude  (real  number) 


(A-10) 


augnmented  stiip  response 


m(j)  =  0  for  ,i  =  1 ,2 ,3,7 

1  for  ,i  -  4,5,6 

(  )*  means  complex  conjugate 

i  wt 

This  presumes  the  wave  motion  is  given  by  c  ,  where  il^  is  real  and 
the  ship  response  is  i).  e"^  . 

The  coordinate  system  presumed  for  ship’s  motions  and  forces  is  a 
righthand  cartesian  system  with  its  origin  at  the  ship's  center  of 
gravity,  its  X-axis  positive  aft,  its  Y-axis  positivt*  stairhoard,  and  its 
Z-axis  positive  upward.  Tae  angular  convention  for  the  relative  heading 
between  the  ship  and  the  waves  assumes  the  following; 


Wave  Hejrding 


Description 


0®  Following  seas 

90®  Beam  seas  with  waves  traveling  from 

port  to  starboard 

180®  Head  '■eas 

270®  Beam  seas  with  waves  traveling 

from  starboard  to  poi t 

The  differences  between  the  wave  heading  convent  ion  and  the  ship's 
coot di nates  shoulil  be  noted. 


Tho  Ship  Motion  Kilo  is  organized  in  logioai  records.  Tlie  specific 
contents  of  each  record  will  be  described  below.  I'liere  are  seven  distinct 
record  types.  The  first  two  records  contain  data  which  are  independent 
of  wave  heading  or  wavelength.  Record  types  3  through  7  are  dependent 
on  heading  and  wavelength  and  are  repeated  in  a  nested-loop  fashion. 

The  overall  form  is; 


Record  1 
Record  2 


Record  3 

Record  4  (presented  tnly  on  first  loop 

,,  ...  per  wavelength  in  the  NCEL  file) 

Record  5 

Record  6 
Record  7 


The  wave  headings  are  assumed  to  be  listed  in  decreasing  order  with 
+180**  being  the  largest  allowed.  The  interpolation  routines  assume  the 
values!  given  for  +180°  will  be  used  for  -180°,  therefore,  data  for  -180° 
need  not  be  given . 

The  wavelengths  are  assumed  to  be  listed  in  decreasing  order  (i.e., 
increasing  frequency  order). 

The  individual  records  of  the  'lie  are  de.scribed  in  terms  of  the 
FORTRAN  read/write  lists  associated  with  each  record. 

RIX0W)_1  NAME  1,  NAME  2,  NAME  3 

Three  Hollerith  variables  providing  identifying  data. 

RECORD  2  (TITO(I),  1=1,  12),  WORD,  WORD  2,  WORD  3,  ELL,  BEAM,  DRAFT,  TVOL, 
TMAS,  TPST,  ZG,  CBV,  NOB,  (FN(I),  1=1,  NOB),  NOH,  (HDGl(I), 

1=1,  NOH).  NOK,  (BAM(I),  1=1,  NOK) ,  VNY,  GRAV,  NRV,  (RANG(I), 

1=1,  NRV),  ((GMU(I,J),  J=l,6),  1=1,6),  ((DC(I,J),  J=l,6),  1=1,6) 

TITO  =  Hollerith  title  consisting  of  12  six-character  words 
WORD  =  length  unit  label  (six-character  Hollerith) 

WORD  2  =  force  unit  label  (six-character  Hollerith) 

WORD  3  =  moment  unit  label  (six-character  Hollerith) 

ELL  =  ship's  length  (L) 

BEAM  =  beam  (L) 

DRAFT  =  draft  (L) 

TVOL  =  ship's  volume  is  obtained  from  (ELL/2) ^  •  TVOL 
TMAS  =  ship's  mass  (FT^L  *) 

TPST  =  longitudinal  distance  from  c.g.  to  forward  most 
station  is  obtained  from  (ELL/2)  *  TPST 
ZG  =  vertical  distance  from  water  line  to  c.g.,  (+  up)  (L) 
CBV  =  vertical  distance  from  water  line  to  center  of 
buoyancy  (+  up)  is  obtained  from  ELL  •  CBV 
NOB  =  number  of  ship  speeds  (SEADYN  expect.s  only  one) 


on 

Wave 

Length 


Loop 

on 

Wave 

Heading 
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F  N  ( ] ) 


=  the  Fronde  ixinnber.s  for  eacli  speed  (only  one 
expec'Led) 

NOH  ■■  liiiniher  of  wave  headings 

lll)(;i(I)  =■  l.he  wave  headings  listed  in  decreasing  order 
starting  with  180“  and  proceeding  no  further 
than  -179“ 

NOK  -  number  of  wavelengths 

BAM(I)  =  uondimcnsional  wavelength  in  decreasing  order, 

A  =  F.LL  .  BAM(T) 

VNY  fluid  viscosity  (L^T)  _2 

GRAV  =  gravitational  acceleration  (LT  ) 

NRV  =  number  of  roll  angles 

RANG(l)  =  the  values  of  roll  angles  (radians)  listed  in 
increasing  order 

GMIJ(I,J)  -■  the  nondimensional  mass  matrix 

DC(I,J)  =  (lie  non-dimensional  hydrostatic  restoring  matrix 

RECORD  3  MM,HDGI(MM),  JJ,  FN(JJ),  U,  BAM(LL) 

MM  -•  heading  number 
JJ  ==  speed  number 

RECORD  4  ((DA(I,J),  J=l,6),  1=1,6),  ({DB(1,J),  J=l,6),  1-1,6) 

DA(I,J)  =  the  nondimensional  added  mass  matrix  for  that 
wavelength 

DB(1,J)  =  the  nondimensional  wave  damping  matrix  for  that 

wavelength 

RECORD  5  (B0D(1\  B0D(I+3),  BEV(l),  BEV(I+3),  1=1,3) 

BOD,  BEV  =  the  nondimensional  wave  force  coefficients 

RECORD  6  (B44S(I),  1=1,  NRV) 

BA4S(I)  =  the  nonlinear  roll  damping  terms  that  are 
added  to  the  linearized  damping  matrix 
depending  on  the  size  of  the  roll  angle 
(nondimens ional ) 

RECORD  7  ((TX(1,J),  J=l,7),  1=1,7),  ((TY(I,J),  J=l,7),  1=1,7),  ((TM(I,J), 

J=l,7),  1=1,7),  (TP(I),  1=1,7) 

The  nondimensional  drift  force  coefficients. 
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Appendix  B 
SHIP  LOAD  DATA  FILE 


The  purpose  of  the  ship  load  data  file  is  to  provide  static  loads 
for  ships  and  other  rigid  bodies  for  wind  and  surface  currents  at  various 
headings.  The  file  can  be  constructed  and  saved  as  a  library  of  static 
load  functions;  Appendix  F  describes  the  procedure  used  in  DSSM  to  scale 
loads  between  similar  ships  using  this  library  of  files.  Input  is 
provided  in  a  FORTRAN  rigid  format.  Each  ship  load  data  set  consists  of: 

SHIP  LOAD  TITLE  RECORD 
UNIT  LABEL  RECORD 
SHIP  PARAMETERS 
WIND  RECORD 
WIND  HEADING  RECORD (S) 

WIND  FORCE  RECORD(S) 

SURFACE  CURRENT  RECORD 

SURFACE  CURRENT  HEADING  RECORD(S) 

SURFACE  CURRENT  FORCE  RECORD(S) 

Each  of  the  data  records  is  described  below. 

The  free*fie]d  input  routine  will  automatically  read  these  rigid 
format  records  when  they  are  placed  between  the  (  )  delimiter  records 

(see  Section  5.0).  These  data  can  appear  anywhere  in  the  input  stream 
after  the  first  title  record  set.  V^en  SEADYN  encounters  NSFILE>0, 
these  rigid  format  data  are  processed  into  the  ship  load  data  file. 

This  is  to  be  done  only  once  in  any  job  deck.  It  is  not  possible  to 
stack  new  problem  cases  using  NEW  which  intend  to  define  two  distinct 
ship  load  data  files. 
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SHIP  LOAD  TITLE  RECORD  (12A6)* 

Variable 

Columns  Name  Description 


1-72 

SHPCAP 

Any  descriptive  title 

UNIT  LABEL 

RECORD  (A6,  4X, 

A6,  4X,  A6,  4X,  A6) 

Columns 

Variable 

Name 

Description 

1-6 

WLBL 

Wind  force  label  (e.g.,  "TONS,"  "POUNDS  ) 

11-16 

CLBL 

Current  force  label 

21-26 

LLBI. 

Length  label  (e.g.,  "FEET,"  "METERS") 

31-36 

VLBL 

Velocity  label  (e.g.,  "KNOT,"  "FT/SEC  ) 

NOTE:  These  labels  are  output  with  the  ship  data  as  a  reminder  of  the 
units  used.  They  are  used  for  no  other  purpose. 

SHIP  PARAMETERS  (8E10.0) 

Columns 

Variable 

Name 

Description 

1-10 

TSLT 

Total  ship  length  (L) 

11-20 

TSAE 

End  projected  wind  area  (L*) 

21-30 

TSAS 

Side  projected  wind  area  (L^) 

31-40 

TSWL 

Water  line  length  (L) 

41-50 

TSB 

Beam  at  midships  (L) 

51-60 

TSD 

Draft  at  midships  (L) 

61-70 

TSDSP 

Volume  displacement  (T.^) 

71-80 

TSAP 

Propeller  projected  area  (L^) 

*FORTRAN  format  specification. 


WIND  RECORD  (215,  6E10.0) 


Variable 


Columns 

Name 

Description 

1-5 

NWIND 

Number  of  wind  velocity  tables  (max  is  5) 

6-10 

NTHETW 

Number  of  headings  in  each  table  (max  is  20) 

11-20 

SCALE 

Test  scale  factT  (A  means  1/A^^  scale) 

21-30 

WNDVEL(l) 

First  wind  velocity  (smallest)  (I.T 

61-70 

WNDVEL(5) 

Fifth  wind  velocity 

WIND  HEADING  RECORD(S)  (8E10.0) 


(Repeat  as  required  to  get  NTHETW  entries) 


Columns 

Variable 

Name 

Description 

1-10 

WNDHED(l) 

First  wind  heading  (degrees) 

11-20 

etc. 

NOTES 

1.  Headings  should  be  between  0°  and  360**,  listed  from  the  smallest  to 
the  largest. 

2.  If  the  largest  value  is  180°,  the  loading  functions  are  assumed  to  be 
symmetric  about  180°  for  the  end  forces  and  skew-symnetric  for  the  side 
forces  and  yaw  moments. 

3.  The  angle  is  measured  relative  to  the  ship's  local  coordinate  system, 
which  is  illustrated  below: 


WIND  FORCE  RECORD (S)  (3E10.0) 


Columns 

Variable 

Name 

Description 

1-10 

WNDCOEd ,  1  d) 

End  force  for  I*^^  heading  and 

velocity  (F) 

11-20 

WNDC0E(I,2,J) 

Side  force  for  I*^  heading  and  J 

velocity  (F) 

21-30 

WNDC0E(I,3,J) 

Moment  for  heading  and  velocity  (FL) 

(I  varies  before  J) 

SURFACE 

CURRENT  RECORD(S) 

(215,  6E10.0) 

1-5 

NCRNT 

Number  of  current  tables 

6-10 

NTHETC 

Number  of  headings  in  each  current 

table 

11-20 

TDEPTH 

Test  water  depth  (L) 

21-30 

CURVEL(l) 

First  current  velocity  (smallest) 

(lt"S 

61-70 

CURVED (5) 

Fifth  current  velocity 

SURFACE 

CURRENT  HEADING  RECORD(S)  (8E10.0) 

(Repeat  as  required 

to  get  NTHETC  entries) 

Columns 

Variable 

Name 

Description 

1-10 

CURHED(l) 

First  Current  Heading  (Degrees) 

11-20 

etc . 

(See  notes  for  WIND  HEADING  RECORD) 

SURFACE  CURRENT  FORCE  RECORD (S)  (3E10.0) 

(One  record  for  each  heading  repeated  for  each  velocity) 

1-10  CURCOEd  ,1  ,J)  End  force  for  heading  and  velocity  (F) 

11-20  CURC0E(I ,2,J)  Side  force  for  I^^  heading  and  velocity  (F) 

21-30  CURC0E(I ,3, J)  Moment  for  I^^  heading  and  velocity  (FL) 

(I  varies  before  J) 
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The  ship  load  file  contains  one  logical  record  for  each  ship  cata¬ 
logued  on  the  file,  and  it  is  written  in  a  binary  form  with  the  following 
FORTRAN  statement; 

WR1TE(10)NW1ND ,NTHETW,WNDVEL,WNDHED .VWDCOE .SCALE ,NCRNT .NTHETC , 

CURVE! , CURHED , CURCOE , TDEPTH , TBLOCK , TSLT , TSAE , TSAS , TSWL , TSB , TSD , 

TSDSP , TS AP , SHPC AP , WLBL , C  LBL , LLBL , VLBL 

A  number  of  the  items  in  the  list  are  arrays,  and  they  are  written  in 
their  entirety  using  the  implied  DO-LOOP  feature  of  FORTRAN  I-O  statements. 

A  description  of  each  item,  including  the  dimensions  of  the  arrays,  is  given 
below: 


VARIABLE 


DESCRIPTION 


NWIND 

NTHETW 

WNDVEL(5) 

WNDHED(20) 

WNDCOE(20,3.5) 


SCALE 

NCRNT 

NTHETC 

CURVEL(5) 

CURHED  (20) 

CURCOE  (20, 3, 5) 


TDEPTH 

TBLOCK 


Number  of  wind  velocity  tables 
Number  of  headings  in  each  wind  table 
Array  of  wind  velocities 
Array  of  wind  headings 

Array  of  wind  load  coefficients  giving  values  for  up 

to  20  headings  for  end  force,  side  force,  and  yaw 

moment  for  up  to  five  wind  velocities 

Scale  for  wind  load  tests  (A  means  l/A^*'  scale) 

Number  of  current  velocity  tables 

Number  of  headings  in  each  current  table 

Array  of  current  velocities 

Array  of  current  headings 

Array  of  current  load  coefficients  giving  values  for 
up  to  20  headings  for  end  force,  side  force,  and  yaw 
moment  for  up  to  five  current  velocities 
Water  depth  for  test 
Ship's  block  coefficient 
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TSLT 

Total  ship  length 

TSAE 

End  projected  wind  area 

TSAS 

Side  projected  wind  area 

TSWL 

Water  line  length 

TSB 

Beam  at  midships 

TSD 

Draft  at  midships 

TSDSP 

Volume  displacement 

TSAP 

Propeller  projected  area 

SHPCAP(12) 

Title  of  12  six-character  Hollerith  words 

WLBL 

Wind  force  label,  six-character  Hollerith  word 

CLBL 

Current  force  label,  six-character  Hollerith  word 

LLBL 

Length  label,  six-character  Hollerith  word 

VLBL 

Velocity  label,  six-character  Hollerith  word 
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Appendix  C 

BUILT-IN  DRAG  COEFFICIENTS 


In  the  event  that  fluid  loading  is  required  and  the 
are  not  given  by  the  subroutine  DRAGCO,  the  program  will 
drag  coefficients  from  a  set  of  built-in  functions.  The 
are  used  whenever  INDRAG  =  0  (PROB  record,  Section  6.3), 
the  drag  coefficient  codes  on  the  cable  materials  and/or 
bodies.  If  INDRAG  ^  0,  the  default  coefficients  will  be 
components  that  have  zero  drag  coefficient  codes. 


drag  coefficients 
select  default 
default  coefficients 
regardless  of 
the  lumped 
used  for  those 


The  default  coefficients  are: 

Spherical  Bodies 

Reynolds  number,  R^  = 

C„  =  0  for  R  g  0.1 

D  e 

C„  =  0.004  =  13.46/(R  )''•  '  for  0.1  <  R  S  1000 

D  e  e 

C_  =  0.47  for  1000.  <  R  ^10^ 

D  e 

C_  =  0.12  for  R  >  10^ 

D  e 


V  d  _  velocity  x  body  diameter 
U  kinematic  viscosity 


(C-1) 


(C-2) 


Cylindrical  Bodies  and  Cabl e  Elements 


R  = 
e 


''n 


normal  velocity  x  body  diameter 
kinematic  viscosity 


eT 


V  d 

T  _  tangential  velocity  x  body  diameter 
o  ~  kinematic  viscosity 


(C-3) 


(C-4) 


C-1 


(C-5) 


s 

0 

for 

R  ^0.1 
e 

s 

£ 

0.45  +  5.93/(R 

e 

for 

0.1  <  K  %  400 

6 

s 

= 

1.27 

for 

400  ^  R  10^ 
e 

c„ 

0.3 

for 

R  >  10^ 

N 

e 

s 

= 

0 

for 

0.1 

s 

= 

for 

0.1  <  R^.J.  ^  100 

s 

= 

0.062 

for 

R^.j,  >  100.55 
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Appendix  D 

MOORING  COMPONENT  INVENTORY 


The  mooring  component  inventory  contains  data  tables  for  the  following: 
Anchors : 

Navy  standard  stockless 
NAVSHIPS  lightweight 
NAVFAC  STATO 

Biwys : 

Bar  riser  chain  type 
Chain; 

Steel  stiid-link 
Hawsers : 

Samson  braids  --  2-in-l^  Nylon 

2-in-lg  Power  Braid 
2-in-l®  Stable  Braid 
12  Strand  Blue  Streak** 

The  inventory  lists  weights,  buoyancies,  and  strengths  in  pounds. 

Lengths  and  buoy  dimensions  are  in  feet.  Hawser  and  chain  sizes  are  in 
inches.  These  units  may  be  converted  to  those  needed  in  the  analysis  by 
providing  the  appropriate  conversion  factors  on  the  INVE  data  record 
(see  St;ction  6.9).  The  contents  of  the  inventory  can  be  obtained  by 
setting  the  print  flag  in  Word  (A)  of  the  INVE  record.  The  listing  of 
the  present  inventory  is  presented  below. 
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c  0  M  F  0  N  !■:  N  T 


inventory 


ANCHOR  TYPE  =  NAVY  STD  STOCKLESS 


WEIGHT 

FED.  STOCK  NO. 

.hold.  POWER 

. JOOOOE+04 

(FIRM  SAND) 

C2040-516-7758 

.21000E+05 

.500UUE+04 

C2040-516-7757 

.35000E+05 

.60000E+04 

C2040-516-7756 

.42000E+05 

. 70000E+04 

C2040-516-7755 

.49000E+05 

.90000E+04 

C2040-516-7754 

.63000E+05 

. lOOOOE+05 

C2040-272-2244 

.70000E+05 

. 13000E+05 

C2040-272-2245 

-91000F+05 

. 14500K+05 

C204 0-2 7 2-2246 

. 10I50E+06 

.  180001. +05 

C2040-516-7753 

.  1260<''>0tj 

.20000E+05 

C2040-272-2247 

. 14000E+U6 

.25000E+05 

C2040-272-2242 

.  17500E+06 

.30000E+05 

C2040-272-2243 

.21000E+06 

.40000E+05 

C2040-277-2423 

.28000E+06 

ANCHOR  TYPE 

=  NAVSHIP(LWT) 

WEIGHT 

FED.  STOCK  NO. 

HOLD.  POWER 
(FIRM  SAND) 

. lOOOOE+03 

H2040-377-8600 

.28374E+04 

. 15000E+03 

H2040-377-8601 

.39565E+04 

.20000E+03 

H2040-377-8602 

.50091E+04 

.30000E+03 

H2040-377-8603 

.69848E+04 

.50000E+03 

H2040-377-8604 

.  10619Et-05 

.75000E+03 

H2040-377-8605 

. 14807E+05 

. lOOOOE+04 

H2040-377-8606 

.I8746E+05 

.  15000E+04 

H2040-377-8607 

.26140E+05 

.20000E+04 

H2040-377-8608 

.33095E+05 

.25000E+04 

H2040-377-8609 

.39740E+05 

.30000E+04 

H2040-377-8610 

.46I48E+05 

.40000E+04 

H2040-377-8611 

.58426E+05 

.50000E+04 

H2040-378-5633 

.70157E+05 

.60000E+04 

H2040-378-5634 

.8)470E+05 

.  lOOOOE+05 

H2040-377-8612 

.12385E+06 

.  13000E+05 

H2040-377-8613 

.  15358E+06 

ANCHOR  TYPE  = 

NAVFAC  STATO 

WEIGHT 

FED.  STOCK  NO. 

HOLD.  POWER 
(FIRM  SAND) 

.20000E+03 

2CF2040-800-9659 

.40000E+04 

.30000E+04 

2CF2040-702-7864 

.60000E+05 

.60000E+04 

2CF2040-702-6785 

.  12000E+06 

.90000E+04 

2CF2040-702-6786 

.18000E+06 

. 12000E+05 

2CF2040-702-6787 

.24000E+06 

. 15000E+05 

2CF2040-801-7938 

.30000E+06 
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STEKL  STUD-LINK  CHAIN 


SIZE  STRENGTH  WEIGHT/LENGTH 


.7500 

.8750 

1.0000 

1.1250 

1.2500 

1.3750 

1.5000 

1.6250 

1.7500 

1.8750 
2.0000 

2.1250 

2.2500 

2.3750 

2.5000 

2.6250 

2.7500 

2.8750 
3.0000 

3.1250 

3.2500 

3.3750 

3.5000 

3.6250 

3.7500 

3.8750 
A. 0000 
^.1250 

4.2500 

4.3750 

4.5000 


48550. 

65280= 

84500. 

106080. 

130070. 

156330. 

185060. 

216030. 

249210. 

284540. 

322000. 

361530. 

403100. 

446660. 

492190. 

539620. 

588930. 

640070. 

693000. 

747680. 

804070. 

862130. 

921810. 

983080. 

1045900. 

1110210. 

1176000. 

1234200. 

.  1311790. 
1381330. 
1452930. 


5.5556 

7.7778 

9.4444 

12.2222 

15.0000 

17.7778 

21.1111 

24.4444 

28.3333 

32.7778 

36.6667 

41.1111 

46.6667 

51.6667 

57.7778 

63.3333 
70.0000 

76.6667 

83.3333 

91.1111 

98.3333 

106.1111 

114.4444 

122.7778 

131.1111 
140.0000 

148.8889 

158.8889 
170.0000 

183.8889 

198.3333 
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BUOY  DATA 


BUOY  TYPE  =  BAR  RISER  CHAIN 


0.  u. 

HEIGHT 

WEIGHT 

NOM.  BUOYANCY 

MAX.  BUOYANCY 

FEU.  STOCK  NO. 

6.53125 

4.03125 

2200. 

3562. 

4300. 

C2050-223-365? 

7.03125 

5.03125 

2500. 

5835. 

7518. 

C2050-264-4497 

9.50000 

5.00000 

7700. 

7420. 

10445. 

C2050-223-3665 

10.50000 

6.50000 

9600. 

14414. 

20879. 

C2050-223-3662 

10.50000 

7.50000 

JOIOO. 

17608. 

25921. 

C2050-264-44981 

D-A 


2-IN-l  NYLON  2-IN-i  ?OWES  BRAID  2-IN-l  Sl'ABL  SilAID  12  ST.  .BLL'ESTRLAK 
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Appendix  E 

RESTART  FILE  STRUCTURE 


The  SEADYN  program  creates  up  to  three  restart  files  (one  each  for 
the  DEAD,  LIVE,  and  UYN  subanalyses).  Multiple  selections  of  a  subanalysir, 
type  simply  extends  the  file  unless  a  rewind  is  signaled  on  the  SAVE 
data  record.  A  counter  is  provided  for  each  of  the  files  to  keep  track 
of  how  many  restart  records  have  been  written.  The  FORTRAN  file  codes 
used  are: 


01  -  DEAD 

02  -  LIVE  (and  TSSS) 
03  -  DYN 


Each  time  the  file  is  rewound,  the  counter  for  that  file  is  set  to  zero, 
and  a  label  record  is  written.  The  write  statement  is: 


WRITE  (NFILE)  (TITLE  (1),  I  =  1,  NHED) ,  NINA,  NPRECZ 


where:  TITLE 

NINA 

NPRECZ 


NHED 


page  heading  title  card  for  the  run  (Each  word  is 
assumed  to  have  10  characters.) 
size  of  unlabeled  common  when  the  file  is  saved 
precision  number  for  floating  point  numbers 

1  =  single  precision 

2  =  double  precision 

the  number  of  words  in  the  title  =  8 


Each  restart  save  operation  uses  the  following  write  statement: 


WRITE(NTAPE)  (A(I) ,1=1 , NINA) , (B(l) ,1=1 , NINE) , (C(T) ,1=1 ,NINC) , 
+(RL(I),1=1 ,NINRL),(P(I),1=1,N1NP0),(T(1),I=1,N1NT),(SH(I),I=1,NINS 
+HP) , (DP(I) ,I=1,N1NDSP),(STM(I),I=1,NINSTM)  , NFILE 
.  ,(IAABU(I),I+1 ,NINB1),(IAACA(1),1=1,NINCI) 

.  ,(IAACL(I) ,1=1  ,NINCLI) ,(1AADS(I),I=1,NINDSI) 

.  ,  (I.AAPOd  )  ,1+1  .NINPOI)  ,  (lAASKP(I)  ,1=1  .NINSHI  ) 

.  ,(IAASTM(I),1=1,NINSTI),(IAATIM(I),I=1,NINTMI) 

+ ,  DLD ,  WLD ,  DYN ,  CHECKR ,  NOVEL  ,NOITER ,  NOFLUD ,  NOLOAD ,  FEF.DBK ,  POUT , 

+  REFUP.STEPUP 
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where  the  arrays  are  defined  by: 


COMMON/ ACOM/  A(l) 

COMMON/BUOYS/  B(l) 

COMMON/CABLE/  C(l) 

COMMON/CONTRL/  RL(1) 

COMMON/DSPCON/  DP(1) 

COMMON/PAYOUT/  P(l) 

COMMON/TIMED/  T(l) 

COMMON/SHIPS/  SH(1) 

COMMON/ STRUM/  STM(l) 

COMMON/ IBUOYS/  lAABU(l) 

COMMON/ I CABLE/  lAACA(l) 

COMMON/ ICNTRL/  lAACL(l) 

COMMON/ IDSPCN/  lAADS(l) 

COMMON/ 1 PAYOT/  lAAPO(i) 

COMMON/ ISHIPS/  lAASHP(l) 

COMMON/ ISTRUM/  lAASTM(l) 

COMMON/ 1  TIMED/  lAATIH(l) 

COMMON/LOGIC/DLD.WLD ,DYN , CHECKR, NOVEL, NOITER .NOFLUD , NOLOAD, 

1  FEEDBK,POUT,REFUP,STEPUP 

The  sizes  of  the  arrays  are  given  by: 

COMMON/SIZE/  NINA,KINB,NINRL,NINDSP,NINP0,NINSHP,N1NSTM,NINT, 

1  N I NC , NCOM , I F I LE ( 4 ) , NPRECZ , N1 NBI , N I NC 1 , N INC  L: , NI NDS I , 

2  NINP01,NINSH1,N1NSTI,N1NTM1 

These  sizes  are  identified  in  Reference  2,  and  count  the  number  of 
single-precision  words  to  be  read/written.  This  count  is  adjusted  for 
double-  nr  single-precision  conditions.  The  arrays  are  always  to  be 
treated  as  single  precision  in  the  RESTART  routine  even  thouglj  they  can 
contain  inixed-iiouble  precision  and  fixed-point  data  in  the  rest  of  the 
program.  The  actual  contents  of  the  labeled  common  blocks  are  defined 
in  the  calling  program  with  the  appropriate  word  format. 
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Appendix  F 

SHIP'S  LOADING  FUNCTIONS 


This  Appendix  describes  the  tabular  approach  used  in  SEADYN  to 
obtain  static  loads  on  a  ship  subjected  to  winds  and  surface  currents. 

The  ship  loads  are  assiamed  to  be  applied  through  the  center  oi 
gravity  of  the  ship.  Three  sources  of  loading  are  considered:  wind 
loads,  surface  current  loads,  and  point  loads  representative  of  working 
loads.  The  point  loads  are  specified  through  the  normal  loading  options 
(see  Section  7.1.7,  LOAD  record).  It  is  assumed  that  the  point  loads  do 
not  change  their  magnitude  or  global  directions  as  the  ship  moves  to  a 
new  position.  Specification  of  the  wind  and  current  loadings  is  somewhat 
more  complicated.  The  SEADYN  program  provides  two  approaches  to  defining 
these  loads.  The  first  approach  is  in  the  form  of  loading  tables  that 
give  loads  versus  ship's  heading  relative  to  the  flow.  The  second 
approach  uses  approximate  analytical  expressions,  which  are  described  in 
Appendix  G. 

The  tabular  approach  is  based  on  the  procedures  given  in  NAVFAC's 
Design  i'aniial.  26  (Ref  17).  The  DM-26  approach  utilizes  experimental 
measurements  for  the  forces  and  moments  for  various  headings  of  wind  and 
current,  for  a  set  of  "representative"  vessels.  Similarity  scaling  is 
then  applied  to  get  loading  values  for  ships  other  than  the  test  model'j. 

The  DH'26  procedure  begins  with  a  set  of  load  measurements  obtained 
from  snbscale  tests  on  a  representative  ship's  model  or  any  other  available 
source.  The  measurements  give  values  for  the  lateral  and  longitudinal 
forces  and  yaw  moment  versus  flow  heading  and  flow  velocity.  These 
measurements  represent  the  combined  effects  of  such  pheiiomena  as  profili: 
and  friction  drag,  lift-induced  side  forces,  and  shifts  in  the  center  of 
pressure.  Tables  of  these  measurements  can  be  specified  as  either  input 
to  the  program  or  as  a  special  ship  loading  file  previously  generated 
and  saved  for  subsequent  referencing  by  SEADYN.  This  saved  file  is 
assumed  to  be  on  logical  unit  03  (see  Appendix  J). 

Given  the  headings  of  wind  and  surface  current  relative  to  the 
ship,  the  leads  are  obtained  by  linear  interpolation  in  the  tables.  In 
the  event  that  there  are  tables  provided  for  more  than  one  velocity,  the 
table  for  the  velocity  nearest  the  one  specified  in  the  analysis  will  be 
used.  This  is  determined  by  comparing  the  squares  of  the  velocity 
ratios . 

After  the  load  coefficients  are  obtained  from  the  tables  for  the 
given  heading,  they  must  be  scaled  to  account  for  differences  in  the 
conditions  modeled  in  the  test  and  for  those  being  analyzed.  The  scaling 
accounts  for  differences  in  flow  velocity,  water  depth,  and  ship  geometry. 
The  formulas  for  adjusting  for  these  effects  are  given  in  DM-26  and  are 
restated  here  for  completeness. 
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C.  V'"  F  (A  /A^  ) 
f  ms  s  ts 


C,  V*"  F  (A  /A^  ) 
t  me  e  te 


C  M  (A  /A_)  (L/L.  ) 
m  m  s  ts  t 


(F-1) 

(F-2) 

(F-3) 


=  lateral  force  on  ship 
=  longitudinal  force  on  ship 
-  yawing  moment  on  ship 
=  lateral  force  on  model 
=  longitudinal  force  on  model 
=  yawing  moment  on  model 
=  wind  velocity 

=  side-projected  area  above  the  water  line  of  ship  being 
analyzed 

=  side-projected  area  above  the  water  line  of  modeled  ship 

=  end-projected  area  above  the  water  line  of  ship  being 
analyzed 

=  end-projected  area  above  the  water  line  of  modeled  ship 
=  length  of  ship  being  analyzed 
=  length  of  modeled  ship 

(F-4) 

(F-h) 

=  linear  scale  of  the  model  (e.g.,  50  to  1;  S  =  50) 

=  wind  velocity  used  in  model  test 
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CURRENT 


^2  “  ''l 

'l  "  ''2 
’'2  =  ■'l  *2'^ 
”2  “  ”1 


where:  h 
V 

hi 


depth  of  water 

velocity  of  current 

water  line  length  of  vessel 


F  =  lateral  or  longitudinal  resisting  force 

A  =  displacement 

M  =  yaw  resisting  moment 


(F-6) 

{F-7) 

(F-8) 

(F-9) 


Subscript  1  denotes  the  full-scale  vessel  for  which  the  model  test  was 
made,  and  subscript  2  denotes  the  vessel  being  analyzed. 

When  the  velocity  from  Equation  F-7  does  not  correspond  to  one  of 
the  tables  given  for  the  model  test,  then  the  forces  and  moments  must  be 
selected  from  the  tables  corresponding  to  the  velocity  nearest  the  value 
of  V,  in  Equation  F-7.  It  will  then  be  necessary  to  adjust  the  values 
by  the  square  of  the  ratio  of  the  V.  velocity  and  the  velocity  represented 
in  the  tables, 

It  is  quite  likely  that  the  depth  at  the  proposed  mooring  site  will 
not  be  the  same  as  that  obtained  for  h„  in  Equation  F-6.  In  that  event, 
a  correction  for  depth  is  required.  Dn-26  suggests  that  the  correction 
be  made  assuming  an  inverse  relationship  with  the  side  resistances  at 
the  two  depths  in  question.  The  curves  given  in  Graph  124  (EC-2)  of 
DM-26  are  used  along  with  Equation  F-6  for  this  purpose.  The  data  are 
given  in  tabular  form  and  the  side  resistances  are  obtained  by  logarithmic 
interpolation.  The  resistance  for  a  depth  greater  than  that  in  the 
table  will  be  the  last  value  in  the  table. 

The  adjustments  for  current  velocity  and  depth  are  summarized  by 
the  following  equations; 


F' 

s2 


s2 


(F-10) 


F-3 


(F-11) 


F' 

e2 


1 


h  ^  2 

tl  J 


F  -  C  A 
I  el  2  ^  p  tl 


*TP  S  *  ''2^ 


(F-12) 


where: 


the  depth  scaling  factor 

the  velocity  at  which  the  test  data  was  obtained 
the  propeller  projected  area 
the  propeller  drag  coefficient 
fluid  density 


The  primes  indicate  the  value  has  been  adjusted  to  the  desired  condi¬ 
tions  for  the  mooring  site.  Equation  F-11  reflects  the  adjustment  in 
the  longitudinal  force  recommended  by  DM-26  with  the  assumption  that 
(1/2)  P  C  “  2.88  (with  V  in  knots).  Assuming  the  specific  weight  of 
seawater  is  64  Ib/ft^  and  the  acceleration  due  to  gravity  is  32.2  ft/sec^, 
then  C  =  1.00.  The  form  using  (1/2)  p  C  rather  than  2.88  is  required 
to  inakS  the  procedure  dimensionally  indepindent. 
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Appendix  G 

BUILT-IN  LOAD  FUNCTIONS  FOR  SHIPS 

This  Appendix  describes  the  analytical  approach  for  obtaining 
static  loads  on  a  ship  subjected  to  winds  and  surface  currents.  The 
approximate  analytical  expressions  for  ship's  loading  are  based  pri¬ 
marily  on  the  work  of  Hughes  (Ref  18),  standard  Naval  architectural 
formulas  (Ref  19),  and  Altman  (Ref  20). 

The  wind  loading  is  given  by: 

F  =  K  p  V^(A  sin^  9  +  A  cos^  8)  cos  (a-6)  (G-l) 

d  S  £ 

where;  K  =  constant,  0.6 

F  =  resultant  wind  force 
=  mass  density  of  air 
V  =  wind  velocity 

0  =  wind  heading  relative  to  the  bow 

a  =  heading  of  the  resultant  wind  force  relative  to  the  bow 

A  =  side  projected  area  of  ship  above  water  line 

s 

A^  =  end  projected  area  of  ship  above  water  line 

The  heading  of  ^|ie  resultant  wind  force,  a,  is  approximated  as  a 
function  of  0  in  a  7  order  polynomial  as  follows: 

a  =  0.0715608  +  7.95A38]  0  -  0.3254561  0^ 

+  0.0073131  0^  -  9.3966  x  lo'^  0^ 

+  6.85008  X  lO’^  0^  -  2.6323  x  lo'®  0^ 

+  4.1453  x  lO’^^  0'  (G-2) 

In  Equation  G-2  both  0  and  a  are  measured  in  degrees. 


G-l 


The  distance  between  the  ship  forward  perpendicular  and  the  center 
of  wind  pressure,  X  ,  can  be  approximated  as  a  polynomial  function  of 
the  wind  direction, ^6.  This  relationship  is: 

X  .5 

=  0.2004112  +  0.0048641  6  -  4.52442  x  10‘^  0“^ 

+  5.45736  X  lO"^  0^  -  3.78789  x  lo"^  0^ 

+  1.02881  X  10“^^  0^  (G-3) 


Here,  as  above,  G  is  measured  in  degrees.  Also,  L  is  the  overall  ship 
length.  The  yawing  moment  clue  to  wind  is  then  approximated  by: 


M 


W 


F  L  sin  0 


(tt-4) 


Analytical  expressions  for  the  resistances  from  cuirent  effects 
utilize  the  .approach  presented  by  Altman  in  Reference  20.  These  expressions 


arc  summarized  below; 

F 

s 

^  ^ - f - ] 

\  (h/H)^  -  1./ 

=  0.215  p  L  H  sin  0 

•^w  w 

(0-6) 

F 

e 

=  TP»  ^  \  ‘=p>  « 

(0-7) 

M  : 

■  ^s  ^CP 

(G-a) 

where :  F 

s 


F 


SOD 


F 

e 
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C 

s 


C 

r 


lateral  current  force  at  the  specified  water  depth 
lateral  current  force  in  deep  water 
longitudinal  current  force 
yaw  current  moment 
displaced  volume 

wetted  surface  coefficient,  input  on  SHIP  record 

hull  resistance  coefficient,  input  on  SKIP  record 

or  calculated  as  C  +  C,  +  O.OOOS  (G-9) 

r  f 

residuary  resistance  coefficient  (see  following 
discussion) 


G-2 


=  frictional  resistance  coefficient, 


0.456 


(losio  R,) 


2.58 


1,700 

R 


R  ^  .5  X  10“ 
e 


(G-IO) 


0.002 
C  V  L, 


R  <  5  X  10“ 
c 


W 


Reynolds  number  for  the  hull  (based  on  iongitudinal 
component  of  flow  and  ship  length) 


(G-ll) 


"CP 


90 


distance  troin  niidsliips  to  hull  centcM'  of  pressure 

H  0,0022.6  (0  -  90“ )1  for  0“  S  0  <  180" 

+  0.00226  (6  -  270"))  for  180"  SOS  360" 

ratio  of  distance  to  center  of  pressure  at  0  =  90" 
to  the  distance  to  tlie  center  of  hull  side  area 


(G-12) 


h 

H 


P 


w 


L 


w 


C 

P 


water  depth 
ship  draft 

mass  density  of  water 
waterline  length  of  ship 
propeller  projected  area 
propeller  drag  coefficient 


Several  of  these  terms  re<juire  further  discussion.  The  hull  resistance 
coefficient,  C„,  represents  the  sum  of  various  coef ficient.s  for  different 
sources  of  hull  resistance.  This  coefficient  can  be  input  or  calculated 
in  the  computer  program.  When  no  input  is  given  for  C  ,  it  will  be 
calculated  as  tiie  sum  of  a  residuary  resistance  coefficient,  a  frictional 
resistance  coefficient,  and  a  foul ing/snt face  effect  coefficieuL.  The 
fouling/surface  effect  coefficienl  is  given  an  arbitrary  value  of  0.0005. 

The  frictional  coefficient  is  calculated  from  K<piatioti  G-10  and  the 
residuary  coefficient  is  obtained  from  linear  interpolation  of  a  digitized 
form  of  Figure  38  of  Reference  20.  This  method  of  obtaining  C  is 
limited  to  low  flow  velocities  since  wave-making  resistances  are  ignored. 

The  longitudinal  location  of  the  center  of  pressure  for  a  hull 
skewed  with  re.spect  to  the  flow  is  estimated  by  Equ.uioii  G-12.  This 
requires  an  estimate  of  the  ratio  of  the  distance  to  the  center  of 
pressure  and  the  center  of  area  for  beam  flow,  •  This  factor  is 

estimated  by  linear  interpolation  between  the  values  for  the  ship's 
Dl)-692  and  EC-2  using_the  block  coefficicMit  as  a  reference.  Reference  20 
gives  the  values  for  for’  the  l]l)-692  and  EC-2  as  0.056  and  -0.138, 

respectively.  (Negative  means  aft  of  midships.) 
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It  should  be  emphasized  that  these  analytical  expressions  are  to  be 
viewed  as  a  convenient  alternative  to  the  DM-26  experimental  curve  proce¬ 
dure.  It  remains  to  be  demonstrated  that  they  are  capable  of  giving 
reliable  approximations  of  the  ship's  loading. 
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Appendix  11 

USER-SUPPiJED  SUBROUTINES 


The  SEADVN  program  allows  additional  modeling  flexibility  by  providing 
for  three  user-defined  subroutines.  They  are:  DRAGCO,  TFNUSR,  and 
CURUSR. 

The  DRAGCO  subroutine  provides  for  definition  of  lino  and  lumped 
body  drag  eoef f icicnts .  The  routine  is  called  each  time  the  coefficient 
is  required.  The  necessary  subroutine  parameter  definitions  are: 

SUBROUTINE  URAGCO  (INDRAG,  IBOD,  IDR,  RE,  RET,  CN,  CT) 
where:  INDRAG  =  option  call  number  from  PROB  record 

IBOD  =  1  -  spherical  buoy 

2  -  cylindrical  buoy 

3  -  cable 

IDR  =  drag  coefficient  munber  specified  by  lumped  body 
(BODY  record)  or  cable  (MATE  record)  input  data 

RE  -  Reynolds  number  based  on  the  normal  component  of  the 

relative  fluid  velocity 

RET  =  Reynolds  number  based  on  the  tangential  component  of 
the  relative  fluid  velocity  (not  given  for  spheres) 

CN  ~  return  variable  for  the  calculated  normal  drag 

coefficient 

CT  =  return  variable  for  the  calculated  tangential 

drag  coefficient  (not  used  for  spheres) 

INDRAG,  IBOD,  and  IDK  can  be  used  as  indices  to  select  the  drag  coe  iicients 
from  user-defined  catalogs  of  functions  if  that  is  desired. 

The  TFNUSR  provides  a  single-valued  function  in  time  that  defines 
the  time  variation  of  loads,  currents,  motions,  payout,  etc.  The  necessary 
subroutine  parameter  definitions  are: 

SUBROUTINE  TFNUSR  (T,  F,  N,  TPARM) 

DIMENSION  TPARM  (1) 
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whoro:  T  =  ciirrcnt  time,  T  S  0 

F  -  tlio  returned  value  of  the  time  funet  iou 

TPAKM  ^  a  single-dimensioii^^  iii'iuiy  of  user-defined  iniuit  ])iii'amel er 
piovided  for  the  N  funelioii.  The  inaxiniuin  luiiiilx'r  of 
paramoters  is  20. 

The  CURllSK  suhroutiiu'  is  cal]<-d  to  define  the  fluid  velocity  at  ail 
node  points  in  the  structure  when  the  Fl.OW  lihr€iry  indi<'ates  a  user- 
defined  routine.  Unless  signaled  otherwise,  the  routine  is  called,  at 
every  iteration  or  stop  of  the  subaria  1  ys  i  s  .  DYN  and  TSSS  subar.a  1  yses 
can  c.il  1  CllRl.ISR  to  get  the  spare-dependent  flow  components  and  nsi'  the 
TFUN  library  to  dot  i  no  lime  variation  or  c.'iii  re<|niri'  UURMSR  lo  give  both 
lime  and  space  variation.  (See  the  Fl.OW  library  data.)  The  neress.ary 
subroutine  parameter  definitions  are: 

SUBROUTINE  UURUSR  (T,  N,  NN,  X,  V,  FI, PAR) 

DIMENSION  X(;i,l),  V(;),l),  Fl.PAR(l) 


where : 

■P 

=  current  time,  T  r?  0.0 

N 

=  f.'ow  field  number,  N  >  0 

NN 

=  niimhor  of  nodes 

X  -  nodal  positions  (X,  Y,  /,  position  lor  each  noile) 

V  =  nodal  flow  vector  (X,  Y,  Z  velocity  for  each  node). 
These  values  for  all  NN  nodes  are  Lo  be  returned  at 
each  ca 1 1 . 

FLPAR  =  a  single-dimensioned  array  ol^jiser-detined  input 
parameters  provided  for  the  N  flow  field.  I'lie 
maximum  numher  of  parainelers  is  10. 
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Appendix  I 

OPTIONAL  DEBUG  OUTPUT 


The  SAO  data  allow  for  requests  of  extra  output  data.  In  general, 
this  r  utput  was  intended  for  use  in  debugging  and  checking  the  program. 

It  nas  also  been  foand  useful  in  finding  difficult  input  errors  and  ana- 
;j:’:ing  ..he  efficiency  of  various  choices  of  analysis  methods  and  parameters. 

'  itu55i  i:=  • ’:  u^-ually  requested,  and  little  effort  has  been  expended 

in  gencr-'-  n  identify  the  output.  It  is  assumed  that  when 

the  hiy  ••  cis  ot  output  detail  are  reqiie^tod,  the  user  is  familiar 
r-'Quo*  .vith  ''*^Oi>ram  structure  and  the  analysis  methods  to  be  able  to 
•'  .'ize  the  dc't  irtellige-r ly .  Indiscriminant  requests  for  extra 
.)!  t  i  ut  a.-  .'y  disco. -raged  since  they  greatly  add  to  the  output 

.  ,i.iime  and  increase  run  costs. 

Tte  ^ORTRAN  ■  ■  iabZe  I?-'  •_  wiic  input  used  to  control  the  extra 
outy Generalii,  .ow  levels  of  output  are  obtained  by  small  values 
for  IB(.'  When  IBG  ic  less  than  7  only  a  small  amount  of  specific  output 
is  usually  produced.  Valuer,  of  IBG  greater  than  6  will  give  significant 
amounts  of  output  on  iterative  analyses.  The  following  is  a  list  of  the 
values  of  IBG  recognized  in  che  subanalysis  and  a  brief  description  of 
the  output  they  generate  (frequency  codes:  I  -  each  iteration,  S  -  each 
step,  F  =  each  wave  frequency,  H  =  each  wave  heading): 

IBG  FREQUENCY  EXTRA  OUTPUT 

S  The  number  of  iterations  required  to  get  convergence 

of  iterative  solutions  (DEAD,  LIVE,  DYN). 

A  message  whenever  the  reference  configuration,  R^, 
is  updated. 

I  For  VR  Method:  standard  output  record. 

2  I  A  message  whenever  an  element  is  in.  compre.ssion  with 

no  compression  stiffness  (i.e.,  slack). 

A  message  whenever  a  buoy  or  anchor  changes  its  status 
relative  to  the  surface  and  bottom  limits. 

The  values  of  the  net  vertical  forces  applied  to  buoys 
and  anchors. 

Best  Available  Copy 


I-l 


IBG 

3 


4 

5 


FREQUENCY 

S 


EXTRA  OUTPUT 

The  values  of  nodal  displacement,  velocity  ..nf’ 
acceleration  in  (see  Ref  1);  the  nodal  displace¬ 
ment  and  acceleration  at  the  incremental  reverence 
state;  the  nodal  positions  and  forces  in  t.,  for  each 
node. 

Payout  update  data  (not  implemented  yet) . 

Ship's  static  loads  in  local  and  global  coordinat.'.s . 

For  MNR  Method:  five  values  are  printed  -- 
residual  norm,  displacement  norm  (the  one  used),  and 
the  values  for  the  three  alternative  displacement 
norms . 


1 


) 


F 


I 

1 


For  VR  Method:  current  numerical  damping  data 
displacement  and  velocity  norms  current  control 
parameters . 

For  the  DIM  Dynamic  Method?  six  values  are  printed  -- 
the  iteration  number,  the  displacement  norm,  the 
maximum  acceleration  increment,  and  the  values  for 
three  alternative  displacement  norms. 

For  the  FREQ  Subanalysis:  the  values  of  the  s  f ady- 
state  solution  vector,  and  the  response  amplitude 
operator  data  for  displacements  and  tensions. 

For  the  MNR  Method:  A  message  when  the  tangential 
stiffness  matrix  is  calculated.  A  message  when  the 
ID  search  is  initiated  and  terminated. 

Element  tensions. 

The  values  of  DtlU,  BETA,  and  the  modified  diagonal 
terms  of  the  tangential  stiffness  matrix  when 
numerical  damping  is  used. 

A  message  when  alternating  estimates  are  detected 
and  the  next  "guess"  is  between  them. 

The  current  positions  of  the  nodes. 

For  VR  Method:  current  displacements  and  positions. 

For  the  Dynamic  Analyses:  A  message  whenever 
alternating  estimates  are  detected  and  the  next  "guess" 
is  between  them. 
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EXTRA  OUTPUT 


IBG  FREQUENCY 

S  For  the  RFB  Dynamic  Method:  A  message  when  tangent 

stiffness  matrix  calculated. 

Element  tensions 

Effective  loads 


I 

U  1 


I 


I 


F 

S 


Incremental  displacements 

For  the  FREQ  Subanalysvs:  the  values  of  the  ship's 
roll  angle. 

For  the  HNR  Method;  current  values  of  the  residual 
norm,  the  components  of  the  residual,  the  force  vector, 
the  total  displacement,  and  the  incremental  displacement. 

IBG  is  reset  to  7  after  five  steps. 

For  the  DIM  Dynamic  Method:  current  values  of  the 
force  residual,  acceleration  increment,  and  displace¬ 
ment  increment. 

For  the  FREQ  Subanalysis:  the  values  for  the  wave- 
induced  ship's  forces  (local  coordinates),  the  force 
vector  (global  coordinates),  and  the  steady-state 
response  vector. 

Details  of  dynamic  tension  calculation. 

For  the  RFB,  \*R  and  SLI  Static  Methods:  the 
incremental  cr  effective  forces  and  the  displacement 
increments . 


^9  Once 


^10  H 

I 


For  RFB  Dynamic  and  HNR  Methods:  the  terms  of  the 
tangential  stiffness  matrix  in  compact  storage  format. 

IBG  set  to  8  after  one  printing. 

For  FREQ  Subanalysis:  the  dimensionalized  form  of 
the  ship's  mass  and  restoring  matrices  (Mg  and  Kg) 

Diagonal  terms  of  the  system  matrix,  +  iu»C) . 

The  parameters  used  ir.  interpolating  the  ship  motion 
File*. 


The  ship's  added  mass  and  damping  matrices  (M^^,  Cg) 
and  loading  deta. 
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IBG  FREQUENCY 


EXTRA  OUTPUT 


11  I  The  values  of  fluid  velocity  and  fluid  induced  loads 

at  each  node  in  the  system. 

The  local  and  global  forms  of  the  ship's  static 
loads  (same  as  IBG  =  5) 

The  local  to  global  transformation  matrix  for  ships. 

^11  H  The  mass  matrix  in  compact  form  for  the  mooring 

system  in  a  FREQ  subanalysis. 

^12  Once  The  global  stiffness  matrix  in  compact  form. 

S13  Once  Each  element  stiffness  matrix  and  local-to-global 

transformation  data. 

The  element  contributions  to  the  residual.  IBG  set 
to  11  after  one  printing. 

14  Once  Same  as  IBG  =  13  except  IBG  is  set  to  zero  after 

one  printing. 
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Appendix  J 

PROGRAM  SIZE  LIMITATIONS  AND  STORAGE  REQUIREMENTS 


An  attempt  has  been  made  to  make  the  SEADYN  program  flexible  in  siae. 
There  are,  however,  some  specific  limitations  imposed  by  FORTRAN  dimension 
statements.  These  have  been  chosen  large  enough  to  accommodate  most  models 
and  can  be  increased  by  modifying  the  specific  dimension  statements 
and  size  variables.  This  requires  program  recompilation  and  assembly, 
the  description  of  which  is  beyond  the  scope  of  this  manual. 

The  specific  size  limitations  are: 


Maximum  number  of:  bodies  in  BODY  record  50 

body  locations  given  by  BLOC  record  50 

limit  conditions  in  LIMI  record  50 

limit  locations  given  by  LLOC  record, 

etc.  50 

lines  connecting  to  a  limited  node  10 

cable  materials  in  MATE  record  10 

entries  in  any  tension/strain  table  20 

entries  in  FLUI  record  2* 

catenary  lines  of  nodes  generated  by 

LINE  record  20 

ship/platform  rigid  bodies  defined  by 

SHIP  record  5 

ship/platform  rigid  bodies  in  FREQ  SAO  1* 

payout/ reel. -in  ends  defines  in  DYN  or 

TSSS  SAO  5 

moved  nodes  defined  in  DYN  oi  TSSS  SAO  5 

lines  connecting  to  a  node  where  an 

anchor  holding  power  CHEK  is  made  20 

strum  strings  defined  by  STRUM  record  30 

elenients  in  any  strum  string  20 

flow  fields  defined  by  FLOW  record  10 

parameters  associated  with  any  flov' 

field  10 

time  functions  defined  by  TFUN  record  20 

parameters  associated  with  any  time 

function  20 

load  variation  sets  defined  by 

LOaD/LVAR  records  3 

PROB  +  REST  data  sets  in  any  run  50 

rigid  format  data  sets  in  any  run  1* 


’’'Program  logic  limitation 


wave  headings  on  ship  motion  file  30 
wave  lengths  on  ship  motion  file  30 
roll  angles  on  ship  motion  file  8 
wind  velocities  on  ship  load  file  5 
wind  headings  on  ship  load  file  20 
current  velocities  on  ship  load  file  5 
current  headings  on  ship  load  file  20 


The  size  restrictions  related  to  the  number  of  nodes  and  line  elements 
that  can  be  included  in  any  model  are  not  rigidly  defined  by  dimension 
statements.  A  form  of  variable  dimensions  is  used,  which  takes  a  given 
block  of  storage  (labeled  common  ACOM)  and  partitions  it  according  to 
the  problem  size  and  specific  net;ds  of  each  analysis  option.  The  main 
program  for  SEADYN  is  simply  a  routine  that  defines  the  size  of  common 
and  calls  the  controlling  routine.  The  minimum  storage  required  for  a 
problem  is  given  by  the  variable  NBASE  (printed  with  PROB  data  set 
output) : 


NBASE  =  (31  *  NE)  +  (66  *  NN)  (J-1) 

where;  NE  =  number  uf  elements 
NN  =  number  of  nodes 

It  is  possible  to  reduce  the  size  of  NBASE  for  problems  that  will  not 
use  stor:^ge  for  velocities  and  accelerations.  At  present,  the  logic  to 
do  this  ib  aol  implemented,  although  a  flag  to  request  it  is  provided  in 
the  PROB  data  set. 

Additional  storage  is  required  by  each  of  the  analysis  options. 

The  formulas  used  to  calculate  storage  needed  for  each  SAO  type  are; 

DEAD,  LIVE,  TSSS,  DYN: 


NEED 

=  NBASE  +  NF3*IB 

(J-2) 

MODE: 

NEED 

=  .5*(3*NF3*KF3  +  5'W3) 

(J-3) 

FREQ: 

NEED 

=  NBASE  +  2*NF3'>IB 

(J-4) 

CHFK: 

NEED 

=  NBASE 

(J-5) 

where;  IB 


equation  half  bandwidth  (0  for  DYN  -  DIM  solutions) 


NF3  -  3*(NN  -  NSLAVE) 

NSLAVE  =  niuiiber  of  slave  nodes 
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Storage  size  checks  are  made  at  the  beginning  of  each  SAO  to  deter¬ 
mine  if  enough  space  is  available.  If  not,  a  message  is  printed  to 
indicate  the  space  needed,  and  the  run  is  aborted. 

The  major  users  of  storage  space  are  simultaneous  equation  solvers 
and  the  Jacobi  eigenvalue  solver.  These  analyses  are  done  entirely  in 
core  storage  assuming  sparse  and  banded  matrices  that  are  symmetric. 

This  poses  only  minor  difficulties  on  virtual  memory  machines.  Real 
memory  machines  (even  with  an  extended  core)  can  place  severe  restrictions 
on  problem  size.  Use  of  an  extended  core  (LCM  on  CDC)  requires  special 
modifications  that  go  beyond  the  scope  of  this  manual.  Get  a  programmer's 
help!  Further  relaxation  of  these  restrictions  are  possible  but  require 
significant  program  modifications. 

An  example  of  the  size  limitations  (number  of  elements  and  nodes) 
for  a  representative  real  memory  machine  is  given  below: 

Machine:  CYBER  76  (CDC7600)  SCOPE  2.1  Operating  System 
Maximum  Available  Small  Core  Memory:  I60,000g 

Maximum  Program  Size  (Segmented  form  -  full  program  w/o  ACOM):  126,000g 
Available  for  ACOM:  32,000g  or  13312^^ 

Assume  half  bandwidth  of  10  and  N£  =  NN  ,  then 


for  DEAD,  LIVE,  TSSS; 

(31  +  66)  NN  +  30  NN  =  13312 


for 


for 


NN  = 


DYN  - 
NN  = 


13312  _ 

127 

DIM: 

13312 

97 


MODE: 

27  NN2  +  15  NN  - 


104 


137 


(13312)(2) 


0 


NN 


-15  ±  225  +  (8)(27)(13312) 

54 


31 


for  FREQ: 


(31  +  66)  NN  +  60  NN 


NN 


13312 

157 


84 


13312 
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SEADYN  uses  the  following  files: 


FORTRAN 

File  Name 

Format 

01 

Bina ry 

02 

Binary 

03 

Binary 

04 

Bina ry 

05 

BCD 

06 

BCD 

08 

Binary 

09 

Binary 

10 

Binary 

11 

Binary 

12 

Binary 

13 

Binary 

15 

Bina ry 

16 

BCD 

20 

Binary 

Use 

DEAD  restart  save  file 
LIVE,  TSSS  restart  save  file 
DYN  restart  save  file 
Alternate  restart  input  file 
System  input  file 
System  output  file 
Ship  motion  file 

Scratch  file  for  temporary  data  storage 
Ship  load  data  file 
FREQ  steady-state  response  solutions 
Storage  of  FREQ  response  data  (RAOs) 

Scratch  file  for  FREQ  wave  heading  data 
Deciphered  record  images  of  frce-form  input 
Scratch  file  for  rigid  format  input  data 
Mode  shape  output 

Scratch  file  for  mode  shape  calculations 


21 


Binary 


